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We consider strong-field effects in laboratory and astrophysical plasmas and high intensity laser and cavity sys- 
tems, related to quantum electrodynamical (QED) photon-photon scattering. Current state-of-the-art laser fa- 
cilities are close to reaching energy scales at which laboratory astrophysics will become possible. In such high 
energy density laboratory astrophysical systems, quantum electrodynamics will play a crucial role in the dynam- 
ics of plasmas and indeed the vacuum itself. Developments such as the free electron laser may also give a means 
for exploring remote violent events such as supernovae in a laboratory environment. At the same time, super- 
conducting cavities have steadily increased their quality factors, and quantum non-demolition measurements are 
capable of retrieving information from systems consisting of a few photons. Thus, not only will QED effects such 
as elastic photon-photon scattering be important in laboratory experiments, it may also be directly measurable 
in cavity experiments. Here we describe the implications of collective interactions between photons and photon- 
plasma systems. We give an overview of strong field vacuum effects, as formulated through the Heisenberg-Euler 
Lagrangian. Based on the dispersion relation for a single test photon travelling in a slowly varying background 
electromagnetic field, a set of equations describing the nonlinear propagation of an electromagnetic pulse on a 
radiation-plasma is derived. The stability of the governing equations is discussed, and it is shown using numerical 
methods that electromagnetic pulses may collapse and split into pulse trains, as well as be trapped in a relativistic 
electron hole. Effects, such as the generation of novel electromagnetic modes, introduced by QED in pair plasmas 
is described. Applications to laser-plasma systems and astrophysical environments are discussed. 
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!. INTRODUCTION 

Nonlinear effects, in which a given phenomenon affects its 
own evolution and dynamics, are prominent components in 
a large variety of physical, chemical, and biological systems. 
The examples range from optical and ner ve fibers, au tocat- 
alytic chemical reactions, to ocean waves dScottl 120031) . The 



field of hydrodynamics lias been especially important for the 
development of nonlinear physics, both concerning analyti- 
cal and computational tools, since there the nonlinear effects 
can play a major role in systems with important applications, 
e.g. meteorology. The subject of plasma physics is a natural 
generalization of the field of hydrodynamics, since it builds 
on the fluid or kinetic equations, while adding the electro- 
magnetic interaction. The plasma state of matter is promi- 
nent in large regions of the Universe, such as our closest star, 
the Sun, accretion discs, and even interstellar clouds. It has 
since long also been noted within the field of plasma physics 
that both nonlinear effects and collective interactions can give 
rise to important new physical effects, such as the pondero- 
motive force concept and Landau damping ( Hasegawai 1197,'^ . 
The low-frequency ponderomotive force, which arises due to 
nonlinear couplings between high-frequency electromagnetic 
fields, plays a central role in the physics of laser-plasma inter- 
actions. This force in an unmagnetized plasma is expressed 
as the gradient of the electromagnetic field intensity, which 
pushes electrons locally and thereby creating a hugee space 
charge electric fields and the plasma density cavitaties. Due 
to the radiation ponderomotive force, one has the possibility 
of many interesting nonlinear pheneomena in plasmas, e.g. 
the generation of intense wakefields, stimulated scattering of 
electromagnetic waves off plasmo ns and phonon s, localiza- 
tion of electromagnetic fields, etc. (Eliezer, 2002). The mo- 
mentum and energy transfer from the laser field to the plasma 
parti cles can be har nessed in, e.g. inertial confinement fu- 
sion lElieze?, "2002'). Moreover, the intense electromagnetic 
radiation generated in state-of-the-art lasers can be used to 
model certain astrophysical plas ma conditions in a laboratory 
environment ( Remingtonl l2005l) . Questions of astrophysical 
interest that can be approached within the field of high en- 
ergy density laboratory astrophysics range from the equations 
of s tate of pl anetary interiors to supernova shock formation 
rsee lHEDLAl (12005) for an overview). In the next generation 
laser-plasma systems the influence of quantum electrodynam- 
ics will become important, and fundamental questions related 
to the nonhnearity of t he quantum vacuum can be approached 
in laboratory systems ( lMourouetalll2005l) . 

Currently, lasers are capable of reachin g intensit ies of 
ipzi 1022 w/cm^ llBahk etalll2004|M purou et all Il998i 
|2005'; Tajima and Mourou, 2002; Tajima, 2003). At such high 
field strengths, the quiver velocity of the electrons is highly 
relativistic, and the radiation pressure, manifesting itself as 
a ponderomotive force term in the evolution equations for 
the plasma, gives rise to local electron expulsion. More- 
over, at these intensities, the nonlinear relativistic dynamics 
of the laser-plasma system gives rise to a number of other in- 
teresting pheno mena as well , such as soliton formation and 
pulse collapse ( Sh ukla et all '1986). The latter could be of 
interest when using laser-plasma systems to generate electro- 
magn etic field intensities approaching the Schwinger intensity 
limit (Bingham, 2003; Bingham et al., 2004; Bulanov et al., 
[2003: Cairns et al.„ .2004.: Mourou et al.. .2005.: .Shukla et aU 
I2004ell2005h . 

The event of future ultra-short (in th e femto-second range) 
intense (lO^^ - 10^^ W/cm^) lasers dMourou et all 1 19981 



l2005t iTaiima and Mouroul l2002t iTaiinial 1200 3") could gener- 
ate new physics within the next few years (see Fig. This is 
based on the development of chirped pulse amplification, and 
the evolution of laser power is predicted to continue evolv- 
ing for quite some time ( Mourou e t al., 1998, 2005). The X- 
ray free electron lase rs (XFEL) under construction at SLAC 
dSLAC LCLSL l2005h and DESY (DESY XFEL, 2005) wiU 
be a major source of experimental data not achievabl e with 
today 's systems, ranging from r nolecular pr operties dPatell 
l2004h to astrophysical conditions ('C henLl2003h . such as super- 
nova shocks ( Woolsev et al. , 2004). The intensity at the XFEL 
focus is expected to reach intensities making the quantum vac- 
uum dkectly accessible for observations cRingwald. .200 1 adbl 
l2003[ see Fig. |3}- Moreover, combined effects of laser pulse 
collapse and ponderomotive force electron expulsion would 
be able to create plasma ch annels in which ultra-high intensity 
field strengths are reached dYu et al.LlT98l . such that the non- 
linear vacuum effect of elastic photon-photon scattering could 
become important (.Shen and Ya..2003:.Shen et aL.2003i) . 
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FIG. 1 The evolution of laser intensity (reprinted with permission 
from Ffaiima and Mouroill2002h . 

A majority of studies have not taken into account the in- 
fluence of electron-positron pair creation or elastic photon- 
photon scattering on the dynamics of laser-plasma sys- 
tems (there are however important exceptions, see e.g. 
Ifiulanov et alJ (|2005)). Effects of this kind wiU be of the ut- 
most importance when laser compression schemes approaches 
the critical field strength 
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(1) 



as the nonlinearity of the quantum vacuum becomes pro- 
nounced. Here me is the electron rest mass, c is the speed 
of light in vacuum, e is the magnitude of the electron charge, 
Ae = h/rrieC is the Compton wavelength, and h is the Planck 
constant divided by 27r. Thus, for such extreme plasma sys- 
tems, the concept of photon-photon scattering, both elastic 
and inelastic, has to be taken into account. 
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FIG. 2 The connection between different high-energy regi ons in 
physics and experiments (reprinted with permission from IChenL 



The interaction of high-intensity laser pulses with plas- 
mas has applications to other fields of science, e.g. table-top 
particle accelerators (iBinghaml Eool . Also, achieving field 
strengths capable of producing pair plasmas in the laboratory 
could facilitate a means of producing anti-matter on a more 
or less routine basis, as they currently are at high energy ac- 
celerators. However, it should be emphasized that the produc- 
tion of pairs from intense lasers requires that severe technical 
constraints can be overcome, such as phasing of the two in- 
teracting short electromagnetic pulses. Even if routine pair 
production via laser systems is not to be reached within the 
near future, the related possibility of directly detecting elas- 
tic p hoton-photon scattering is indeed a fascina ting possibil- 
ity ( Brodinetal . , 200 1 ; Solia cic and Segevt !2000b ) . Further- 
more, the creation of multi-dimensional high intensity elec- 
tromagnetic pulses, using guiding structures (suc h as plasma 
boundaries), could result in pulse collapse ( Brod in et all 
12003: Shukla et al.. 2004b) . Such pulse collapse would give 
rise to intensities close to the Sch winger critical field ([l). 
Thus, the combination of laser-plasma interactions and QED 
effects, such as pair production and photon-photon scattering 
could spark new methods for producing conditions reminis- 
cent of astrophysical environments in future experiments (see 
Fig. 12). In fact, some of the most pertinent in today's funda- 
mental physics research , such as the question of dark matter 
(e.g. through the effects of light pseudoscalar fields, such as 
the axion field, o n QED interac t ions and light propagation, se e 
[Bernard' ("1999); Bradley et al. ('2003); Dupays et al. (2005)), 
cosmic accelerators (s uch as t hroug h laboratory plasma wake- 
field accelerator tests llChenIl2003h) , and possible new high- 
density states of matter ("Remington", "20051, are related to the 
high energy events for which laboratory astrophysics would 
yield valuable insight. Thus, it is of interest to study such 
high energy scenarios, e.g. photon-photons scattering in the 
context of laser-plasma systems (Bulanov, 2004), as these are, 
in the near future, likely to yield the right conditions for such 
events to take place. 



A. Nonlinear quantum electrodynamics 

In classical electrodynamics, as described by Maxwell's 
equations, photons are indifferent to each other as long as 
there is no material medium present. This is not so in quan- 
tum electrodynamics (QED). Due to the interaction of pho- 
tons with virtual electron-positron pairs, QE D offers the pos- 
sibility of photon- photon scattering (iHeisenberg and Euleil 
[l936; Schwinge^ Il95ll) . This is commonly expressed 
through the effective field theory appro ach represented 
by the Heisenberg-Euler Lagrangian (Fradkin et al.|, Il9^ 
i Greiner et all I1985L iGrib et all Il994t [Se isenberg and EukJ 
Il936t ISchwinged Il95lt IWeisskopfi fm&i . neglecting dis- 
persive effects. T his Lagrangian [see Eg. dp below] and 
its generaUz ations dDittrich and Giesl l2O()0tlDunn3 l2004t 
IVallurietalll2003h . gives nonlinear corrections to Maxwell's 
vacuum equations, similar to the self-interaction terms en- 
countered injionfinear_oPtics_du^^ presence of a Kerr 
medium dAgrawall 1200 ll iBloembergenl 11996.) . and higher 
order corrections can easily be incorporated by taking into 
account higher vertex order diagrams. Since the lowest 
order effective self-interaction term is proportional to the 
fine structure constant squared, the field strengths need to 
reach appreciable values until these effects becomes impor- 
tant, see Eg. (|l) (Fradkin et al., 1991; Greiner et al., 19851 
iGrib et all The corrections give rise to both single 

partic le effects, such as closed photon paths ("Novello et aU 
I2OOII) . vacuum birefringence ( Hevl and Hernquist, 19972^, 
photon splitting (lAdleil I1971I) and lensing effect i n strong 
magnetic fields (see, e.g. iDe Lorenci et al.l ( 120001) : iHardind 
([199 1)), as well as collective e f fects, li ke the self-focusing 
of beams (Soliacic and Segev, '2000b') or the formation 
of light bullets (Brodin et al... .2003) . Recently, it has 
also been shown, using analytical means, that these ef- 
fects give rise to collapsing structure in radiation gases 
(jMarklund et al. , 2003), results that ha ve been extended and 
confirmed by numerical simulations dShukla and EfiassoiiL 
l2004aft . Possible detection techniques, as well as phys- 
ical implications, of the effects of photon-photon scat- 
tering have attract ed a rather constant interest since first 
discus s ed (e.g. iBialvn icka-Birula and B ialvnicki-Biruk 
(|1970 |) : iDicus et al.l (Il998l); iDing a nd Kapla nI (Il989h: lErbei 
('l966');'Fradkin et al.'('l99l"); Greiner et al. ( 1985); Grib et al 
(.1994): Kaplan and Ding (2000); Latorre et al. (1995) ; Tsai 
(ll974albl)T and the concept of self-trapping of photons due 
to vacuum nonlinearities was discussed independently by 
Rozanov (1993, 1998) and Soliacic and Segev (2000b) in the 
context of the nonlinear Schrodinger equation. 

The above studies assume that the dispersive/diffractive ef- 
fect of vacuum polarization is negligible, and this, of course, 
puts constrai nts on the allowed space and time variations 
of the fields (ISoliacic and Segevt |^00^. In the context of 
pair creation, rapidly varying fields have been analyzed, since 
when individual photons pass the pair creation energy thresh- 
old 2meC^, real electron-positron pairs may be created from 
the vacuum by a "down-conversion" process of photons. Sim- 
ilar processes are thought to be of importance in the neigh- 
borhood of strongly magnetized stars, where the magnetic 
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field induces photon splitting ('Adler et al.Vl970VAdler' 


'l97l', 


Adler and Shubert!, 'l996; .Baring and Harding, 2001; 


Erben 


ip66), and may effectively absorb the photons ( Duncan 
Hevl and Hernauist, il997al). It has been suggested t 


I2002t 


lat the 



nontrivial refractive index due to photon-photon scattering 
could induce a lensing eff ect in the neighbourhood of a mag- 
netar dShaviv et allll999ft . 



The physics of elastic photon-photon scattering has inter- 
ested researchers for a long time, and several suggestions 
for ways to detect this scattering in the laboratory have been 
made during the last decades ( Alexandrov et al... 1985; Dewari: . 

and the recent strong increase in available laser intensi- 
ties have stimulated various schemes ("Ding and Kaplan',' '19921 
[Rozanov, 1993, 1998). It has been suggested by Brodin et aL 
{2001 , 2002) that the effect of photon-photon scattering could 
be detected using fields significantly weaker (10 MV/m) than 
state-of-the-art laser fields. 

Next, we present some physical systems in which the ef- 
fects of photon-photon interactions may either be of impor- 
tance (magnetars), or become important in the near future 
(such as state-of-the-art laser-plasma systems). 



1. Intense field generation 

a. Electromagnetic cavities High performance, i.e. large elec- 
tromagnetic fields combined with low dissipative losses, can 
be found in superconducting cavities, which among other 
things are used for particle acceleration (Graber, 1993|). The 
waves that can be sustained within such a cavity can have a 
field-strength i? ~ 10 MV/m, i.e. close to the maximum that 
can be tolerated by the walls without field emissions. For such 
cavities, the different high intensity wave modes can act as 
pump waves for the quantum vacuum. Through the interaction 
between these waves and virtual electron-positron pairs, new 
modes with well-defined frequencies and wavenumbers will 
be generated. Those satisfying the dispersion criteria for the 
given cavity could then also reflect within the cavity with very 
small losses, thus yielding a method for detection of the quan- 
tum vacuum nonlinearities. For example, for a cavity resis- 
tance R ^ 1 nil, corresponding to superconducting niobium 
at a temperature 1.4 K and a frequency uj ^ 2x 10^" rad/s of 
the wave mode generated via the nonlinear quantum vacuum, 
one finds that the saturated energy flux P3 of the generated 
mode is of the order of 10~^ W/m^ (see Sees. | III.A.4l and 
lIV.A.St ( Brodi n et all l200lt lEriksson et all l2004 . This en- 
ergy flux is above the detection level by several orders of mag- 
nitude. However, one should note the importance of the super- 
conducting walls for the output level of the excited mode. For 
copper at room temperature, the cavity resistance increases by 
a factor ^ 10^ as compared to the above example, and conse- 
quently the energy flux of the excited mode falls by a factor 
~ jjjg latter case, it is questionable whether the 

excited signal can be detected. The concepts of cavity mode 
interactions and cavity experiments will be further discussed 
in Secs. lIII.A.4l and lIV.A.5l 



b. Laser development The event of ultra-short (in the 
femto-second ra nge) in tense ( 10^^ — 10^"^ W/cm^ ) lasers 
llMourou et al 1 ll99Ri l2005t iTaiima andMourouI l2002t 
'Taiima','2003') holds the promise of generating large amounts 
of new physics within the next few years. This promise is 
based on the development of chirped pulse amplification, 
and the increase of laser power is predicted to continue for 
quite some time (Mourou et al., 1998, 2005). There are 
two ways for reaching high intensities within laser systems. 
The method most common is to shorten the pulse duration 
(< 100 fs), while keeping the energy content in each pulse 
rather modest (~ 1 — 10 J). Such pulse generation techniques 
can have high repetition rates, which can be advantageous in 
certain experiments where a large number of shots is needed. 
The other route is to increase the pulse contents while keeping 
the pulse duration of the order 0.5 — 1 ps. Such systems has 
the advantage of providing a high signal to noise ratio for 
some experiments. The Nova Petawatt laser at the Lawrence 
Livermore National Laboratory,USA, used this principle, 
and each pulse, which had a duration of ~ 500 fs, had an 
energy contents of ^ 500 J. Similar systems are operating 
at ILE/Osaka, Japan, (ILE/Osaka, 2005) and the Rutherford 
Appleton Laboratory, U.K., JCCLRQ 12005). The OMEGA 
EP laser under construction at the University of Rochester, 
USA, will also work according to the high energy principle, 
and have pulse energies 1 — 2.6 kJ with durations 1 — lOps 

(Iomega 

EP, 2005) . Apart from being a tool for practical 
use, such as inertial confinement fusion and material science, 
intense laser facilities are now of international interest for 
basic research (such as at, e.g. the National Ignition Facility 
at the Lawrence Livermore National Laboratory (USA), 
the Laboratory for Laser Energetics at the University of 
Rochester (USA), the Advanced Photon Research Center 
(Japan), the Institute for Laser Engineering at Osaka Univer- 
sity (Japan), LULI Laboratoire pour 1' Utilisation des Lasers 
Intenses (France), LIL/Laser Megajoule at CEA (France), or 
the Central Laser Facility, Rutherford Appleton Laboratory 
(UK)), and the increased laser output power also gives 
the opportunity, for the first time, to obtain astrophysical 
energy scales in a controlled laboratory setting llChenll2003t 
HEDLA, 2005; Remington, 2005). 

The generation of high (electromagnetic) field strengths is 
at the heart of understanding a variety of phenomena, such as 
astrophysical shocks and jets, in a laboratory setting, and fur- 
thermore forms the basis for a number of applications, e.g., 
table-top plasma accelerators. Thus, schemes and mecha- 
nisms for generating such high fields, other than the direct 
laser pumping, will be of importance to the development of a 
wide range of scientific areas. 



c. The free electron laser The x-ray free electron laser (XFEL) 
is an alternative to the current laser generation techniques, 
and it has as its base the particle accelerator, where high en- 
ergy electrons are generated to obtain high-frequency radia- 
tion (DESY XFEL, 2005; SLAC LCLS, 2005). Within these 
lasers, a large number of coherent photons are generated (10 
orders of magnitude more than regular synchrotron sources). 
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The XFEL concept has a wide variety of interesting appli- 
cations, among these the possibiHty to probe the structure of 
large mole cules, com monly found within molecular biology 
systems (Patel, 2004). It is also hoped that X FEL could form 
the basis of elec t ron-positron pair creation jAUdiofer et all 
l200lt iRingwaii l2001al [Roberts et all l2002l) . Both at the 
TESLA collider at DES Y and LCLS at SLAC, the energy den- 
sity at the focus (with a spatial width ^ 10^^° m and on the 
time-scale 10~^^ s) of t he XFEL is expect ed to reach energy 
densities of lO^^ J/m^ ( iRingwaldl 12001 ah . see Fig. This 
corresponds to electric field strengths of the order 10^*' V/m, 
i.e. two orders of magnitude above the Schwinger critical 
field Q, at which pair creation is expected to take place. The 
possibility to 'fuel' the generation of electron-positron pairs 
by nonlinear effects is therefore a very promising prospect, 
and as some authors have noted , this essentia lly amounts to 
'boiling the vacuum' llRingwaldL EoOlbl l2003h . " Moreover, it 
is believed that nonlinear QED effects will be a very impor- 
tant component of the interaction of XFEL generated radiation 
with dense media. Therefore, it is of interest to achieve an un- 
derstanding of the influence of those nonlinear effects within 
the parameter regime obtainable by the XFEL. 
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FIG. 3 Evolution of peak brilliance, in units of pho- 
tons/(smrad^ mm^ 0.1 %xbandwidth), of X-ray sources. Here, 
ESRF stands for the European Synchrotron Radiation Facility in 
Grenoble. (Source: DESY Hamburg, http : / /xf el . desy . de/) 



2. Laser-plasma systems 

High laser powers generate enormous radiation pressures, 
and accelerate particles to relativistic velocities. The fields 
generated by the particles will, therefore, backscatter and cre- 
ate a nonlinear feedback, something which can give rise to, 
e.g. laser pulse compression and electron density cavitation. 



This may play a significa nt role in differ ent propos ed schemes 
of la ser self-fo cusing (iBulanov et all 12003.: Shukla et all 
12005b [ShorokhovetaU l2003l) . completely changing the dy- 
namics of the suggested methods and altering the final results 
in nontrivial ways. This serves as an important example of 
nonlinear effects in the evolution of laser-plasma systems. 

The interaction of high power lasers, reaching intensities of 
10^" — 10^^ W/cm^, with plasmas has long been the back- 
bone in differen t schemes for table-top particle accelerators 
(iBinghamlEool . and is also essential for the concept of iner- 
tial confinement fusion (pliezer, 2002). The ponderomotive 
force generated in high power laser-plasma systems due to 
laser intensity gradients may give rise to a plethora of phe- 
nomena, such as laser pulse self-focusing and filamentation, 
soliton formati on, parametric instabilities, and magnetic field 
generation (see lEliezeJ ( l2002l) and references therein). 

Laser-plasma systems have been suggested as sources of 
high intensity radiation. Due to the laser ponderomotive force 
plasma electrons will be pushed out of the path of the laser 
pulse, trapping and compressing the laser pulse, such that fur- 
ther electrons are pushed out. The plasma can sustain high 
field strengths, and the pulse compression c an therefore reach 
appreciable intensity values dMourou et al In fact, us- 

ing a Langmuir wave as a plasma 'mirror' for the laser pulse, 
pulse intensities c ould reach, and even s urpass, the Schwinge r 
critical field Q (Bu lanov et all l2003t iMourou et all l2005l) . 
Thus, the interaction between intense lasers and plasmas is an 
intrig uing tool for generati ng pulse intensities above the laser 
limit dMourou et allll998h . 



3. Astrophysical and cosmological environments 

Astrophysical environments and events display truly enor- 
mous energy releases. Supernova explosions, black hole ac- 
cretion, magnetar and pulsar systems are a few examples 
of such extreme situations. Moreover, the energy scales 
in the early universe are equally immense, or even greater, 
and our understanding of the origin of the universe is ham- 
pered by the fact that the energy density scales are so far 
from anything that can be generated in a laboratory, except 
perhaps in relativistic heavy ion collisions (RHIC, 2005). 
It is therefore not surprising that these environments can 
often act as laboratories for phenomena that we currently 
do not have technology to reproduce in earth based labo- 
ratories. Quantum electrodynamical nonlinear vacuum ef- 
fects have received a fair amount of attention in strongly 
magne t ized s ystems, such as pulsars (Beskin et al., 1993; 
Curtis, 1982) and magnetar environments (Har dingt 1 1 99 1 : 
Kouveliotou. J_998). The magnetic field strengths of mag- 
netars can reach energy levels comparable, or even surpass- 
ing, the energy corresponding to the Schwinger critical field 
strength 10^^ V/m, thus making the vacuum truly nonlinear 
(nonlinear QED effects in the magnetized vacuum are de- 
scribed in a large number of publications, and for a representa- 
tive but incomplete list, see Adler et al. ( 1970); Adler ( 1971); 
Adler and Shubert ( 1996); Baier et al.l996; Chistyakov et al] 
(.1998.) ; .Erber. (.1966.) ; .Mentzel et al.. a994.) ; .TsaL (1974^ 
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and references therein). The effect of a nonhnear vacuum may 
even be of crucial importance for our understanding of these 
objects (Baring and Harding, 2001). 



II. EFFECTIVE FIELD THEORY OF PHOTON-PHOTON 
SCATTERING 



mal vacuum, in contrast to the non-thermal one, two-loop cor- 
rections will dominate over the one-loop effects ( Gies..200d) . 
However, treating all these studies (we have by no means ex- 
hausted the list of papers in this short expose) in detail is out- 
side the scope of the present paper, and since we moreover are 
interested in the problem of collective effects, the treatment of 
thermal effects, although of interest, is left for a future review. 



The development of quantum electrodynamics was the re- 
sult of a long and collective effort, and paved the way for 
an understanding of the weak and strong forces as well. It 
has, since its advent, been confirmed to an unprecedented ac- 
curacy, compared to any physical theory. It solved some of 
the long-standing conceptual problems of relativistic quantum 
theory, as proposed by Dirac and others, and it furthermore 
changed the way we look at the elementary interactions be- 
tween particles and fields. The theoretical proposal, due to 
Dirac, of anti-matter as a result of relativistic quantum theory 
was put on a firm foundation with QED. 

The quantization of the vacuum has lead to some remark- 
able insights and discoveries. Consider for example the 
Casimir effect l iCasimir and Poldeilll94a lc asim m 
which two parallel conducting plates with area A separated 
by a distance d. Due to the different boundary conditions 
between and outside the plates, there will be a net attractive 
force F cx A/d* between the plates. In a heuristic sense, 
the vacuum between the plates is 'emptier' than outside, since 
fewer states are allowed due the finite distance between the 
plates. R elated to this is the much debated Scharnhorst effect 
iBarton, '1990"; "Barton and S charnhorstL Il993t IScharnhorst', 
[l990, 1998), at which the phase (and group) speed exceeds the 
speed of light c in vacuum. As will be demonstrated later, the 
opposite occurs in the electromagnetic vacuum, i.e. the phase 
and group velocities decrease due to the electromagnetic in- 
fluence on the quantum vacuum. 

In conjunction with any description of photon-photon scat- 
tering, it should also be mentioned the large amount of lit- 
erature and interest in finite temperature effective field the- 
ory effects. Thermal effects generalize the classical results 
of Schwinger in t he weak field lim it jHeisenberg and Eule^. 
L1936V Schwinger, 1951; Weisskopf, 1936). It was pioneered 
by Di ttrich. tl979) who investigated the thermal effects in 
combination with an external magnetic field, and later a com- 
prehensive study using the real time formalism in the case of 
general electromagnetic field ba ckground was performed 



an general electromagnetic Held ba cKgrouna was periorme d 
by lElmfors and SkagerstamI ( 1 19951) (see also iGiesI ( Il999al) ). 
The dispersion relation, including dispersive effects, was dis- 
cussed by Gies. C 1999b.) . and it was later shown that in a ther- 



' Casimii considered particle-pailicle and particle-plate 

iCasimrr and Poldel ll94Sl). and plate-plate systems (Casirnil, 194^ . 
since the problem stemmed from research on colloidal solutions, but 
is most clearly represented by the parallel plate example. The Casimir 
effect has since bee n confirmed by man y different experime n ts (see, e.g., 
iBordag et alj <200ll):lBressi et alj l2003):lHarber etail<200; ): lLamoreauJ 
<1998l) : iMostepanenko and Trunovl Il997l) : ISukenik et alJ il993l) and 
references therein). 



A. The concept of elastic scattering among photons 

Photon-photon scattering is a non-classical effect arising 
in quantum electrodynamics (QED) due to virtual electron- 
positron pairs in vacuum, see Fig. |3 In the low energy 
limit, i.e. hui <C rrieC^ the magnitude of photon-photon 
can be descibed in term s of the differential cross-section 
jBerestetskii et allll982i) 



dn 



324007r^ 



(3- 



(2) 



where u is the photon frequency in the center-of-mass sys- 
tem, is the classical electron radius, and a = e^/Aireohc ~ 
1/137 is the fine structure constant. Integrating Q gives the 
total cross-section 



a. 



_ 973aV2 / huj 
'^'^ ~ 10125 



0.7 X 10"^^ 



fku 
TeV 



(3) 

We note that the cross-section decreases very fast with de- 
creasing photon energy. In the high energy limit, the cross- 
section on the other hand goes like w^^. The cross-section 



reaches a maximum of a-, 



2 X 10 cm "' for photon en- 



ergies hu! ~ rrieC^ jBerestetskii et all Il982h . indeed a very 
small number 

Instead of a microscopic description, the interactions of 
photons may be described by an effective field theory. For- 
mulated in terms of such an effective field theory, using the 
Heisenberg-Euler Lagrangian (valid i n the long wavelength 
and weak field l imit, see Eq. @) ( iHeisenberg and Euleil 
11936; SchwingeJ (19511) . this results in nonlinear corrections 
to Maxwell's vacuum equations, which to lowest order in the 
fine structure constant are cubic in the electromagnetic (EM) 
field. These correction takes the same form as in nonlinear op- 
tics, where the material properties of, e.g. optical fibres, gives 
rise to cubic nonlinea r terms in Maxwe l l's equations, so called 
Kerr nonlinearities jAgrawall l200lt iKivshar and Agrawall 
I2003h . Since the effective self-interaction term is proportional 
to the fine structure constant squared, this means that the field 
strengths under most circumstances need to reach values close 
to the critical field ^ until these effects becomes important 
(Uradkinetal., 1991; Greiner et al., 1985; Gribetal., 1994). 
With this at hand, we now continue to focus on the concept of 
photon-photon scattering. 
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FIG. 4 The Feynman box diagram for the lowest order photon- 
photon scattering at the one-loop level. 



B. Weak field limit 

We first derive the general dispersion relation in the low 
photon energy and weak field limit, assuming (see Eq. Q) 



w < We = mf,(? jh and |E| < 



(4) 



where Scrit is given by Q, and cjg w 8 x 10^"rad/s is 
the Compton frequency. When these constraints are valid, 
pair-creation, both due to single photons and collective ef- 
fects, will be unimportant, and the effective Lagrangian may 
therefore be treated solely in terms of its real part. It should 
be remembered that the second of these constr aints comes 
from the pair creation probability of Schwinger dSchwing eil 
[1951), which was derived for a pure electric field and may 
therefore not strictly be applied to the case of a radiation 
gas. Thus, this investigation goes far beyond the descrip- 
tion of the thermodynamics of nonlinearly interacting inco- 
herent photons with a photon gas in a plasma environment 
ffsintsadze and Mendonca, 1998). 

Photon-photon scattering is a second order effect (in 
terms of the fine structure constant a), and for constant 
or weakly varying fields it can be formulated in stan- 
dard notation using the Euler-Heisenberg Lagra ngian density 
(iHeisenberg and EuleilfT936tlSchwingeiill95lh 



where 



45m4c5 gOTreo-B^i, 3xl029j/m3 



(5) 



(6) 



Moreover, the field invariants are defined in terms of the field 
tensor Fab according to 

^ = iF.bF'^'' = i(c2B2-E2), <S = i^^ab-F"'' = -cE-B, 

pah ^ ^o-bcdp^^^j2^ and ^nd are the lowest order QED 
corrections. We note that ^ = = in the limit of parallel 
propagating waves. The latter terms in (|5jl represent the effects 
of vacuum polarization and magnetization, and the QED cor- 
rected Maxwell's vacuum equations take the classical form, 
using 



D = eoE + P 



M 



(8) 



where P and M are of third order in the field amplitudes E 
and B (see Eqs. ( II 1> and ( I12> below), and /io = l/c^eo. Fur- 
thermore, they contain terms ^ and such that P = M = 



in the limit of parallel propagating waves. It is therefore nec- 
essary to use nonparallel waves in order to obtain an effect 
from these QED corrections. 

From the constituent relations (|8|l we can deduce the gen- 
eral wave equations for E and B according to 



1 



-V^E : 



-Mo 



c2V(V-P) 



|(V X M) 



(9) 



and 



1 a^B 



V^B 



Mo 



V X (V X M) + —(V X P) 



(10) 

Furthermore, the effective polarization and magnetization, ap- 
pearing in (|8|i, in vacuum due to photon-photon scattering in- 
duced by the exchange of virtual electron-positron pairs can 
L agrangia r 

e.g., Soliac ic and SegevI (l2000bl) ') 



be obtaine d from the Lagrangia n (|5jl and are given by (see 
: ic and SegevI (l2000bl) ') 

2Ke^ [2(£;2-c2b2)E + 7c2(E-B)B] , (11) 



and 

M = 2Kelc^ {-2{E^ - c^B^)B + 7(E • B)E] . (12) 

Equations (I9l-( I12> offer the starting point for the study of 
a weakly nonlinear electromagnetic vacuum in terms of the 
classical field strength vectors E and B. 

The corrections in the Lagrangian (|5jl is the power series 
expansion in the field strengths of the ful l one-loop coii'ec- 
tion given bv (iHeisenberg anid Euleil[l936HSchwingeilfT95lt 
IWeisskonfifl936i) 







n ab 








L cnt 




where 





dz 



cot 



(13) 



E„ 



z2 (a2 - 62) 



(jr2 +^2)1/2 



1/2 



^2 



(J? 



cnt 



?2U/2 



1/2 



(14) 

Thus, ^ = (a^— 62)/2and \'^\ = ah. The Lagrangian correc- 
tion ( I13> is the starting point of the effective field theory anal- 
ysis of a strongly nonlinear quantum vacuu m, such as used in 
some studies of photon splitting (see Dit trich and GiesI ( Il995 
l200(]i) and references therein), and defining some of the prop- 
erties o f a strongly nonlinear gas of photons (Marklund et all 
I2005bh . 



C. The dispersion function 

One may find the dispersion relation of photons in 
an arbitrary constant, or weakly varying, electromagnetic 
background ( Bialvnicka-Birula and Bialvnicki-Birula, 197(J 
iDe Lorenci et all 120001: iThomat |2000). Starting from the 
Lagrangian (|5|l, introducing the four-potential A'^ such that 
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Fab = daAb — dbAa, the Maxwell equations resulting from 
the variation with respect to the four-potential becomes 



(15) 



where we adopt the convention (—1,1,1,1) for the metric rjab, 
used for raising and lowering four-indices a, 5, . . . = 0, 1, 2, 3. 
Next, assuming that Fab = fab + 4>ab, where fab denotes the 
varying background field, and (pab fab) is a weak field 
propagating on this background, we find that the background 
satisfies daf"'^ — 0, and daF"'^ = is identically satisfied 
due to the definition of Fab in terms of the four-potential . 

Linearizing Eq. J15> with respect to 0, and Fourier de- 
composing perturbations according to (i>ab{x) — {ka€b — 
ki,ea) exp{ik ■ x) + c.c, where ea is the polarisation vector, 
k ■ X = kax"' and c.c. denotes the complex conjugate, we 
obtain the following algebraic set of equations for the polar- 
ization vector 

Mab'^'' = [k'^gab - kah - keq (SaQfl;, + Maaab) 1 e'' = 0, 



D. Corrections due to rapidly varying fields 

As we saw in the previous section, it is possible to derive a 
dispersion relation for photons moving in a given background 
field. However, the field variations were neglected, and in 
order to take them into account, a modified weak field La- 
grangian must be used. 

It is well-known that the weak field theory of photon- 
photon scattering can be formulated using the effective La- 
grangian density = + J^ue + -^d, where Jfo is the 
classical free field Lagrangian J^he is the Heisenberg-Euler 
correction as giv en in Eq . The derivative corrections are 



I tq 

given by iMamaev et all[l981i) 

= aeo [{daF''''){d,F%) - FabOF'^'] , (21) 

where □ = ^a^", and cr = {2/15)ac^ /ojI ^ 1.4 x lO^^s j^2^ 
As we have seen in the dispersion relation ( I19I I. the param- 
eter K gives the nonlinear coupling. Here, we find that the 
parameter a gives the dispersive effects in the polarized vac- 
uum. Physically, setting the parameter a ^ corresponds to 
to taking correction due to rapidly varying perturbation into 
account. Since the Compton frequency is 10'^°rad/s, we 
see that w/we <C 1 in most applications. Thus, the dispersive 
term is normally a small correction. 

In the previous section, the requirement was assumed 
to be satisfied. Here, even though we include effects of the 
rapidly varying fields, we require that there is no electron- 
positron pair creation, neither by single photons nor by collec- 
tive effects, i.e. the conditions should stiU hold. Further- 
more, the dispersive/diffractive effects must be small, other- 
wise the limit of weak fields would imply unphysical branches 
in the dispersion relation (Rozanov, 1998). 

As in the previous section, we set Fab — daAi, — dbAa, 
and obtain the field equations from the Euler-Lagrange equa- 
tions db\d^/dFab\ = 0. Thus, we have feozanovl 1 1991 
LShuklaetalll20Q4d) 



(16) 

where ab = fbck'^ and cib = fbck'^ has the properties k^ab — 
k''abO. From this it follows that A/"f,fc'' = 0, and that the po- 
larization may therefore, without loss of generality, be taken 
to obey ef,fc'' = 0, corresponding to the Lorentz gauge. In or- 
der to simplify the analysis, it is assumed that the background 
is slowly varying in spacetime compared to the perturbation. 
Using Gb and cib as the polarization eigenvectors gives two 
equations 

k^ = Sneofab^'^k^'k,, and k^ = Uneofabf^'k^k,, (17) 

from Eq. ( I16> for the polarization a° and a°, respectively. 
Since we can decompose fab = UaEb — UbEa + ^abcB^ for 
an observer with four-velocity w", we have 

= -(k-E)2-c2(k-B)2+cj2E2+c2k2B2-2cwk-(ExB) l{l+2(jU)daF 



ah 



(18) 

and w a?. Equations ([T3 can b e written in the clear 
and compact form CBialvnicka-Birula and Bialvnicki-Birulal 
IT97Q) 



e^nda 



{F,aF''')F-'' + UF,dF^^)F" 



c^«c|k|(l-iA|Qn, 



(19) 



where A is 8k or 14k, respectively, depending on the polariza- 
tion state of the photon, and 

IQp = eo|k X E + ck X (k X B)p (20) 
= eo [e^ + c^B^ - (k • E)2 - c2(k • Bf - 2ck • (E x B 



Here, the hat denotes the unit vector The expression ( I19> 
is valid for arbitrary, slowly varying background fields. It is 
straightforward to show that expression ( I20t vanishes in the 
case of a self-interacting plane wave field. The two different 
possibl e polarization directions can be given in a similar man- 
ner (see lBialvnicka-Birula and Bialvnicki-BirulaCl970.) ). 



(22) 

Equation (I22> describes the nonlinear evolution of the electro- 
magnetic field through the nonlinear dispersive vacuum. We 
note that when ct, k 0, we obtain the classical Maxwell's 
equations, as we should. 

Repeating the procedure leading up to Eq. il6\ . we find the 
corresponding expression in the case of a dispersive vacuum 

Mabe'' - [A-C - 2fTfc2 {k'^gab - kakb)] - 0, (23) 

where Af^j, is given by Eq. ( I16> . With and as the princi- 
pal polarization direction, Eq. (I23> yields 



and 



(1 - 2(7fc2)fc2 ^ 8Keo fabf'^k'k,, (24a) 



(1 - 2afc2)fc2 ^ 14Keo/ah/"'fc'fce, (24b) 



for the two different polarizations Of, and cib, respectively. 
When (7 = 0, Eqs. i24l reduces to Eqs. ( I17> . which yield the 
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dispersion relation ( I19> . With cr 7^ 0, we may use Eq. ( I19> 
in the dispersive term of Eqs. J24> . Thus, th e final dispersio n 
relation is of the form |Rozanovl|l998; Shu kla et alll2Q04d) 



1_ iA|Q|2(i + 2aA|Q|2|kp) 



(25) 



where | Q P is given by Eq. ilOi . Thus, we see that the effect of 
the dispersive parameter is, as expected, to make uj a nonlinear 
function of k. 



E. Special cases of weak field dispersion 

1 . Magnetized background 

In the case of a backgrou nd magnetic field Bp, the disper- 
sion relation (0^ b ecome s jAdler et all 1 1970t lAdleii 11971 



l AdlerandShuberl Il996t iDittrich and Giesi Il998t Wbei 
Il966h 



c|k| (1 - iAeoc^lBopsin^ 



(26) 



where 9 is the angle between the background magnetic field 
Bq and the wavevector k. Thus, photons propagating par- 
allel to the background magnetic field will not experience 
any refractive effects, while a maximum refraction is ob- 
tained for perpendicular propagation. This dispersion rela- 
tion will be relevant for photon propagation in pulsar mag- 
netospheres and in magnetar environments, where for exam- 
ple the vacuum becomes birefringent ( Hevl and HernauisJ 
Il997at iTsai and Erbeit Il975h . This in turn may affect 
the optical depth of neutron star thermal emission, and 
thereby also influence the interpretation of pulsar observations 
jHevl and Shaviv, '2 0021 iHevl et all l2003t iLodenaual et"all 



1 974; .Ventura. 197S 

For a magnetar with surface field strength |Bo| = 
10" T (lKouvehotoulll998h . the field has the energy density 
eoc^lBop « 8 X 10"j/m3. Since A - 10k w 1/(3 x 
10^8 J/m^) (see Eqs. (|6} and ([19}), we find from ^ that 
the phase velocity v — w/|k| satisfies w/c « 1 — 0.13 sin^ 6. 
However, the magnetar field strength does not qualify as a 
weak field, and one may question if it is appropriate to use 
Eq. (I26> in this case (see Sec. III. EH . 



2. Random photons in a magnetic field 

From the previous example, we saw that the effect of a mag- 
netic field on the thermal or random distribution of photons 
could be observationally important. Thus, for a random en- 
semble of photons in a strong magnetic field we have a direc- 
tion independent dispersion relation 



c|k|(l-iAeoc2|Bo 



(27) 



We note that the value of the effective action charge A still 
depends on the polarization of the thermal photons. As the 
crust of magnetars are subject to enormous str esses due to 
the immense field strengths (^ 10^*^ — 10^^ T (iKouveliotoui 



ll998|)), it will suffer from crust quakes, at whi ch bursts of low - 
frequency random photons are released (.Kondratvevl l2002ft . 
In such a scenario, the above dispersion relation may be of 
relevance. 



3. Random photons in a plane wave field 

Analogously, we may treat the case of incoherent pho- 
tons on an intense plane wave background Ep. Then, in the 
equilibrium state of the radiation gas, the propagation direc- 
tions of the photons in the gas are random and the EM pulse 
is a superposition of uni-directional plane waves such that 
Bpkp X Ep/c. Thus, we obtain the direction independent dis- 
persion relation 



c|k|(l-|Aeo|E,p). 



(28) 



Equation (I28> is the proper dispersion relation to use in some 
laser-plasma interaction applications, where a large num- 
ber of incoherent photons are produce d fBi neham. .20031 
iBingham et all 12004 ICalrns et allEool . 



4. Radiation gas background 

Consider a single photon transversing a dense background 
radiation gas with energy density S". Then the disper- 
sion relation can be written as (.Dittrich and Giesl l2000t 
iMarklund et'alll2003h 



:|k| 1 



A^ 



(29) 



As one considers higher redshifts z, the cosmic microwave 
background will increase in energy density, since (0(2) = 
(1 + z)*S'o (we note that radiation decouples from matter 
at a redshift - lO'"^ (Peacock, 1998)), where = aT^ de- 
notes the current energy density, and a ~ Sn^kg/lBh^c^ w 
7.6 X 10^^^ J/m'^K'* is the radiation constant. Thus, Eq. ( I29l l 
gives v/c ~ 1 — (2/3)Aa(l + z)*Tq , for the phase velocity v. 
Using To = 2.7 K, a correction of 10% to the phase velocity 
in vacuum is obtained for a redshift Zc ^ 10^°, i.e . rough ly at 
the time for neutrino-matter decoupling ( |Peacocklfl99l . 

Including the dispersive correction, as given by Eq. J25t . 
the dispersion relation for a background of incoherent photons 
takes the form 



k| [1 - fAf? (1 + |(TA<f |k|' 



(30) 



i.e. high-frequency pulses may suffer spectral dilution when 
propagating through a radiation gas. 



5. Other field configurations 

Similar dispersion relations can be found for other back- 
ground field configurations, e.g. plane wave backgrounds 
or partially coherent electromagnetic fields, as is relevant in 
ultra-high intensity laser applications. Of special interest for 
detection purposes is the configuration of photon propagation 
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perpendic ular to a co llection of co nstant elec tric and magnetic 
fields (iBakalov et alj..l99&.Rikken and Rizzo. .2000- 2003 ) 

We should also note that in all the cases above the group 
and phase velocities of the test photons are subluminal, as 
expected, since we have excited the quantum vacuum by 
using electromagnetic fields, analogous to a normal disper- 
sive material medium. This can be contrasted with the 
Scharnhorst effect (Barton, 1990; Barton and ScharnhorsJ 
[1993; Scharnhorst, 1990. ll998h . for which we obtain super- 
luminal phase and group velocities between two conduct- 
ing plates. This can be traced back to the Casimir effect 
JCasimir and Poldeilll948t lc asim hi ll 948ft . where the quanti- 
zation between two conducting plates allows fewer states than 
for field with boundary conditions at infinity. Thus, in this 
sense, the vacuum between the plates is 'emptier' than out- 
side, giving rise to superluminal velocity. 



F. Ultra-intense fields 

The dispersion relations treated so far have used the weak 
field expansion of the general Heisenberg-Euler correction 
il3i . However, both from an application point of view and due 
to theoretical issues, the inclusio n of fully nonlinear vacu um 
effects deserve attention (see also lDittrich and GiesI (l200fll) for 
a thorough discussion of the strong magnetic field case). 

Within astrophysical and cosm ological settings, such as 
neutron stars and magnetars ( Ko uveliotoul the strong 

field conditions can be met. Even in laboratory environ- 
ments, such conditions could be encountered in future high 
energy laser configur ations. While today's lasers can produce 
10^^ - 10^^ W/cm^ dMourouet aL 1998.) . it is expected that 
the next generation laser-plasma systems could reach 10^^ 
W/cm^ (Bingham, 2003; Bingham et al., "2007; ' Cairns et all 
E004) . where field stren gths close to the Schw inger critical 
value could be reached ( iBulanov et all l2003l) . Thus, non- 
linear effects introduced by photon-photon scattering will be 
significant, and the weak field approximation no longer holds. 
In terms of Feynman diagrams the discussion to follow will 
consider the full one-loop correction 



(31) 



The general dispersion relation for a test photon in a vac- 
uum dressed by a strong electromagnetic field is given by 
(t?ittrich andOies, 1998) 



[1 + iAeo(E2 + c^B^ 



2Aeock • (E X B 



-Aeo 



i(E2 + c^B^) - (k • E)2 - c^(k ■ 



V 

c 

- 1, 



where v — u!/\k\ = c/n is the photon phase velocity, and n is 
the refractive index. The effective action ch arge A is no longer 
a constant, but instead defined through (iDittrich and GiesI 
[1998) 



A = 



1 



32 



(33) 



In many cases, the part of the denominator of ( I33> stemming 
from the nonlinear QED correction ( I13> can be neglected, 
since it will be much smaller than the remaining terms, even 
for fields _Eciit- 



1 . Pure magnetic field 

The case of pure magnetic fields enables a simplifica- 
tion of the evaluation of the Lagrangian ( I13l l. The refrac- 
tive index for the strongly magnetized vacuum can be de- 
termined in terms of special functions. For a pure magnetic 
field Bq, we immediately obtain using Eq. WW that a = 
c|B n| and b = 0, respectivel y. Dittrich and Gies (1998) (see 
also lDittrich and GiesI (l2000l) X starting from the Lagrangian 
(I13> . devised the general expression (I33> for the effective 
action charge A (as also can be found in the work of, e.g. 
iBialvnicka-Birula and Bialvnicki-Birulal ( 1197 0)). 

Using a = c|Bo| and 6 = and the approximation A « 
(5^ + d^)^/2eo (see Eq. (|33}), the effective action charge 
takes the form 



27reoc2|Bo|2 



-4a;lnr(x) + 2a;ln27r+ i + 4C'(-l,4a;) + — 
6 Dec 



, (34) 



where x — i?crit/2c|Bo|, ip is the logarithmic derivative of 
the r function, and (' is the derivative of the Hurwitz zeta 
function with respect to the first index. Thus, the refractive 
index, given in Eq. i32\ . becomes 



1 - AeoC^IBopsin^e > 0, 



(35) 



where higher order terms in the fine structure constant a has 
been neglected, and 9 is the angle between the background 
magnetic field Bq and the wavevector k. 

The refractive effects of a super-strong magnetic field is of 
interest in neutron stars and in magnetar environments, since 
they generate extreme conditions in terms of the field strength. 
Comparing with the case presented in Sec. lll.E.TI we see that 
given the magnetic field strength |Bo| ~ 10^^ T, we have 
Aeoc^lBop ~ 15a/n w 0.03. Thus, the effect of the mag- 
netic field on the refractive index is weaker than predicted by 
the lowest order calculation. 



2. Crossed field background 



For a crossed field configuration, i.e. |E| = c|B| and E • 
B = 0, it immediately follows that a ~ b ^ Q. For a test 
photon belonging to an ensemble of random photons (such as 
in a photon gas), we obtain ( Marklun d et al.. 2005b.) 



l + Aeo|E|2 
l-iAeo|Er 



1/2 



(36) 
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from Eq. i32\ . where the effective action charge in i33\ 
is a constant due to the random properties of the test pho- 
tons, = (45/22)(47r/a)eo-Ecrif Note that this is the 
same charge as the geometrical average of the coefficient ob- 
tained from the polarization tensor in the wea k field limit 
(iBialvnicka-Birula and Bialvnicki- Rimla[ll97nl) . I.e. an av- 
erage over polarization states. The refractive index diverges 
as eo|Ep 3A^^. As these field strengths are reached, it 
is not correct that the test radiation gas is in thermodynami- 
cal equilibrium, and the assumptions behind the derivation of 
(I36t are no longer valid. 



3. Incoherent radiation bacl<ground 

We may characterize single photons in terms of plane elec- 
tromagnetic waves. Thus, for an electromagnetic wave mov- 
ing in an isotropic and homogeneous medium with the refrac- 
tive index n we have |B| = n|E|/c, = — cE -6 = 0, and 
^ = — l)|Ep > 0. Here we have assumed that n > 1, 
which implies a = [(n^ — l)J'/eo]^/^ ^ Oand6 = 0, while for 
superluminal velocities, we have a — and b ^ 0, whi ch al- 
lows for spontaneous pair production dSchwingej. 1 1 95 ll) . The 
n > 1 assumption is consistent with elastic photon-photon 
scattering. The vacuum is now treated fully nonlinearly so we 
have to take into account the backreaction of the random pho- 
tons onto themselves, an interaction mediated by the refractive 
index. The effective action charge ( I33> will therefore depend 
on both the field strength and the refractive index. Since the 
refractive index itself depends on the effective action charge 
A, it will be nonlinearly determined via j32> . For incoherent 
photons the Poynting flux in the gas rest frame vanishes, and 
we may uniquely characte rize the gas by it s energy density S". 
From (|33 we then obtain jMarklund et al.U2005bl) 



1 l-|A<f +./l-2A<r + i(A^)2 
— = (37) 

while the effective action charge takes the form 



A 



4^eoa2 2 + {a I S^)[F {a / E,,^) + G{alE,,^)] ' 



(38) 



Here 



F{a/E,n,) ^ lim (9i + d?.)^ 

1 f^°° dz p f 1 — z coth z 1 
= — / _e--f^™^/'^ 1- -zcothz 



(39) 



47r 



sinh z 3 



and 



GialE^iA = — lim [-29^^^ - a^d%^c]- (40) 

Note that the latter function only gives a small correction to 
the effective action charge, and can in most cases safely be ne- 
glected. Moreover, the function F may be expres sed in terms 
of special functions, see Eq. ( I34t (.Dittrich and Gies,.200Q,) . 



e/E , 



FIG. 5 The phase velocity v in units of c, as given by Eqs. <37t and 
Bsl . plotted as a function of the normalised en ergy den sity S /Sa\t 
(reprinted with permission from lMarklund et all 120051J) ). 



The we ak field limit of (l37li and J38> tak e s the 
form M9\ ( Bialvnicka-Birula and Bialvnicki-Bi rulaL Il970t 
.Dittrich and G ies. 1998: Marklund e t al., 2003). On the other 
hand, in the ultra-strong field limit (S'/epE^ ^ oo), we ob- 
tain the asymptotic constant phase velocity i Marklund et all 
l2005b) 



c/V5w 0.45c, 



(41) 



valid in the low frequency approximation. Thus, for very high 
radiation densities, we expect the phase and group velocities 
to be approximately half that of the speed of light in vacuum. 
In this limit, the one-loop radiation gas essentially evolves by 
free streaming, and it is therefore likely that higher-order loop 
corrections are important (Ritus, 197^. The phase velocity as 
a function of intensity is depicted in Fig. |5] 

Under most circumstances, the contribution proportional to 
a in the denominator of Eq. ( I38t is small and may be ne- 
glected. When (o > epE^^,, we obtain F » o/BE'crit from 
Eq. (|39|l, and we have ( IMarklund et"aill2005bh 



;F{a/E,nt) 



1 



1 



(42) 

This expression for the effective action charge can be used 
in conjunction with (I37> to analyse the propagation of single 
photons in a radiation gas where the energy density is in the 
intermediate range. 



III. NONLINEAR COLLECTIVE PHOTON INTERACTIONS 

In the preceding section, we presented the dispersion rela- 
tions for single photons interacting nonlinearly with the vac- 
uum, for a number of special cases. The single-photon picture 
contains many interesting physical phenomena, such as pho- 
ton splitting and the birefringence of the quantum vacuum. 
However, in many appl ications, collectiv e effects among pho- 
tons may be dominant dMendoncalEool . 
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A. Coherent field interactions 

The formulation of the interaction between coherent elec- 
tromagnetic waves and possible background field configura- 
tions is a starting point for the discussion concerning pos- 
sible detection techniques of elastic photon-photon colli- 
sions. Experiments for detecting elastic photon-photon scat- 
tering are important tests of QED, and constitute a new 
type of tests of, e.g. Lorentz invariance in extensions 
of the Standard Model such as supersymmetric field theo- 
ries (Colloday and Kosteleckv, 1998; Jackiw and Kosteleckil 
[1999; Lipa et al., 2003; Nibb elink and Pospelovt 120051) . 

The Maxwell equations that results from the weak field 
Heisenberg-Euler corrected electromagnetic Lagrangian (|5j 



V-E 
V B 

ot 



(p-V.P)/eo, 
0, 



0, 



Mo J 



- + VXM 



(43a) 
(43b) 

(43c) 
,(43d) 



where the vacuum polarization and magnetization are given 
by the expressions il 1> and il2i . respectively, and p and j are 
the charge and current densities, respectively. From these, it 
is straightforward to derive the wave equations (|9} and dlOt . 
These may in turn be used to derive the dispersion function 
for the appropriate wave field on a given background. 

As noted above, the interaction between waves in parallel 
propagation does not yield any interaction due to the vacuum 
dispersion function I?vac = o;^ — |kpc^ appearing as an 
overall factor in the wave equations (|9} and ( I10> . Thus, dis- 
persive effects need to be introduced in the wave propagation. 
This can be done in a multitude of ways, such as crossing 
light beams (Soliacic and Segev, 2000b), photon propagation 
on a constant coherent field bac kground ( Bakalov et aU 
|l998^■'Ding and Kaplanl ll98lll992tlKaplan and Ding..200a 
Hikken and Rizzo, 2000 . l2003t iRozanovl 11 993L 1 1 998h ^ cav- 
ity fields (Brodinet al., 2001, 2002), wave guide prop a gation 
fBrodin et al.. 2003: .Shen and Yu. .2003; .Shen et all 12003 
Shukla et all l2004bl). plasma in teracti o ns (iMarklund et al 



2004b[ l2 005d' S hen and"^ l2003t IShen et all 12003 
Stenflo et al. . 20_05), and interaction between coherent and 



incoherent photons (Mai'klu nd et all l2003l l2004ald l2005d: 
IShukla and EhassonLi2004ai;lShukla et alll2004cldl) 



1. Nonlinear vacuum magneto-optics 

We have seen that the propagation of photons on a magne- 
tized background can be expressed according to the dispersion 



^ Here it should be noted that the paper by 'Bakalov et alj Il998h is a 
progress report to one of the few actual experimental setups within photon- 
photon scattering, and their d etec tion techniques are ba sed on the work by 
llacopini and Zavattinil 1 19791) and lBakalov et alj <1994) . 



relation ( I26> . We may also start from the constituent relations 
(|8} together with the expressions il l\ and (I12t for the polar- 
ization and magnetization, respectively. This was first done by 
Klein and Nigam ( 1964a b), and later for arbitrary intensities 
by I Icvl and Hernguist ( 1997a) . Denote the slowly varying 
background magnetic field by Bq, and the perturbation fields 
by E and B. Then, we find that Di — CijEj and Hi — HijBj, 
where the quantum vacuum electric permittivity and magnetic 
permeability are given by 



IJ-ij = Mo [Sij 



(44) 
(45) 



respectively. Here, we have introduced the background mag- 
netic field energy density = I Bo P//io and the background 
magnetic field direction b — Bo/|Bo|. Thus, the permittivity 
and permeability are diagonal when using the magneto-optical 
axis as eigen-direction. Denoting this direction by z, we have 
<^xx = <^yy = eo(l - 4k^o) and e^^ = eo(l + IOk^o), while 
fJ'xx = Myy = Mo(l + 4k^o) and fi^^ = fio{l + 12k3§o). 

For an electromagnetic wave propagating perpendicular to 
Bq, there are essentially two different polarization states, and 
we may write eij — eSij, fiij — ^Sij . When E ± Bq, we 
have e = e^, /i = /i^ according to 

e± = eo(l - 4k^o), and = /io(l + 12k^o), (46) 

while if E||Bo, we find e = e||, /i = where 



eo(l + IOk^o), 



and 



=/io(l + 4K/^o)- (47) 



Thus, we see that for strong magnetic fields, there is a 
significant difference in the behavior of the two polariza- 
tion r nodes. This has been expl o ited in various sc e narios 
(e.g. iBakalov et all (Il994i Il998t) : iDing and Kaplar^ ("19891 
1992h:lHevl and Hernauistl(ll997ah : Kaplan and Ding (2000); 
Rikken and Rizzo (2000, 2003)). The procedure is straight- 
forward to perform for other uni-directional background field 
configuration. 



2. Nonlinear self-interactions 

Based on these results, the self-action of an electromag- 
netic pulse on a given background can also be consid- 
ered. Taking into account the lowest order cubic nonlin- 
ear terms of the (complex) pulse amplitude E, and employ- 
ing the slowly varying envelope a pproximation ( Hasegawa 
1 19751 iKivshar and Agrawall l2003l) . one can derive a non- 
linear Schrodinger equation (NLSE) for E. Letting E = 
J Ek exp[z(k — ko) ■ r — i{uj — ujQ)t] dk, and expanding the 
frequency around the background values (denoted by 0), we 
obtain 



a; « cjQ + 



duj 



doj 
dh 



[h - fcoi) + 



1 d^uj 



2 dk.jdkj 



5|Q| 



(IQP-IQoP), 



(48) 
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where |Qp is given by Eq. 
satisfy the NLSE 



d 



EOb . Thus, the envelope will 



1 



I OKOj Z 

(49) 

where Vg is the group velocity and /(|Qp) = /i?fe(|Qp 



be polarized in the x-directio n. The lowest order soli- 
tary wave solution is given by ( iKivshar and Agrawall l2003^ 
A{x) — Aq sech (\/^2r z) where we have denoted the eigen- 
value r — — fc^KEo^o/^- This one-dimensional soliton so- 



lution is stable, as opt 


)osed to higher dimensional solitons 


dBerse and RasmusseJ 


[19961 


Kiyshar and Aarawal, 


2003). 


Furthermore, Soliacic and Segev] 


(l2000b) also suesested the 



QoH exp[i(k - ko) • r - i{u; - LOo)t]dk is the nonhnear re- possibility of hig her dimension al soht on formation, e.g 



neckla ce sohtons llSoljacic et al 1 119981: ISoliacic and Segevl 
I2000al). 



sponse. We note that the term containing |Qo|^ represents a 
phase shift, and can be removed by a transformation. When 
high-frequency corrections are added, the NLSE will attain a 
second order derivative along the propagation direction. The 
group velocity Vg — doj/d'k on an arbitrary background 
can be written as (Bialvnicka-Birula and Bialvnicki-Birula, 
= ck- (cAeo/2)[|E|2 + c^\B\^ + (k • E)^ + c^k ■ 



Il970h ■ 



B)2]k - cAeo[(k • E) E + c2(k • B) B + c E X B] in the weak 
field case. We note that \vg\ ~ uj/k < c. We will now discuss 
some special cases. 



a. Constant background Gelds For a constant background con- 
figuration, the back reaction of the photon propagation may, in 
some cases, be neglected. Then, the nonlinear contribution to 
the self-interaction of the pulse will come through a coupling 
o f higher or d er in t he parameter A. 

iRozanovl (|ip98') considered the perpendicular propaga- 
tion of high intensity laser pulses on a background Eg — 
Eqx, Bo = Boy. By choosing the polarization directions of 
the laser pulse in the direction of the background fields, one 
obtains the NLSE^^ 



d 



d 



dt 



dz 



1 



(50) 



where v'g — c{\ — A|QoP/2)/fco is the group velocity disper- 
sion, e = fcoceoA3|Qo|V8, and IQop = ^^{Eo - cBof. The 
NLSE (I50> could be of interest for laboratory applications, 
when studying high intensity laser pulse propagation in given 
background electromagnetic fields. 



b. Crossing beams When the background is given by the 
source itself, the nonlinear self-interaction term will be of first 
or der in A, thus requ iring weaker background conditions. 

ISoliacic and Segev. (,2000b) derived a NLSE for the dynam- 
ics of the envelope A{x) of the interaction region due to cross- 
ing laser beams. By symmetry arguments concerning the 
QED corrected Maxwell's equations, they reduce the problem 
to a ID stationary NLSE 



fA 



(51) 



where k is the wavenumber of the laser beams, F is the 
eigenvalue of the equation, and the beams are assumed to 



^ We note that ^Rozano^l . ll998^ obtained a the NLSE 1501 with the dispersive 
correction. 



3. Propagation between conducting planes 

Similar to the case of a Casimir vacuum, one of the sim- 
plest geometries where dispersive effects makes the presence 
of QED vacuum nonlinearities apparent, is given by two par- 
allel conducting planes. They are the first example where 
multi-dimensional photon configurations can self-compress 
to reach intensities above the lase r limit ( Br odin et al., 20031 
iMourou et all (19981 IShnkIa et allE004hl). 



a. Variational formulation Consider the propagation between 
two parallel conducting planes with spacing xq of one TE^q- 
mode (J, = 1, 2, ...) given by 



. flTlx\ 
sm 



A sin ( \ exp[i(fcz — wt)]y + c.c. 



(52) 



in the radiation gauge (0 — 0). The linear dispersion relation 
is w^/c^ — A;^ — ^^tt^/xq — 0. From Maxwell's equations ( l43l 
a nonlinear dispersion relation can be derived by inserting the 
linear expression for the fields and separating into orthogonal 
trigonometric functions. The coefficients in the NLSE can be 
found from the resulting equation. One may also start from the 
Heisenberg-Euler Lagrangian (|5Jl, and minimize the resulting 
expression for the act ion. This may appe ar more elegant and 
gives the sa me result (lBrodi gjL^lJ, E003l) . 

We follow iBrodin et al.l(l2003l) . Let A A{t, y, z) and as- 
sume A to be weakly modulated so that |9A/9t| <C \ojA\, 
\dA/dz\ <^ \kA\. To lowest order, the nonlinear terms and 
the slow derivatives in ^ are omitted. Averaging over the 
plate spacing, xq, shows that this lowest order Lagrangian is 
identically zero due to the dispersion relation. To the next 
order of approximation in the Lagrangian, first order slow 
derivatives are included. After variation of the corresponding 
action, this leads to an equation where the envelope moves 
with the group velocity. The next order and final approxima- 
tion includes second order slow derivatives. After performing 
the averaging between the plate inter-spacing, the final expres- 
sion for the Lagrangian is 



'dA 
:.So I -^A 



dA* 
~dt 

dA 



A 



., 2 i'dA dA* 
oz oz 



A 



dz 



\A\' 



(53) 
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Variation of the action due to the Lagrangian i53\ with respect 
to A* leads to the NLSE 



d_ 
dt 



d 



S'A 



2lj dy^ 



Vg'd^A 

2 dz^ 



L^\A\'^A = 0, 
(54) 

where Vg and v'g ~ dvg/dk follow from the linear disper- 
sion relation, and = (3£'^c'^tt'^ neo) / (ujxq). The nonlin- 
ear correction in Eq. i54\ is due to the self-interaction of the 
TE^o-mode i52\ via the quantum vacuum. In one space di- 
mension, i.e. d^A/dy^ = 0, Eq. (I54> reduces to the cubic 
Schrodinger equation which admits an envelope soliton solu- 
tion (^vshar and Agrawal, 2003). 

Changing to a system moving with the group velocity while 
rescaling the coordinates and the amplitude according to t = 
u;t/2, V = ujy/c, C = {uj /v'gY/'^{z - Vgt), and a = V2LA, 
and assuming cylindrical symmetry, Eq. takes the form 



.da Id 



^ dr pdp \ dp 



da 



+ \a\^a = 0, 



(55) 



where a = a{t, p), and = + 

Equation ( I55> is a 2-dimensional radially symmetric NLSE, 
to which exact solutions are not available. However, an 
accurate analytical approximation of the dynamics of the 
pulse-like solutions of Eq. i55\ can be obtained by means 
of Rayleigh-Ritz opti mization based on su it ably chosen 
trial functions (see e.g. Ande rson et al.l ('1999a); Desaix et al. 
(^^91j and references therein). An accurate approximate so- 
lution, mimicking the solitary behavior as wel l as capturing 
the collapse properties of Eq. ( I55> is given by ( iDesaix et all 
1199 1.) 



arir, p) = F{t) sech 



fir) 



exp [i&(r)p^ 



(56) 



where^ /(r) = [/^(O) + 7 (1 - h/Ic) t^Y'^, showing the 
instability of the stationary solution /o = /c, either collapsing 
to zero width in a finite time when /q > /c, or diffracting 
monotonously towards infinite width when /o < la- 
in the next section, a perturbation analysis will show that 
the exact equations produced unstable solutions. 



b. Instability analysis Following IShukla et all ( l2004bh . a 
rescaUng of Eq. i54\ gives the dimensionless equation 



dt 



n ^ 9 \ , l&'A 13^ d^ A , ^,2 , „ 



di 



2 9y2 



2 dz 



(57) 

where (3 = Ittc/xquj, and the time is scaled by oj^^, 
the spatial variables by c/oj, and the vector potential by 
(«eo/2c.2)V2. 



* The complex amplitude F{t) an d the phase function h jr) can be expressed 
in terms of the pulse width /(r) IDesaix et alill99ll) . Here 7 = 4(ln 4 + 
1)/(27C(3)) « 0.29, /(r) = /2(r)|F(r)|2 = /2(0)|F(0)P = 7o, 
and /c = (21n2 + l)/(41n2 - 1) 1.35. 



Conditions for the modulational and filamentation instabil- 
ities can be obtained as follows. Let A — [Aq + Ai exp(«0) + 
A2 exp(— 10)] exp(— iwot), where = K • r — i7t is a phase, 
bjQ is a constant frequency and the constants Aq, Ai and 
A2 are the complex amplitudes of the pump wave and the 
two electromagnetic sidebands, respectively. The wavevector 
and the frequency of modulating perturbations are denoted by 
K ~ yKy + and SI, respectively, where y and z are the 
unit vectors along x and y axes, respectively. Following the 
standard proced ure of the modulational/filamentational insta- 
bility ( Anderson et all fT999b: Shu klaet al.. 1986) . the non- 
linear dispersion relation becomes 

+ [\Ao\' - i {K^g + P^Kl)] {Kl + P^Kl) - 0. (58) 



Letting fi = + 17 in ( I58> . one obtains the mod- 

ulational instability growth rate 



7 



[\A,\^~\{Kl+p^Kl)\ 



,1/2 



P^Klf^. (59) 



The instability grows quadratically with the amplitude Aq, 
and attains a maximum value at a critical wavenumber. Values 
of j3 different from unity make the instability region asymmet- 
ric with respect to Ky and K^- On the other hand, the spatial 
amplification rate F = iKz of the filamentation instability in 
the quasi-stationary limit (viz. VL = {)) and for 0^r'^ <C Kg is 



r = {\A, 



lr^2Nl/2 Kg 



(60) 



'Shu kla et alJ ( l2004b l) performed a numerical study of Eq. 
( I58> showing the instabilities indicated by the approximate 
solution ( I56> . see Fig. |6l Using the normalized Eq. ( 15 7> 
an initially Gaussian pulse was shown to col lapse or dis- 
perse in accordance to the collapse criterion in iBrodin et alJ 
(2Q03\ see Fig. Q The collapse is unbounded in the 
weakly nonlinear model given by ( 15 7> . As the collapse pur- 
sues, the intensity of the pulse will reach values at which 
the weakly nonlinear theory breaks down and higher or- 
der effects (Bialvnicka-Birula and Bialvnicki-Birula, 197(J 

trklund et al., 2005b), and possible pair creation processes 
hwingert il951 ), has to be taken into account. For the lat- 
ter, a significant energy dissipation into the electron-positron 
plasma will take place. 

In practice, the trapping of an electromagnetic pulse could 
be achieved by two highly conduction layers. As an exam- 
ple, we consider the generation of two-dimensional plasma 
chann els due to the intera ction of a plasma with high-intensity 
lasers dShen et ail l2003h . In this case, a vacuum will be cre- 
ated within the plasma due to the complete evacuation of the 
electrons by the ponderomotive force of intense laser beams. 
The resulti ng plasma waveg uides can sustain very high field 
intensities (Bingham', '2003'), and with future laser systems 
(Mourou et al., 1998, 2005) the intensities could surpass even 
the th eoretical laser hm it (Mourou et allll998HShen and Yul 
I2OO2I: Ishorokhov et all l2b03i,) . The trapping of intense laser 
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FIG. 6 The growth rate 7 as given by <59> of the modulational insta- 
bihty for an initially homogeneous radiation field as a function of the 
wavenumber {Ky, Kz), fo r different value s of /3 and pump strength 
^0- (Reprinted from Shu kla et alj j2004d) . Copyright (2004), with 
permission from Elsevier.) 
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and'Br odinetalJ (l2002l) were able to derive a NLSE for the 
self-interaction of a single mode in a rectangular cavity. 

Calculations of the three-wave coupling strength between 
various eigenmodes can be made including the nonlin ear po- 
larizat ion (II 1> and magnetization (I12> . see e.g. Brodin et alJ 
(l200lh . However, a m ore convenien t and e legant approach, 
which was pioneered bv' Brodin et alJ ( l2003h and which gives 
the same result, starts directly with the Lagrangian density (|5}- 

A general procedure for finding the cavity eigenmode cou- 
pling and the saturated amplitudes of the excited mode can be 
formulated (B rodin et al., 2001, 2003; Eriksson et al., 200l 
IWeiland and WilhelmssonLll977^ : 

1. Determine the linear eigenmodes of the cavity in terms 
of the vector potential . 

2. Choose resonant eigenmodes fulfilling frequency match- 
ing conditions for modes 1, 2, and 3, such as 



W3 = 2lJi — UJ2- 



(61) 



3. Assume a slowly varying amplitude of the vector poten- 
tial eigenmode amplitudes and minimize the effective action 
obtained from the Lagrangian (|5} and follow steps 1 and 2. 
The lowest order linear terms vanish, since the dispersion re- 
lation of each mode is fulfilled. 

4. In the absence of dissipation, the mode coupling equa- 
tions imply steady growth of mode 3, until the energy of that 
mode is comparable to that of the pump modes. A damp- 
ing mechanism, such as finite conductivity of the cavity walls, 
may be inserted on phenomenological grounds. This saturates 
the amplitude at a level depending on the the mode-coupling 
growth versus losses. 



FIG. 7 The time-evolution of jj4| as given by Eq. ^ for an ini- 
tially Gaussian shaped electromagnetic pulse. The pulse slowly self- 
focuses and finally collapses. Here /3 = 0.5 and a Gaussian pulse 
A = 10~^ exp[-(y^ + z^)/{2 x 10^)^] is used. In the upper right 
and lower left panels, the pulse self-compresses and in the lower right 
panel the field strength of the pulse has reached a critical limit where 
Eq. O is no longer valid. (Reprinted from Shukla et aL (■2Q04a) . 
Copyright (2004), with permission from Elsevier.) 



fields could yield the right conditions for the electromagnetic 
modes to self-interact via the nonlinear quantum vacuum, giv- 
ing rise to pulse evolution according to Eq. ( I57> . 



4. Cavity mode interactions 

Wave-wave interactions can give rise to a host of in- 
teresting phenomen a, well-k nown in optics and plasma 



physics fAgrawal', '2001"; 'Kivshar and Agrawaf, l2003t 
[Weiland and Wilhelmsson, 1977). As cubic nonlinearities act 
within a cavity environment, they will produce wave-wave 
couplings, and given certain resonance conditions, a new 
mode wil l be generated, that will satisfy the cavity dispersion 
relation. iBrodin et alJ (|2001) showed that these new modes 
could reach detectable levels within state of the art cavities. 



a. Rectangular cavities For a rectangular prism cavity with di- 
mensions {xq, uo, zq), choosing the radiation gauge, the pump 
modes have vector potentials of the form 



Aj — Aj sin 



exp{—iujjt)y + c.c, (62) 



where j — 1, 2, £j = 1, 2, 3... are the mode numbers the pump 
waves, and xi ^ x, X2 ~ y, xio — xo, and 2:20 = Uo- The 
dispersion relations are 



The mode excited due to the QED nonlinearities is given by 



(63) 



k-3 



where 



A^ sm I — I sm 

yo 



Zo 



exp(— ict;3t)x + C.C., (64) 



(65) 



Followit ig the scheme as gi ven by iBrodin et alJ ( 1200 ll 
l2003h and lEriksson et alJ (2004), with the resonance condi- 
tion = 2uji — UJ2, one obtains the evolution equation for 
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mode 3 



dt 



' K A"^ A* 



2 2 



111 



0^2^3(- 



-W2) 



(66) 



where the dimensionless coupling coefficient K^^c is 



(67) 



-0.4 y 0.4 0.8 1.2 



(a) 



The different signs in the expression (|^} for the coupling 
strength correspond to the mode number matchings 2£i — ^2 + 
£3 = 0, 2£i + £2 - £3 = 0, and 2li - £2 - £3 = 0, respectively, 
that must be fulfilled for nonzero coupling. The coupling co- 
efficient for specific mode numbers and geometries can thus 
be evaluated. If a saturation mechanism is included in ( I66I 1. 
one may solve for the steady-state value of ^3. 



b. Cyli ndrical cavities As was shown by lEriksson et al.l 
(l2004ft . the efficiency of the mode conversion can be slightly 
improved by the choice of a cylindrical cavity. The results 
can be obtained along the lines of the previous example. For 
TE-modes with no angular dependence, the vector potential 



£ttz 

A = AJi{pj3/a) sin ( ) exp(— ici;t)(^ + c.c. 



(68) 



gives a complete description of the fields. Here a is the cylin- 
der radius, zq the length of the cavity, Ji the first order Bessel 
function, t is the mode number, and /3 one of its zeros. The 
cylinder occupies the region < z < centered around 
the 2-axis. We have here introduced cylindrical coordinates p 
and z as well as the unit vector ip in the azimuthal direction. 



The eigenfrequency is given by w 



^[(/?/a)2 + (£7r/zo) 



for all modes uj — ^1^2, 3- From the matching condition 
ijj'i — 2uji — UJ2 it follows that all the eigenmodes cannot have 
the same order of thek respective /3, and one thus introduces 
(3 — /3i,2,3- Proceeding along the lines of the previous section, 
one obtains 



dt 



" J ic , A-^ A* 



(69) 



for the mode number matching £3 = 2ii + £2- Here the cylin- 
drical coupling coefficient A'cyi is defined in terms of integrals 
of Bessel functions and can be found in Eriksson et al. (2004). 
As in the case of a rectangular geometry, the linear growth of 
A3 as dictated by Eq. ( I69l l will be saturated by dissipative 
mechanisms. The intensity of the generated field amplitudes 
and the pump field amplitudes is shown in Fig. |8] 



B. Incoherent field interactions 

Above we have considered the interaction via the nonlinear 
quantum vacuum between coherent electromagnetic waves. 



-0.4 , 0.4 o.a 



(b) 
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FIG. 8 The m ode structure for the filtering geometry (figure adopted 
from lEriksson et al.. l2Q04)). The small region to the right is the en- 
trance region, the large middle region is where the interaction takes 
place, and the region to the left is the filtering region. Here In li?! is 
plotted using an arbitrary normalization. We have the cavity distance 
z and the cylindrical radius p in units of meter, (a) Pump mode 1 : 
We see that spatial exponential decay has diminished the amplitude 
by a factor 10~® in the filtered region to the left, (b) Pump mode 
2: As for mode 1, exponential decay diminishes the amplitude in the 
filtering region, here by a factor 10^*. (c) Excited mode: The ampli- 
tude of the mode excited by the QED vacuum nonlinearities is almost 
unaffected when passing into the filtering region. 



However, in many situations where the Heisenberg-Euler La- 
grangian is important, such as in astrophysical and laser- 
plasma applications, there can be intense incoherent electro- 
magnetic fields present. We will below study two scenarios. 
First, a plane wave pulse propagating on a vacuum dressed 
by an intense gas of incoherent photons is analysed, and, sec- 
ondly, the effects on a radiation gas by the quantum vacuum 
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excited by an intense electromagnetic (EM) pulse is investi- 
gated. These two results are then used in conjunction to ob- 
tain the relevant equations governing the nonlinear interaction 
between the pulse and the radiation gas. 



1 . Coherent pulse interaction with incoherent photons 

The dispersion relation i30\ represents the propagation 
of test photons in an intense incoherent radiation back- 
ground. Using standard methods for slowly varying envelopes 
iHasegawa, 1975) in conjunction with (13 0> . we obtain the spe- 
cial case of (I49> . f or a p ulse in an inte nse photon g as back- 
ground (Marklund etaill2003L l2004at IShukla et all l2004d) 
(see also iRozanovl d 19981) for a similar result using a differ- 
ent strong background field) 



+ v 



d_ 



En 



_9_ 

2fco 



vi 



a2 



(70) 

where 5S' — — (fg is a radiation gas perturbation due to 
the pulse propagation and S'q is the unperturbed background 
radiation energy density. Here we have adapted the coordi- 
nates such that ko = fcoz, ko being the background wavevec- 
tor of the pulse, included the high-frequency correction repre- 
sented by Pz, and denoted the (complex) pulse amplitude by 
Ep. Moreover, Vg = (9a-'/9k)o is the group velocity on the 
background, Pz — (32/3)(7A^(fQ^fcQ is the vacuum dispersion 
parameter, /i = {2/'i)ckQ\[l + (16/3)(TA(oofco] nonlin- 
ear refraction parameter, = — (ko • V)^, and A = 8k 
of 14k depending on the polarization of the pulse (see Sec. 

pro. 



evolution of the radiation gas ^. Taking the moments of the 
kinetic equation M2\ (Marklund et al.^ ■2004a) . we obtain the 
energy conservation equation 



'dt 



V • (^fu + q) 



n dt ' 



(73a) 



where (S'{r,t) — Jhutfdk is the energy density, and 
q(r, t) = J hjj-wf dk is the energy (or Poynting) flux. Here 
we have introduced f = u + w, where / w/ dk — 0. Thus w 
represents the random velocity of the photons. Equation ( 17 3 al l 
is coupled to the momentum conservation equation 



'dt 



u«) n 



-Vn, 



(73b) 



where U = J Kkf dk is the momentum density and P = 
Jw (g) (ftk)/(ik is the pressure tensor. It follows immedi- 
ately from the definition of the pressure tensor that the trace 
satisfies Tr P = J hkw • k/ dk = ("^/c) / hujw ■ k/ dk. For 
an observer comoving with the fluid, i.e. a system in which 
(u)o = (the denoting the comoving system), Eq. ( 17 1> 
shows that (w • k)o = (n)o/c, so that the trace of the pres- 
sure tensor in the comoving system becomes (Tr P)o = (<?')o- 
For an isotropic distribution function, the pressure can be ex- 
pressed in terms of the scalar function P = Tr P/3, satisfying 
the equation of state P = S'/3. We will henceforth adopt the 
comoving frame, in which u = 0, and the equation of state 
Pij — PSijSijS' /3, in order to achieve closure of the fluid 
equations. 



3. Quasi-linear theory 



2. Radiation gas response of pulse propagation 

We saw that the effects of a plane wave on incoherent pho- 
tons could be expressed via the dispersion relation ( I28> . Fol- 
lowing iMarklund et al.l (l2003i) . the response of the radiation 
gas can be determined using a kinetic theory. For a dispersion 
relation ui = cfc/n(r, t), where n is the spacetime dependent 
refractive index, we have the Hamiltonian ray equations 



duj c - • 
r = — = - k, and k 
ok n 



-\/uj = -Vn, 
n 



(71) 



where r denotes the group velocity of the photon, k the force 
on a photon, and the dot denotes a time derivative. 

The equation for th e collective interaction of photons can 
then be formulated as dMendoncalEool 



dt 



-V-(r/(k,r,i))- 



|--(k/(k,r,t))-0, (72) 



where the distribution function /(k, r, t) has been normalized 
such that J /(k, r, t) dk gives the photon gas number density. 
In what follows, we will neglect the dispersive effects on the 



We now assume that the radiation gas is perturbed around 
the equilibrium state (oq = constant and Ho = 0, letting S = 
<oo + where |(5(o'| <C <§q- Then, using Eqs. ( I73t . we obtain 
an acoustic equation 



dt^ 



5S 



2Aeo(?o 



a<2 



\En 



(74) 

to lowest order in 5S'. This gives the dynamics of a ra- 
diation gas due to the pulse propagation ( Marklund et aQ 
2004a ). Equatio n s M0\ and MAI were first presented 
by iMarklund et alJ (l2004ah . and generalize the Marklund- 
Brodin-Stenflo (MBS) equations IMarklund et al., 2003), to 
the case of a dispersive vacuum. The MBS equations are dif- 
ferent from the Karpman equations cKarpmark .1971. .199^ . 
that govern the dynamics of small amplitude nonlinearly in- 
teracting electromagnetic waves and ion-sound waves driven 
by the radiation pressure in an electron-ion plasma, due to the 



' The dispersive con'ection, due to the variations in the photon-field, will give 
rise to higher order effects in the final fluid equation for the energy density 
of the radiation gas. Thus, we may at this stage neglect the dispersive term 
from the fluid equations. 
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difference in the driving term on the right hand side of Eq. 

(|71- 

The dispers i on-free case admits pulse collapse 
llMarklund et all l2003t LShukla and Ehassoni |2004a), and 
similar features appear within the dispersive case, with the 
difference that pulse splitting may occur, resulting in a train 
of ultra-short pulses. If the time response of the radiation 
background is slow, Eq. ilAi may be integrated to yield 
SS" w 2A(?oeo|-E'pp, and from Eq. (^Oj we obtain the standard 
equation for analyzing ultra-short intense pulses in normal 
dispersive focusing media, see Chernev and Petrov ( 19 93); 
lOaeta (2003); Kivshar and Agrawall ( 120031) ; lRothenber3 
(11992) : Zharova et al.. 120031) and references therein. It 
is well known that the evolution of a pulse within this 
equation displays fi rst self-focusing, then pulse splitting 
dChernev and Petrovl fl992: JLothenberg. 1992) . and the ap- 
proximate description of the solutions c an be given as a prod- 
uct of bright and dark soliton solutions dHavata and KoshibaL 
[1993). A modulational in stability can be found as well 
jKivshar and Agrawalll2003l) . 



4. Instability analysis 

a. The two-dimensional case As 0z goes to ze ro in Eq. ilOi . 
we regain the MBS equations dMarklund et a l.. 2003). In this 
case, the dispersion relation for the modulational and fila- 
mentational instabilities of a constant amplitude photon pump 
(ujq, kp) can be found by linearizing the simplified set of Eqs. 
( I70> and ( I74> around the unperturbed state Ep = Eq — 
real constant and SS" = 0. Followin g the standard proce- 
dure of parametric instability analysis ( iKivshar a nd Agrawall 
Eo03: Shuklai, 1992 ) (see also next section), we consider per- 
turbations varying according to exp[i(K • r — f]t)]. Here i? 
and K are the frequency and wavevector of the acoustic-like 
disturbances. The n, Eqs. i7Q\ and (I74|, yield the nonlinear 
dispers ion relation ( fshukla and Eliassonl l2004aHShuklaet all 
I2004d) 



4fc2 



(3i72 - K^c^) 



^^{n'+K^c')X'<^,eoEl 



(75) 



where ^ + + = R'l + K^. Defining 
Q = cKz + i^2D, \12d\ ^ cKz, we obtain the approximate 
modulational instability growth rate 
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(76) 



where X2D - {8/3)cko\''^QeoEl 

Similarly, in the quasi-stationary limit, SI — 0, a filamenta- 
tional instability may occur. For <C Kj_, we obtain 



Thus, filamentation of an intense photon beam on a radiation 

background takes place when eo I i?poP > 3if^/(16fcQA^(fo)^ 
due to elastic photon-photon scattering. 

In the case of a slow acoustic response, we may integrate 
Eq. (I74> in comoving coordinates, to obtain a relation be- 
tween the radiation gas perturbation and the pulse intensity. 
By inserting the relation into Eq. i70\ (with — 0), the 
equation 



.dEp 



\''ckoSo£'o\EXEp 



0, (78) 



is obtained. The collapse properties of Eq. ( I78> can 
be obtained by approximate analytical means. Start- 
ing from an two-dimensional approximately Gaussian 
pulse Ep — A{t) sech{r±/a{T)) exp[ib{T)r'^], where 
+ j7^, a n approximate solution can be found 



dPesaix et all 119911) . The relation |A|/|Ao| = ao/a, 
where the demotes the initial value. Moreover, the 
parameter 7 = (4/3)A^eofSbfcol^opao(^i/-^2), where 
= xsech'^{x)dx = (4 In 2 - l)/6 and h = 
sech^ (x) dx = 9C(3)/8, and C is the Riemann zeta 
function, characterizes the critical behaviour of the solution 
in terms of the initial data. The collapse criteria can be seen 
inFig.|9l 




FIG. 9 The pulse width as a function of normalized time (r ~+ cfcor) 
in the two-dimensional case. Note that the solitary solution (7 = 1) 
is unstable. 

Exact results regarding the two-dimensional modu- 
lational and filamentational instabilities was found by 
Ishukla and Eliasson (2004a), where numerical solutions of 
Eq. ( I75> were presented, see Figs. 1101 andlTTI The growth 
of random seeds on a radiation gas background was also in- 
vestigated and is shown in Fig. [21 



b. The three-dimensional case We will now show the pres- 
ence of modulational and filamentational unstable modes for 
the three-dimensional case, as given by Eqs. ilOi and ( I74l i. 
Fourier analysing Eq. ( I74t according to SS" and \Ep\'^ cx 
exp[i(K • r — nt)], we obtain 



^1 4 2 2 2 



(77) 



6^ = -X^oWeo\Ep\\ 



(79) 
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FIG. 10 The growth rate and real part of the frequency, respectively, 
as a function the orthogonal wave number for different values of 
the parallel wave number (see Eq. <76» . The dimensionless pump 
strength is = A^rgpeoi^n ~ 0.02 (reprinted with permission from 
IShukla and EliassorJ i20043) '). 
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FIG. 12 The two-dimensional evolution, as given by Eqs. |70l 
and <74t with /3z = 0, of the dimensionless pulse energy density 
A^jfoeolSpp for initially random perturbations on a radiation gas 
background. The color bars give the dimensionless pulse energy den- 
sity (reprinted with permission from .Shukla and Eliassoa 12004a) '). 
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FIG. 1 1 The growth rate and the real part of the frequency, respec- 
tively, as a function the orthogonal wave number for different dimen- 
sionless pump streng ths A = X^SoeoEn (see Eg . <76t ) (reprinted 
with permission from lShukla and EUassoj <2004j) '). 



where W = {Q^ + c^K^)/{~Q^ + c^K^/i), with Q 
and K the frequency and wavevector, respect i vely, o f the 
Fou rier component. Nex t , follo wing IShuklal (1 19921) (see 
also iKivshar and Agrawall (I2003I) '). we let Ep = {Eq + 
El) ex.p{i5t), where i5 is the nonlinear phase shift and Eq 
i':^ \Ei\) is a real constant. To zeroth order in Ei we 
have the nonlinear phase shift 6 = —kEq. We let Ei = 
di exp[i(K • r — f2t)] + d2 exp[— «(K • r — Qt)], with di and 
d2 real constants. Linearizing Eq. 1^0} with respect to Ei, we 
obtain a coupled system of equations for di and ^2- EHminat- 
ing_di and ^2, we obtain the dispersion relation ( Shukla et all 
I2004d) 



[Q - K^Vg 



2ko 



2fco 



(80) 



where x = 4/iA(loeo-E'o- Remembering that W depends on 
the perturbation frequency and wavevector, we see that the 



solution to Eq. (jSUJl in terms of Q is nontrivial. 

Letting Q = K^Vg + 17^ in ( I8OI 1. \^m\ K^Vg, we obtain 
the approximate modulational instability growth rate 




Kl + {l + vl/c^)Kl 
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2ko 
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(81) 



Thus, when Vg sa c, we see that, unlike the standard modula- 
tional instability, we have larger growth rate for smaller length 
scale, with the occurring asymptotically for K±^/2Kz, where 
the approximate expression ( 18 U diverges. 

If the perturbations are stationary, Eq. dSOt yields 



±- 



2fco^ 
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1/2 



(82) 



when PzKl ^ K\, and we see that a filamentation instability 
will occur for x > VgK'j_/2kQ. 

Solving ( I80> for the growth rate, one obtains the insta- 
bility regions show n in Fig . 1131 We note that the re- 
sults found by .Shukla et al.l ll2004d) concerning the mod- 
ulational instability are simi lar to the ones obtained by 
'Karpm an and Washimil (11977'), where it was found that the 
largest growth rates are due to parametric instabilities for 
wave vector s obliqu e to the pulse propagation direction. 
IShukla et al.l (l200 4c*) performed a numerical simulation of 
the three-dimensional system of equations i70\ and (I74> . In 
Fig. [2]the collapse of the initially weakly modulated beam 
Ep = (A2(?oeo)"^/^A[l + 0.1sin(407rfcoz)]exp(-r2_/2a2) 
can be seen. Using A — \/0.02 and a — lOk^^, the beam 
interact with the initially homogenous radiation gas back- 
ground. The collapse is seen by the decrease in the beam 
width rj_ 



y2 and the increase in the beam energy 



density. In Fig.^]the evolution of an initially spherical pulse 

Ep = (A^(obeo)^^^^^exp(— r^/2a^) is seen. As the pulse 
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propagated through the initially homogeneous radiation gas, 
collapse ensues, as can be seen by the decrease in the pulse 
width r = ^ x"^ + y"^ + z'^ and the increase in the pulse en- 
ergy density. The values of A and a are the same as in the 
beam case. We note that the pulse in the last panel under- 
goes splitting, and the radiation gas response develops a pho- 
ton wedge, analogous to a Mach cone, through which energy 
is radiated. 



Growth ratev/ck^. A=0.02, p=0.01 




FIG. 13 The normalized growth rate 7/(cfco), as given through the 
dispersion relation <80> . as a function of K±/ko and K^/ko- We 
note that due to cylindrical symmetry, the area of non-zero growth 
rate is really a cone-l ike structure. Here A = X^SoeoEo (reprinted 
with permission from iShukla et alj <2004ct) . Copyright 2004, Amer- 
ican Institute of Physics). 



5. Pulse collapse and photonic wedges 

The approximate analytical results in conjunction with the 
numerical simulations presented above, shows that collapse, 
filamentation and pulse splitting as generic features of the sys- 
tem ilOi and ( 174^ On the other hand, both these methods of 
investigation have shortcomings, and it is therefore of great 
interest to obtain more analytical prediction valid for a wider 
range of parameters. In this section we will show that under 
very general conditions, pulse split and collapse are essential 
features of the photon pulse-acoustic system (I70> and i74l . 

As we have seen, the dynamics of an intense short photon 
pulse and the radiation background response is described by 
the system of equations (|70j and (|74}- 

Following Ishukla et al. ( 2004c), we re-normalize the sys- 
tem of equations (I70> and il4l . by introducing the new vari- 
ables E = {'ikofiXeoS'o/vgy^'^ Ep{t,rj_,z - Vgt), r = 
Vgt/ko, and SS' = {vg/2kon) {\E\'^ - It is implied 

that the wave packet has the cylindrical symmetry and moves 
along the z-axis. From Eqs. (I70> and (I74> we then obtain 

BE 

2i—+ylE-P,—E+\EfE + .yrE = 0, (83a) 
Bt Bz'^ 

and 

(83b) 
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FIG. 14 The three-dimensional evolution of a normalized modu- 
lated beam energy density X^(ffoeo\Epf (left panels) and the nor- 
malized radiation gas energy density A SS (right panels). The color 
gives the intensity. The beam collapses after a finite time, as does 
the inhomogeneity in the radiation gas (reprinted with permission 
from lShukla et all fM)4^ . Copyright 2004, American Institute of 
Physics). 



respectively. 

Representing SS' in terms of \E\'^ and J^, the part of the 
perturbation of the radiation gas that is concentrated in the re- 
gion of localization of the wave packet and takes part in its 
compression, has been isolated. The second part of the rep- 
resentation of SS corresponds to the interaction of the pulse 
with the radiation field This interaction is described by 
the last term in Eq. ( I83al i. As seen by Eq. (|83b} [or Eq. i74l]. 
the velocity Ug w c of the source (the pulse) exceeds the ve- 
locity c/ \/3 of the radiation waves, i.e. Vg > c/a/3, and it is 
therefore expected that the radiation field n will be localized 
behind the pulse. Hence, it is assumed that the pulse and the 
radiation field are localized in different volumes, and in the 
region with a possible overlap the relation ^ ^ 1-^1^ is sat- 
isfied (see also the end of the next section). This inequality 
allows us to neglect the last term in Eq. ( I83at . which means 
neglecting the back reaction of the radiation on the pulse. In 
this approximation, the pulse field E drives the radiation field 
jy through Eq. ( l83bL while the pulse propagation is unaf- 
fected by the radiation field. 

From Eq. (|83aJ we have the conservation of the "field mass" 
parameter of the pulse 




(84a) 
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FIG. 15 The evolution of a normalized spherical pulse energy den- 
sity \^SQeo\Ep\'^ (left panels) and the normalized radiation gas en- 
ergy density A (right panels). The pulse first undergoes compres- 
sion, and later splitting, while the radiation gas develops a high in- 
tensity region, from which photon wedges (similar to a Mach cone) 
are radiated (reprinted with permission from .Shukla et alj (i2004c.') . 
Copyright 2004, American Institute of Physics). 



and the Hamiltonian 



(84b) 



where the Hamiltonian density is given by J^f — jVj^Sp — 
(3,\d,E\^ - \E\y2. 

Some information on the spatio-temporal behaviour of the 
initially localized wave-packet can be found by following the 
evolution of the characteristic widths of the packet in the 
transversal and longitudinal directions. These widths are de- 
fined by 
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N 



and 



N 



rl\E\'^ dr^dz, 



z^\Efdr±dz. 



(85) 
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Straightforward calculations give the following evolution 
equations 
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N 
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dr±dz, 



(87) 




i dr^dz > (88) 



for the widths (I85t and ( I86t . respectively. 

In the two-dimensional (2D) case, the process of self- 
compression of the wave-packe t can be clearly seen 
(Kuznetsov et al., ^99?;' Tskhakava, 1982). A necessary con- 
dition for collapse in this case consists in the predominance of 
the "field mass" N of the given state over some critical value 
Nc, i.e. N > Nc, where Nc is the "field mass" calculated 
for the stationary ground state (Kuznetsov et al., 1995). The 
ground state is described by the (positive) radially-symmetric 
solution 'i> = E{r, r) exp(— iAt) of the equation 



0, 



(89) 



derived from Eq. j83a> and computed with A = 1. In the 2D 
case, the right-hand side of the expression for the mean square 
transverse radius ( I87t is 2H/N, i.e. a combination of the con- 
served quantities. Thus, integrating Eq. ( I87> twice, we obtain 



(90) 



R^{t) ^ R^{Q) + Ct+ (H/N)i 



where the constants i?^(0) and C are defined by the initial 
conditions. Hence, a sufficient condition for 2D self-focusing, 
ultimately leading to complete collapse in a finite time, is H < 
0. Because this is independent of the value of C, there will 
always exists a finite time to for which the transverse radius 
vanishes. In the three-dimensional (3D) case and for f3z < 
(corresponding to anomalous dispersion) the equality ( I90I I 
must be replaced by the inequality 



R^ir) <Co + CiT + iH/Ny 



(91) 



and the sufficient co ndition for the collap se of the wave-packet 
is again H < (K uz netsov et all|l995l) . 

Equation ( 18 3 at with > and ^ = corres ponds to the 
normal dispersion region of the wave -pac ket. InlBerg e et alJ 
( ll995l) : lBer'ge and Rasmusseiil (fl996) and lBeree et alJ ( 1^98^ 
the features of the pulse self-focusing, described by the so- 
lution to Eq. ( I83a> . have been investigated in detail. Their 
conclusions, can be applied directly to the numerical re- 
sults pre sented here concer ning the break-up of the wave- 
packet. iBer ge et alJ (Il995h showed that the characteristic 
length Z'^{t) of the wave-packet localization along the di- 
rection of propagation is bounded from below by a positive 
constant. The relation ( I88t indicates an asymptotic longitudi- 
nal spreading. Thus, a 3D (global) collapse cannot take place 
in media with normal dispersion, since the necessary condi- 
tions for the pulse self-focusing to occur is (dR^/dr) < and 
(dZ^/dr) > 0, which explicitly excludes the 3D case. The 
self-focusing in the transverse direction is accompanied by a 
longitudinal spreading and in the following by the splitting of 
the pulse. This does not, however, immediately exclude par- 
tial (local) collapse scenarios. 

The process of pulse splitting close to the time of self- 
focusing t is described by Berge and Rasmussen ( 1996); 
iBerge et aO ( Il998h employing a quasi-self-similar analysis. 
This approach uses, besides the description of the solution 
of Eq. ( I83a> (neglecting ^), the time-dependent character- 
istic lengths in the longitudinal l^jt) and trans versal I \ (t) 
directions. iBerge and RasmussenI d 19961) and iBerge et all 
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1119981) found that the transverse scale l± (t) exhibits a chang- 
ing behaviour as one passes some critical point z*{t) on 
the z-axis. Inside the localized region z < z*{t), the 
transverse width collapse with the rate l±{t) ^ {to — 
t)i/2/[ln{ln [l/(to ~ ^)]}]^/^, while in the complementary 
(delocalizing) domain, z > z*{t), the transverse pulse width 
spreads out with the rate l±{t) ~ •\/<ln{ln[l/(to — i)]}- 
Hence, the time derivative of l±{t) changes sign around the 
point z*{t). Meanwhile, the self-similar longitudinal scale 
Iz (tHncreases slowly in time. iBerge and R asmussenI ( Il996h 
and^erge et al. ( 1998) found that near the self-focusing time 
t to, lz{t) ~ t, which implies a linear increase of z*{t) 
in time such that z*{to) (~ ^/Plhito))- The presence of the 
coefficient ^/J3^ leads to the decrease of the distance of the 
critical point z*{to) from the origin for small P^. The scale 
l± (t) remains strictly positive at times t ^ tg. Consequently, 
the transversal scale reaches a minimum value at a finite dis- 
tance z*{to). 

Since the wave-packet is assumed to be cylindrically sym- 
metric and also symmetric relative to the origin z = 0, the 
total field distribution during self-focusing must exhibit two 
maxima located at z — ±z*{to), respectively. The wave- 
packet has therefore been split into two identical smaller cells, 
symmetrically placed on each side of the origin z = 0. 

Hence, a wave-packet, propagating in media with 
anisotropic dispersion, will be spread out along the direc- 
tion of the negative dispersion, and split up into smaller cells. 
These a nalytic results have been confirmed by numerical so- 
lutions (E^gEet^ [l998'; Tutheret^ The duration 
At of self-focusing, accompanied by the pulse splitting, can 
be estimated using Eq. j87> as 



At 



(92) 



In the first part of our numerical solution — where the split- 
ting of the wave-packet takes place — coincides with these 
results. The small coefficient which determines the neg- 
ative dispersion, changes the distance (from the origin) along 
th e z-axis, at which the field is l ocali zed after spl i tting. 

iBerge and RasmussenI lll996h and lBerge et alJ lll998h have 
shown that this splitting process can be continued (multi- 
splitting process) if the newly formed cells will possess a 
transverse energy higher than the self-focusing threshold Nc- 
In our case, the wave-packet splits into two cells only, as the 
energy localized in each new cell is below the threshold Nc- 
Furthermore, the wave-packet also loses energy to the radi- 
ation gas during the splitting process. The analysis of the 
formation of these photonic wedges, or the radiation Mach 
cone, can b e analyzed in accordance with the presentation in 
IShukla eTal. ( ,2QQ4 c), such that e.g. the energy loss from the 
pulse can be estimated. 

In conjunction with pulse collapse in a radiation gas it 
should be mentioned that if the field invariant — c^B^ > 
pair creation will occur as the pulse intensity increases, and 
the loss of energy through the photonic wedges will be neg- 
ligible in comparison to the energy radiated into Fermionic 
degrees of freedom. This will give rise to a rich and complex 
dynamical interplay between the pulse photons, the radiation 



gas, and the pair plasma (Bula nov et '2005^. For the case 
of interaction between a pulse and a pure radiation gas, we 
have c|B| > |E|, due to zero dispersion, and pair creation 
would not occur However, since we will in practice always 
have some ionized particles present, pair creation is likely to 
be the result of pulse collapse due to weak dispersive plasma 
effects. 



6. The strong field case 

Our knowledge of the nonlinear refractive properties of 
the radiation gas gives a means for investigating the ef- 
fects of higher-order nonlinear corrections to the standard 
first order Heisenberg-Euler Lagrangian, and to probe the 
significance of higher-order effects for photonic collaps e 
llMarklund et all I2003L l20Q4at LShukla and EUassoni l2004ah . 
The dynamics of coherent photons, travelling through an in- 
tense radiation gas, may be analysed as above, following 
iMarklund et alJ ( l2003ft . We obtain a nonlinear Schrodinger 
equation fo r the slowly varying pulse envelope Ep according 
to (iKivshar and Agrawa li l2003h 
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dt 



vo • V re. 
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(93) 

where the subscript denotes the equilibrium background 
- (ko ■ V)^ J(? = - ifo is a perturbation. 



State, 



vq = v{S'o), HQ = n{S'o), n„i((5^) = I]m=i "o SS'''^/m\, 
nQ™"* = (T^no/dS^, and the refractive index n is given 
through Eqs. j37> and j38l l. 

For a dispersion relation uj — |k|c/n(r, f), the motion of 
a single photon may be described by the Hamiltonian ray 
equations ( 17 U ( Mendonc^i 1200 1|). Since n = n{<S) and 
dn/dS > always hold (see 'M arklund et alJ ( l2005t JI)'). a 
denser region of the radiation gas will exercise an attractive 
force on the photon jPartovi 1 19941) . thus creating lensing ef- 
fects. The single particle dynamics thus supports that photonic 
self-compression is an inherent property of the one-loop radi- 
ation gas, but we note that as the density of a region increases, 
the phase velocity approaches a constant value, given by ( 14 U . 
i.e. V Inn ^ 0. 

Following Marklund et al. ( 2003), the response of the radi- 
ation gas to a plane wave pulse may be formulated in terms of 
an acoustic wave equati on, generalizing Eq. (|74l l to the strong 
field case, according to dMarklund et alil2005bi) 
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If the time response of is slow, Eq. ( I94t gives 
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Using ( I95> and the expression for n^\{5<S), we can write Eq 
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When n'^\E, 



^/2nQ ^ 1, we have a self-focusing nonhnear- 
ity in Eq. ( I96> . but as | £'p | grows the character of the nonhnear 
coefficient changes. The coefficient is positi ve when l-EnP < 
El, = \2n'„/n'^^\, but since nf,' < for all S'o fMarklu nd etall 
^005b), the sign changes as the pulse amplitudes grow above 
the saturation field strength E'sat, making the nonlinearity de- 
focusing and arresting the collapse. The numerical value of 
this turning point is dependent on the background parameter 
(fo- For low intensity radiation gases, « 0, and Eq. i9bi al- 
ways displays self-focusing, i.e. the field strengths can reach 
values above the Schwinger field. When <oo roughly reaches 
the critical value S'^nt, the weak field approximation breaks 
down, and Eqs. ( I37> and Hli can be used to derive an expres- 
sion for Esat- As an example displaying the general character 
of the intense background case, consider (oq = ^fcrit x 10^. We 
find that _E,sat ~ 2 x 10^'' V/cm > Ecm, i.e. the pulse satu- 
rates above the Schwinger critical field. Thus, both the weak 
and moderately strong intensity cases, the latter desribed here 
by Eq. ( I96> . display self-compression above the Schwinger 
critical field. 

This analysis can be generalized to take into account the 
statistical spread in th e coherent pulse, g iving rise to a damp- 
ing of the instabilities llMarklundll2004 . 



7. Other field configurations 

Above, we have seen that the propagation of an electro- 
magnetic pulse through a radiation field gives rise to instabil- 
ities, wave collapse, self-focusing, and pulse splitting. These 
concepts can be carried over to the case of multiple beams or 
pulses propagating through a radiation gas, including pulse in- 
coherence. Marklund et al. (2005c) showed that when several 
pulses are present, they can exchange energy via a background 
radiation gas and instabilities can occur, even if their propaga- 
tion is p arallel. 

IShuklaet alJj2004dl) considered an incoherent non-thermal 
high-frequency spectrum of photons. As will be shown, this 
spectrum can interact with low-frequency acoustic-like per- 
turbations. The high-frequency part is treated by means of 
a wave kinetic description, whereas the low-frequency part 
is described by an ac oustic wave equation with a driver 
(iMarklund etM l2003h which follows from a radiation fluid 
description. The high-frequency ph otons drive low-frequen cy 
acoustic perturbations according to ( IMarklund et allEoOSl) 
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^ + c^V^^ J ha;f{k,r,t)d'k, (97) 



where the constant Sq is the background radiation fluid en- 
ergy density and / is the high frequency photon distribution 
function. This hybrid description, where the high-frequency 
part is treated kinetically, and the low-frequency part is de- 
scribed within a fluid theory, applies when the mean-free 
path between photon-photon collisions is shorter than the 
wavelengths of the low-frequency perturbations. We note 
that the specific intensity Ik — hojf/eo satisfies Eq. i72\ . 
and is normalized such that = J Ikd'^k, where E 

is the high-frequency electric field strength. These equa- 
tions resemble the photo n-electron system in the paper by 
IShukla and Stenfld(ll998h . where the interaction between ran- 
domly phased photons and sound waves in an electron- 
positron plasma has been investigated. 

Next, we consider a small low-frequency long wavelength 
perturbation of a homogeneous background spectrum, i.e. 
/ = /o + fiexp[iiKz - m)], l/il « /o and £ = 
Si cxp[i(ifz — fit)] and linearize our equations. Thus, we 
obtain the nonlinear dispersion relation 



1 = - 



where — jX^c^HSq. 

(a) For a mono-energetic high frequency background, we 
have /o — no'5(k — kp). The nonlinear dispersion relation 
( I98> then reduces to 



(99) 



n-" + K^c^)[Kc + {2fi - ■iKccosOo) cos 6*0], 

where we have introduced cos 6*0 = ko ■ z- This mono- 
energetic background has a transverse instability when Oq — 
7r/2, with the growth rate 



r = 



Kc 



-,1/2 



- 1 



(100) 

where 7^ = —ifl, and vt = {iinokocY^^ is a characteristic 
speed of the system. The expression in the square bracket is 
positive definite. 

In fact, when S'ohckono / E*^^ ^ 1, a condition which is 
satisfied due to (|4}, we have vt ^ c Using the expression 
i99\ . we then have two branches. The branch corresponding 
to f2 « Kc/ V3 is always stable for small vt, while for the 
branch corresponding to J? w Kc cos we obtain the growth 
rate 



r = KvT 



1 — cos fc'o 
1 — 3 cos Oa 



(101) 



In Fig. 



which is consistent with ( I100> in the limit 6q ti/2. 
[T6l the behavior of the growth rate dlOU is depicted. 

(b) The high-frequency photons have generally a spread in 
momentum space. For simplicity, we here choose the back- 
ground intensity distribution as a shifted Gaussian, i.e. 



^0 



^3/2 fc3^ 
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(102) 
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FIG. 16 {F/KvtY, according to Eq. jlOH . plotted as a func- 
tion of qq = cos 00 in the mono-energetic case. (Reprinted from 
IShukla etalJj2004j) . Copyright (2004), with permission from Else- 
vier.) 



where J^o = (l-^oP) is the (constant) background intensity 
and kw is the width of the distribution around ko. Assuming 
that the deviation of ko from the z-axis is small, and that 5 = 
ko/kw ^ 1, we can integrate Eq. ( 19 8> with ( I102> . keeping 
terms linear in S, to obtain 



6 cos 6*0 



772 - 1/3 



8Sr] cos 60 



(103) 



J] — ASt] cos 6*0 I (2 arctanhry — in) 



for < rj < 1. Here rj = Q/Kc and — 
(4/97r3/2)A2eo4J^oexp(-A:g/fc^). Thus, we see that the 
non-zero width of the distribution complicates the characteris- 
tic behavior of the dispersion relation considerably. It is clear 
that the width will reduce of the growth rate compared to the 
mono-energetic case. 

We may also look at the case when the time-dependence 
is weak, i.e. d'^S/dt^ < c^V^^, such that Eq. (E) yields 
S = 2A(ob / hiokf <Pk. Upon using this relation, we find that 
\7uj = -iik\7 J k'f d^k'. Hence, Eq. ^ becomes 
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dt 
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= 0. (104) 



A similar equation may be derived for the specific intensity 
Ik . Equation ( I104> describes the evolution of high-frequency 
photons on a slowly varying background radiation fluid, and 
it may be used to analyze the long term behavior of amplitude 
modulated intense short incoherent laser pulses. The results 
in this section can be generaliz ed to several partially coherent 
pulses ( iMarklund et alll2005d) . 



C. Effects due to plasmas 

Plasma channels are closely connected to both the plasma 
dispersion and the pro pagation of wave mo des in waveguides. 
IShen et alJ (l2003h and lshen and^ ( l2003l) first suggested the 



use of plasma cavitation as a means of fostering conditions 
in which nonlinear quantum vacuum effects, such as photon- 
photon scattering, could take place. As a high intensity elec- 
tromagnetic pulse propagates through a plasma, the interac- 
tion with the plasma may completely evacuate regions giving 
conditions similar to the that of Sees. IIII.A.3l and lIII.A.4l 

Moreover, although in many cases the presence of a plasma 
will swamp the effects due to photon-photon scattering, it can 
under certain circumstances provide a means for the propa- 
gation of non-classical plasma modes. Due to the nontrivial 
dispersion of electromagnetic waves in plasmas, there will 
be a net effect due to photon-photon scattering, such that 
low-frequency modes will be generated. In general, the ef- 
fects of a nonlinear quantum vacuum i s expected to become 
pronounced for next gene ration lasers jBulanov et all I2OOOI 
12004: Mourou et alil2005i) 

Also, for future applications, the combination of photon- 
photon scattering induced pulse compression in conjunction 
with pair creation ( Nitta et al.,, 2(304) could provide interest- 
ing insights both into fundamental properties of the quantum 
vacuum as well as into the prospects of creating high power 
electromagnetic sources. The effects of plasmas within the 
environment of a quantum vacuum therefore deserve further 
investigations. 



1 . Plasma cavitation and plasma ctiannels 

If the power of the laser pulse propagating through the 
plasma surpasses the critical value Pcrit = 17 {lu / ujp)'^ GW , 
where u; is the laser frequency and ojp the electron plasma fre- 
quency, there may be complet e expulsion of p lasma particles 
from the high intensity region ( Ma x et allll974 ). thus forming 
a wave-guide (Shen and Yu, 2003). In such wave-guides, the 
effects of photon- photon scattering could be of importance 
dShen et 'ail l2003h 
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see, e.g. 


Bulanov et al. 


iGoloviznin and Shet) 


IShukla et alJ 


2004bh: 



Taiima and Taniuti' ( 119901) . and iMourou et alJ (|2P05 
Pukhov (2003) for recent reviews). In the context of doing 
fundamental physics and mimicking astrophysical events in 
laboratory environments, the evolution of laser intensities 
has received a lot of attention. Examples of experimental 
suggestions are axion detection (as a dark matter ca n didate ) 
dBernardl Il999t iBradlev et all l2003l iDupavs et a ll. '2005|); 
pair production (see , e.g. jRingwal dl (l2001alb[ 12 003). and 
iBamber et all d 1 999') : 'Burke et al .' (' 1 997) : "Meverhofer { 199% 
for a discussion of the detection of pair production from real 
photons); laboratory calibration of observations, relativistic 
jets, analogue general relativistic event horizon experiments 
(su ch as Hawking and Unruh radiation (Hawking, 
lUnruh. 1976)), a nd probing the q uantum spacetime properties 
IChea .2003l .Chen and Taiimal fl999.) . The possibility of 
reaching extreme power levels with such setups is one of 
the p romising aspects of laser-plasma systems (Bi ngharnI 
l2003h . and also holds the potential of overcoming the laser 
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intensity limit 10^^ W/cm^ (iMourou et all Il998l) . As 
the field strength approaches the cr itical Schwinger field 
Edit ^ 10^'' V/cm (Schwinge3 ll95ll) . there is p ossibility of 
photon-photon scattering, even within a plasma ('Sh en et all 
12003) . as the ponderomotive force due to the inten se laser 
pulse gives rise to plasma channels llYu etallflQSl . Under 
such extreme circumstances, the effects of pair creation will 
be pronounced. Electron-positron plasmas are also produced 
by interactions of matt er with powerfu l multi-terawatt an d 
petawatt laser pulses tGahn et all l2000t iLiang et all Il998h . 
The concept of trident pair-production, as described within 
the framework of perturbation theory, could give a means 
for creat ing electron-posi t ron pa irs by intense laser pulses in 
vacuum llBerezhiani et al.L Il992h. Moreover, the future x-ra' 



•ay 

i 



free electron laser systems (iPatell '2004: Ringwald, 2001a' 
Eb03 ) coul d result in methods for creating pair plasmas in the 
laboratory jAlkhoferetalll200ll) . The possible field s trength 
output could reach E « 0.1 ('Alkh ofereta'ill200l() . Even 
on an experimental level, pair-production due to collisions 
of electron backscattere d photons w i th the original photon 
beam has been observed jBurkeet allll99'1) . Thus, there are 
ample evidence that the investigation of nonlinear interactions 
of pair-plasmas a nd high intensity e lectroma gnetic fields 
deserves attention OFarina and Bulanovt i2001ai: iKozlov et all 
Il979ah . 



a. The effect of re lativistic nonlinearities Here we follow 
IShukla et al.1 ll200-5h . The propagation of a circularly polar- 
ized inte nse laser pulse in an unmagnetized plasma is gov- 
erned by jYuetal.L n"982) {df - c^\/^)A + {ljIN/j)A = 0, 
where A is the vector potential of the laser pulse, ojp = 
(7ioe^/eor7ie)^/^ is the unperturbed electron plasma fre- 
quency, 7 = + e^jAp/m^c^ is the relativistic gamma 
factor including the electron mass variation in intense laser 
fields, and N = rie/np is the ratio of the electron number 
density to the background plasma number density uq. 

At intensities beyond lO^'* W/cm^, the electron quiver 
speed Vosc = 6 x 10^^^c\\/l exceeds the speed of light, 
and hence nonlinear effects in plasmas cannot be ignored. 
Here / is the intensity in W/crn^ and A is the laser wave- 
length in microns. T hus, the relativi stic ponderomotive force 
Cshuklaetal., 1986: lYu et allll982ft F = -meC^V-f of in- 
tense laser pulses will separate charges and thereby would 
create a huge ambipolar electric potential (j) in the plasma. 
At equilibrium, the balance between the relativistic pondero- 
motive force and a slow electric force eV0 will yield (/> = 
{meC^ /e){'-f — 1), which, when substituted into Poisson's 
equation, gives = 1 + AgV^7. Here Ag — c/ujp is the elec- 
tron skin depth, and the ions are assumed to be immobile. The 
electron density will be locally evacuated by the relativistic 
ponderomotive force of ultra-intense nonuniform laser fields. 
The laser pulse localization and compression would then oc- 
cur due to nonlinearities associated with relativistic laser pon- 
deromotive force created electron density evacuation and rel- 
ativistic electron mass increase in the laser fields. This phe- 



nomena can be studied by means of the equation 
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(105) 

where A — eA/mgC, and the time and space variables are 
in units of ujp^ and Ag, respectively. For the propagation 
of a modulated laser pulse along the z axis, we obtain from 
Eq. (I105t after invoking the slowing varying envelope approx- 
imation. 



dt 



dz 



2iuj{—+Vg—\+^'l + I--{l + V'P) = 0,(106) 



where we have set A = (l/2)/(r, z, t){x + iy) exp(— iwi + 
ikz)+ complex conjugate, and denoted P = (1 + /^)^/^. 
Here the normalized laser frequency and the normalized laser 



group velocity are denoted by = + and ^ 



k/uj. 



respectively. In the one-dimensional case [viz. set = 
d^/dz"^ in Eq. ( I106H . we have the localization of intense 
electromagnetic waves in t he form of a lar ge amplitude one- 
dimensional bright soliton llYuetallfl98l . We have numer- 
ically solved Eq. ( I106t in order to study the evolution of a 
cylindrically symmetric modulated laser pulse. The results 
are displayed in Fig. 1171 
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FIG. 17 The variation of /, as given by Eq. <106> . and against 
r (the radial coordinate) and z for an initial laser pulse which ini- 
tially has a Gaussian shape. We observe the pulse compression and 
the formation of a light bullet. (Reprinted with permission from 
IShukla etalJj2005h .) 

Initially, the pulse is assumed to have a Gaussian shape, 
/ = /oexp[— (z^ + r^)/2GG], and we used the normalized 
wavenumber k = 2 and initial pulse amplitude Iq — 1. We 
notice that the compression of the pulse envelope (left pan- 
els), which is correlated with the excavation of the normalized 
electron density (right panels). Our numerical results reveal 
that self-compression of the pulse is more rapid when one ac- 
counts for the relativistic light ponderomotive force induced 
electron density depletion, contrary to the constant density 
case (viz. N — 1). Physically, the enhanced compression 
and self-focusing of an intense laser pulse occur due to the 
localization of light in a self-created electron density cavity. 

The fact that evacuation of the plasma takes place as the 
light intensification due to self-compression occurs means that 
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the situation discussed in lShen et a l.' '("2003'), where the quan- 
tum electrodynamical effect of photon-photon scattering at 
high intensities takes place, could be realized in the next gen- 
eration laser-plasma systems. Moreover, as the intensities in 
the evacuated region increase, the concept of vacuum catas- 
trophic collapse, at which the pulse due to quantum vacuum 
nonlinearities self-compresses, may ensue ( iMarklund et all 
5004a). The intensities that can be reached at this stage, in 
principle, surpass the Schwinger field, but then the process 
of pair creation has to be investigated and removed, since 
this would otherwise quickly dissipate the electromagnetic en- 
ergy into Fermionic degrees of freedom. The problem of self- 
consistent analysis of pair pro duction in a plas ma environment 
has been approached by Bula nov et al.l ( l2005l) where a model 
for incorporating a particle source term was given (see also 
Eq. lfT36li 1. 



b. Self-interaction in electron-positron plasmas Consider the 
propagation of intense light in an electron-positron plasma. 
By averaging the inertialess equations of motion for electrons 
and positrons over one electromagnetic wave period, the ex- 
pressions for the electron and positron number densities 
and Up, respectively, in the presence of relativistic pondero- 
motive force (.Shukla et al., 1986) of an arbitrary large ampli- 
tude laser pulse takes the form 



and V^<I> = Ne — Np, respectively, where Ng = exp[/?e(l — 
VTTpf + $)], Np = exp[/3p(l - v/TTpf ~ $)], 
Pe^p = c^/i'Te.p' ™d VTe,p = {ksTe^p /nie)^/'^ . Here the 
dispersion relation ujq = c^k^ + Ld'^{neo + npQ)/nQ has been 
used, and Vg = (c^/ijJo)k is the group velocity. In the quasi- 
neutral limit Ne = Np, w e hav e $ = (1 - ^1 + |^p)(/3p - 
Pe)/if3p + Pe)- Equation ( II 10> then becomes 
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(111) 

where (3 = 2f]ePp/{Pe + (3p) is the temperature parameter. 

The dispersion relation for the modulational and filamen- 
tational instabilities for an arbitrary large amplitude electro- 
magnetic pump can be derived fro m il 1 1 ) following stan- 
dard techniques dShuklaet al.L 119871 il98 8). From the ansatz 
A ^ (ao + ai)exp(i5r), where oq is real, gq 3> |ai| and 
(5 is a constant nonlinear frequency shift, the lowest order 

solution gives ^ = 1 — exp /3 ^1 — ^1 + a^j + Gq. 

Linearizing Eq. il 1 1> with respect to ai, with the ansatz 
fli = (X + iY) exp(iK • ^ — iSlr), where X and Y are 
real constants and SI (K) is the frequency (wavevector) of the 
low-frequency (in comparison with the light frequency) mod- 
ulations, the nonlinear dispersion relation reads 



no exp 
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where ^ and qp — e are the electron and positron 
charges, /cs is Boltzmann's constant, Te{Tp) is the electron 
(positron) temperature, and 7 — + e^\h.\^ j m\(? for cir- 
cularly polarized light. The ambipolar potential associated 
with the plasma slow motion is found from Poisson's equation 
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which gives the modulational instability growth rate F 
—tfl according to 
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From the particle momentum conservation equations, the 
electron and positron velocities are given by Vg p — 

-qe^pA/me^cp, where 7e,p = (1 - We,p/c^)~^^^ = 7- 

The dynamics of the intense light is obtained from 
Maxwell's equations and reads 



7"o 



Up) A = 0, 



(109) 



where A is the vector potential. 

For a circularly polarized electromagnetic wave A = 
(l/2)A(r, i)(x + iy)exp(ik • r — iujot), using the scalings 

ne,p = noNe,p, t = TUJo/ujp, r = - UgT)/wp, Ug 

{u!o/ujp)\-g/c, A ~ {mec/e)A, and = {meC^ /e)^, Eqs. 
(I108> and ( I109> can be written in the dimensionless form as 



,0A^1 f jV^ 

dr 2 I 



(113) 

The growth rate increases with the larger (3 values (i.e. for 
lower temperature), while for the intensity field, we do not 
necessarily obtain higher growth rates for higher intensities, 
see Fig. ^] This is attributed to an interplay between the rel- 
ativistic particle mass variation and the relativistic light pon- 
deromotive driven density responses. From the expression 
M\3\ one observes a decrease in the growth rate for large 
enough (3. However, this result should be interpreted with 
caution since for large-amplitude fields in a low-temperature 
plasma, electron inertia effects will become important and 
may dominate over the thermal effects. Then, in this case the 
assumption that the electrons (and positrons) obey a modified 
Boltzmann distribution may no longer be valid. 

Since N^ = Np = cxp[/?(l - ^1 + \A\'^)], the increase 
of the pulse intensity will cause an almost complete expul- 
sion of the electrons and positrons from that region, see Figs. 
[T9I and |20l The simulations show the evolution of an ini- 
tially weakly modulated beam ^ = 10"^[l + 0.Q2sin(z/8) + 
A = 0, (110) 0.02 cos(z/4) + 0.02 cos(3z/8)]exp(-rV32), where = 
+ y^, while the electron perturbation is zero initially. Thus, 
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as can be seen in Figs. 1191 and 1201 the modulated beam self- 
compresses and breaks up into localized filaments. This gives 
rise to electron and positron holes, and as the pulse intensity 
grows, the conditions for pure elastic photon-photon scatter- 
ing improve within these holes. 
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FIG. 20 The normalized electron number density = exp[/3(l — 
\/l + \A\^)] as a function of r and z at five different times r (A is 
given by Eq. IIIH '). Note the efficient expulsion of plasma particles 
in the later panels, an d their correlation to the intensity peaks in Fig. 
fT9l (Reprinted from lShukla et alJ <2004el) . Copyright (2004), with 
permission from Elsevier.) 



FIG. 18 The modulational instability growth rate F as given by Eq. 
<1 13> versus the wavenumber K, for (3 — 200 and ao —0.1 (dashed 
lines), P = 100 and ao = 0.1 (solid line), f3 = 100 and ao = 
0.2 (dash-dotted lines), an d for = 100 and ao = 0.05 (dotted 
line). (Reprinted from .Shukla et alj <2004d) . Copyright (2004), with 
permission from Elsevier.) 
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FIG. 19 The light amplitude A in logjQ-scale at five different times 
r (see Eq. <111» . The horizontal axis represents the distribution 
along the radial (r) coordinate, and the axial (z) distribution is on 
the vertical axis. (Reprinted from Shukla et al. l2004e) . Copyright 
(2004), with permission from Elsevier) 



c. Thin-foil amplification As noted in the previous two sec- 
tions, trapping and amplification of laser pulses can take place 
given the right plasma environment. This could be an impor- 
tant tool for stepping up the available electromagnetic inten- 
sities, and could therefore be important for investigations into 
photon-photon scattering. Here we will describe a method 
which could yield high intensity pulses. 

Laser-foil interactions have been use d as a method for pro- 
ton acceleration on table top scales fcMcKenna et all l2005t 



ISilva et all l2004l IZepf et all l2003l) . By letting a high inten- 
sity laser pulse impinge normally on a thin metal foil the foil 
material is ionized, creating a plasma in which protons are 
accelerated up to MeV-energies. The exact mechanism(s) be- 
hind the proton acceleration is still not completely clear al- 
though there exis ts a number of plausible suggestions (see 
ZepfetalJ (l2003h for a discu ssion). It was suggested by 



Shen and Mever-ter-Vehn (2001a) that this could be used to 



create confined high density relativistic electron plasmas. Let- 
ting two counterpropagating laser beams illuminate a thin foil 
normally, a spatially confined high density plasma could be 
created, and be used for, e.g. har monic generation, pair pro- 
duction , and 7 photon generation (|Sh en and JMever-ter-Vehnl 
l2001bh . 

Building on the w ork of lShen and Mever-ter-Vehi] ( 12001 ah . 
IShen and Yul ( l2002h suggested to let two oppositely directed 
laser beams interact via two closely placed thin foils. As 
above, when the high intensity lasers impinges normally on 
the thin foils, the foil mat erial will be ionized and a plasma 
will be produced. IShen and Yu ( 200Z) showed that this may 
le ad to electrornagneti c trapping. 

|Sh en and Yul ( l2002h started with the trapping of a circu- 
larly polarized electromagnetic pulse propagating in the z- 
direction in a positiv e electron density profile, using the sta- 
tionary equations of Shen and Me ver-ter-Vehnl (12001 al) and 
IShen and Mever-ter- Vehn (2001b) 



M = (7^ - l)ce'/uj, 



(114a) 



W=Ul^-iy^ {cj'/Lof + M^ +i7(7-2iV0,(114b) 



where M and W are two constants motion. 



7 



^/l + e^jylp/mgC^ is the relativistic gamma-factor, uj is the 
laser frequency, A is the vector potential cx exp[?a;t + iO{z)], 
Ni ~ ni/ric, n.i is the constant ion density, Uc = 1.1 x 
lO'^^ {\/ fivaY cm^^ is the critical electron density, A is the 
laser wavelength, and the prime denotes differentiation with 
respect to z. It is possible to find solitary solutions of Eq. 
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il 14> representing trapped electroma gnetic pulses between 
parallel high den sity plasma regions sEsirkepov et all Il998t 
iKim et all l2000l) . These analytical soliton solu tions suggest 
the possibility to trap laser light between foils. Shen and 
fc002 ) performed particle-in-cell (PIC) simulations of the sys- 
tem il 14t . Using the foil spacing A = 0.46A they showed 
that such configurations would yield a 100-fold amplification 
of the initial laser pulse intensities, over a trapping time of 26 
laser cycles, which is in the fs range for /im lasers. 

In multidimensional environments, true analytical soli- 
ton solutions are not known, but it is a well-established 
fact, due to approximate and numerical investigations 
llKivshar and Agrawall l2003h . that solitary like solutions ex- 
ists in dimensions > 2. However, these solutions are un- 
stable, and will suffer either attenuation or self-c ompression, 
depending on intensity and pulse width jPesaix et al.. 199ll) 
(see Fig. |9}. Thus, in a two-dimensional thin-foil environ- 
ment, light intensification could take place. The production of 
high-intensity pulses by these thin-foil amplification also has 
the valuable property of being realizable in a relatively small 
scale setting. 



d. Laser-plasmas and relativistic flying parabolic mirrors As we 
have seen above, the propagation of intense electromagnetic 
pulses in plasmas yields interesting nonlinear dynamics, and 
effects such as pulse self-compression can occur. These non- 
linear effects act as a very promising tool for, e.g. produc- 
ing intense ion beams teulanov et al.[ l2004l 'Esirkep ov et all 
[l999, 2 004}, which is of importance in laboratory astro- 
physics dChenl I2OO3I) . With regards to the nonlinear quan- 
tum vacuum, an interesting proposal has been put forward by 
iBulanov et al. (2003). The self-compression of laser pulses, 
towards intensities close to the Schwinger limit, can take place 
by using r elativis tic flying parabolic mirrors. 

Bulanov 'et alJ (|2P03) consider a plasma wakefield in the 
wave-breaking regime. They let a short intense laser pulse 
create a wakefield in a plasma, such that the wakefield phase 
velocity equals the laser pulse g roup velocity (which is clo se 
to c) in an underdense plasma (Taiima and Dawson. Il979l) .^ 
Due to the nonlinearity, the resulting wake field will experi- 
ence wave steepening entering the wavebreaking regime, to- 
gether with a local electron density spike approaching infinity. 
With such a set-up, a sufficiently weak counter-propagating 
laser pulse will be partially reflected from the electron den- 
sity maximum. The relativistic dependence of the Langmuir 
"mirror" on the driving laser pulse intensity will cause bend- 
ing of the surfa ces of constant phase, th us creating a parabolic 
plasma mirror (iBulanov and Sakharovltl99 1 ). This curvature 
of the plasma mirror will focus the counterpropagating (weak) 
pulse to a spot size A/47^ along the mirror paraboloid axis in 
the laboratory frame, where A is the wavelength of the source 



^ We note that the formation of sub-cycle intens e solitary waves could pen- 
etrate into a highly overdense plasma l.Shen et alll2004) . transferring en- 
ergy between low- and high-density regions of the plasma, which could be 
of importance in, e.g. laser fusion. 



of the reflected pulse and 7 is the relativistic gamma factor of 
the wakefield. Similarly, the f ocal spot w idth is A/ 27 in the 
transverse direction. With this. IBulanov et al. (2003) showed 
that the intensity gain will be roughly 64(Z?/A)^7"^, where D 
is the width of the pulse effectively reflected by the mirror. 
The focal spot intensity of the reflected pulse can then be es- 
timated to 




(115) 



where tOd is the frequency of the pulse driving the Langmuir 
wave, u) is the frequency of the source of the reflected pulse, 
and I is the i ntensity of the source of the reflected pulse. 
IBulanov et all ( I2b03i) gave the following example of light in- 
tensification through Eq. dl 151 . A 1 /im pulse generate the 
Langmuir mirror in a plasma where ~ 10^^ cm^^. The 
estimate 7 sa LOd/^p then gives 7 ^ 100. The pulse to be 
reflected is assumed to have / ~ lO^^W/cm^, and D = 
400 /im. Then /focal ^ lO^^W/cm^, to be compared with 
the critical intensity Ic = ceqE^^^^ « 3 x 10^^ W/cm^. The 
estimated focal intensity thus seems to reach the Schwinger 
limit. However, with the given value on / it is likely that the 
backreaction on the Langmuir wave has to be taken into ac- 
co unt, thus altering the estimate. 

Bulanov et all ( l2003h also presented numerical results us- 
ing a fully relativistic code, see Figs. |2] and |22] Using 
a three-dimensional PIC code, an intense laser pulse drives 
the Langmuir wave along the x-axis, while the counterprop- 
agating source pulse for the reflected wave has an intensity 
~ IQis W/cm^ in the /im wavelength range. The value of the 
of the source pulse intensity is chosen in order to avoid degra- 
dation of the Langmuir mirror Figure}^ shows the electron 
density profile of the Langmuir wave. Using a Langmuir laser 
driver with intensity 4 x 10^® W/cm^ and /im wavelength, 
the plasma waves move with with a phase velocity 0.87c, and 
the gamma-factor is 2. The density profile is shown along the 
X-axis, and a steep gradient can be seen. Figurel22l shows the 
electric field components of the source pulse and its reflec- 
tion [y — plane) and the Langmuir driver (z = plane). 
The focusing of the reflected can be seen. The intensity in- 
crease in the focal spot is 256 times the source intensity, i.e. 
/focal ^ 10^^ — 10^^ W/cm'^ for a /tm source laser. This is 
similar to the intensification obtained for the thin foil setup in 
Sec. lIII.C.Lcl 

e. Electromagnetic wave localization The nonlinear interaction 
of high-intensity ultrashort electromagnetic waves with hot 
plasmas is of primary interest for the fast ignitor concept of 
inertial confinement fusion and for the development of high 
power sources of hard EM radiation, as well as for laser- 
plasma particle and photon accelerators, and compact astro- 
physical objects containing intense electromagnetic bursts. 
Recent progress in the development of super strong elec- 
tromagnetic pulses with intensities / ^ 10^^-10^^ W/cm^ 
has also made it possible to create relativistic plasmas in 
the laboratory by a number of experimental techniques. At 
the focus of an ultraintense short electromagnetic pulse, the 
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FIG. 21 Typical electron density profile in the wake field, where 
the isosurfaces represents densities n — O.lSriciit, and ricit denotes 
the density at which the plasma goes from unde rdense to overdense. 
(Reprinted with permission from lBulanov et aljr2003i) .) 
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FIG. 22 The electric field components at d ifferent times. (Reprinted 
with permission from lBulanov et alJ j2003h .') 



electrons can acquire velocities close to the speed of light, 
opening the possibility of simulating in laboratory condi- 
tions, by using dimensionless simulation parameters, phe- 
nomena that belong to the astrophysic al realm. In the past, 
several authors presented theoretical ( Esirkepov et all 11998 ; 
Farina and Bulanov, 2001b; Kaw et al. ni992t iKozlov et al. . 
19791A and particle-in -cell simulation ( Bul anov et all [l999 ; 
Naumova et all l200ll) studies of intense electromagnetic en- 
velope solitons in a cold plasma, where the slow plasma re- 
sponse to the EM waves is modeled by the electron con- 
tinuity and relativistic momentum equations, supplemented 
by Poisson's equation. Assuming beam-like particle dis- 
tribution functions, relativistic electromagnetic solitons in 
a warm quasi-neutral electron-ion plasma have been inves- 
tigate d (Lontano et al., 2003). Experimental observations 
( iBorghesi et al.> .200 2) show bubble-like structures in proton 
images of laser-produced plasmas, which are interpreted as 
remnants of electromagnetic envelope solitons. 

IShukla and Eliassonl ( l2005l) presented fully relativistic non- 
linear theory and computer simulations for nonlinearly cou- 
pled intense localized circularly polariz ed EM waves and rel- 
ativistic electron hole (REH) structures jEliasson and ShuklaL 



l2006h in a relativistically hot electron plasma, by adopt- 
ing the Maxwell-Poisson-relativistic Vlasov system that ac- 
counts for relativistic electron mass increase in the electro- 
magnetic fields and relativistic radiation ponderomotive force 
jBingham et all i2004: Shukla et al. 1986). in addition to 
trapped electrons which support the driven REHs. Such a sce- 
nario of coupled intense EM waves and REHs is absent in any 
fluid treatment (E sirkepov et al., 1998; Farina and Bulano^ 
l2001ht iKaw et all f 19921 iKozlovetall ri979bl) of relativis- 
tic electromagnetic solitons in a plasma. Electromagnetic 
wave localization i s a topic of significant interest in photonics 
jMendoncaL 1200 ll) . as well as in compact astrophysical ob- 
jects, e.g. gamma-ray bursts (Pi rarL,2004) . 

The electromagnetic wave equation accounting for the rel- 
ativistic electron mass increase and the electron density modi- 
fication due to the radiation relativistic ponderomotive force 



jMendoncaL Eool F 



'dj/dz, where 7 = (1 



pi/mi 



e \A\ /m^c ) ' is the relativistic gamma fac- 



tor, are included. Here, pz is the z component of the electron 
momentum, A is the perpendicular (to z, where z is the unit 
vector along the z axis) component of the vector potential of 
the circularly polarized EM waves. The dynamics of nonlin- 
early coupled EM waves and REHs is governed by 
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and 
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fdpz - 1, 



(116) 



(117) 



(118) 



where A is normalized by rriec/e, 4> by ksTe/e, pz by 
rrieVTe and z by td- Here 7 = (1 + a^pj, + |Ap), Vte = 
{kBTe/meY, a = Vte/c, and TD = Vre/wp. In Eq. J116l l. 
we used the Coulomb gauge V • A = and excluded the lon- 
gitudinal (z-) component d^(t>/dtdz — jz, where jz is the par- 
allel current density, by noticing that this component is equiv - 
alent to Poisson's equation ( II 18> ('Shukla and Eliassorll2005h . 

Shukla and Eliasson have discussed the stationary as well as 
time dependent solutions of Eqs. dl 16> -( fTT8l in the form of 
REH which traps localized electromagnetic wave envelopes. 
Typical profiles for the amplitude of the localized EM vec- 
tor potential W and potential and density of the REH, as well 
as the local electron plasma frequency squared (fi^) includ- 
ing the relativistic electron mass increase, are depicted in Fig. 
I23I We observe that for large electromagnetic field the REH 
potential becomes larger and the REH wider, admitting larger 
eigenvalues A that are associated with the nonUnear frequency 
shift. This is due to the relativistic ponderomotive force of 
localized EM waves pushes the electrons away from the cen- 
ter of the REH, leading to an increase of the electrostatic 
potential and a widening of the REH. We see that the de- 
pletion of the electron density in the REH is only minimal, 
while the local electron plasma frequency il is strongly re- 
duced owing to the increased mass of the electrons that are 
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FIG. 23 Large-amplitude trapped EM wave envelope (upper panel), 
the potential (second panel), the electron number density (third 
panel), and the square of the local electron plasma frequency (lower 
panel) for large amplitude EMWwaves with a maximum amplitude 
of Wma = 1.5 (solid lines) and Wmax = 1.0 (dashed lines), and 
as a comparison a REH with small-amplitude EM waves which 
have Wraax 1 (dotted lines). The parameters are: the nor- 
malized speed Ho = 0.7, a = 0.4 and the trapping parame- 
ter 3^= —0.5 (correspondin g to a vortex distribution presented 
bv iBu iarbarua and Schamel ( 198l|) and ISchamel (2000) involving 
an eq uilibrium Synge-Jiittner distribution function ide Groot et aL 
11980)). The selected value of /3 are related to the maximum 
REH potential according to a specific relation similar to one in 
[Suiarbarua and Schamel (1981) and Schamel (2000). (Reprinted 
with permission from Shu kla and EliassonI <2005ll .) 
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FIG. 24 Phase space plots of the electron distribution function (left 
panels) and the modulus of the electromagnetic field (right panels) 
for f = 0, t = 50, t = 125 and t = 162. (Reprinted with permission 
from iShukla and Eliassort t2005»') .) 



accelerated by the REH potential; the maximum potential 
4'max ~ 15 in Fig. 1231 corresponds in physical units to a po- 
tential a^4>max X 0.5 X 10'' 1.2 x lO*" V, accelerating the 
electrons to gamma factors of w 6. 

In order to study the dynamics of interacting solitary struc- 
tures composed of localized REHs loaded with trapped EM 
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FIG. 25 The electromagnetic field (upper left panel), potential (up- 
per right panel), squared local plasma frequency (lower left panel) 
and electron density (lower rig ht panel) for two co l liding REHs. 
(Reprinted with permission from iShukla and EliassonI j2005|).) 

waves, IShukla and EliassonI ( l2005l) numerically solved the 
time-dependent, relativistic Vlasov equation ( II 17l i together 
with the nonlinear Schrodinger equation, that is deduced from 
Eq. dl 16> in the slowly varying envelope approximation. The 
results are displayed in Figs. 1241 and l25l As an initial condi- 
tion to the simulations, Shukla and Eliasson used solutions to 
the quasi-stationary equations described above, where the left 
REH initially has the speed = 0.7 (normalized by c) and is 
loaded with EM waves with W^ax = 1.5, while the right REH 
has the speed i>o = —0.3, and isjoaded with EM waves with 
W^max = 2.5. Furthermore. IShukla and EliassonI (12005) used 
fco = = in the initial condition for A and in the solution 
of the nonlinear Schrodinger equation ( Shukla and EliassoiJ 
120051) . Figure displays the phase space distribution of the 
electrons and the electromagnetic field amplitude at differ- 
ent times. We see that the REHs loaded with trapped EM 
waves collide, merge and then split into two REHs, while 
there are two strongly peaked EM wave envelopes at z w 30 
and z sa 70 remaining after the splitting of the REH. A popu- 
lation of electrons has also been accelerated to large energies, 
seen at z = 100 in the lower left panel of Fig. \^ The time 
development of the EM wave amplitudes, REH potential, the 
squared local plasma frequency and the electron number den- 
sity is shown in Fig. |25l Collision and splitting of the REHs 
can be observed, as well as creation of the two localized EM 
envelopes at z w 70; clearly visible in the left two panels at 
t > 150. 



2. Photon-photon scattering within plasmas 

a. Charged particle effects and Ch erenkov radia tion As pre- 
sented by jPremi iil (l2002h and Mar klund et al . (2005a), sim- 
ilar to a charged particle moving in a isotropic dielectric, a 
charged particle can suffer Cherenkov losses when propagat- 
ing through a intense gas of photons. The main difference 
compared to the case of a regular medium is the frequency 
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spectrum of the emitted radiation. Since the natural cut-off 
in the quantum vacuum is given by the Compton frequency, 
7-rays may be emitted by such a particle. 

In 1934, Cherenkov obse rved th e type of radiation now 
bearing his name (|CerenkoV|,ll934l) . His experimental result 
was explained bv lTamm and FranM ( ll937h . In an isotropic di- 
electric medium, a charged particle in rectilinear motion sat- 
isfying the so called Cherenkov condition, i.e. its velocity 
exceeds the (parallel) p hase speed i n the medium in which 
it moves, will radiate (IChefranovl 120041) . The radiation 
shock-front, called the Cherenkov cone, is analogous to the 
Mach cone formed as objects move with supersonic speeds 
through air. In quantum mechanical terms, the Cherenkov 
condition corresponds to energy and momentum conservation. 
Cherenkov radiation has technological uses, e.g. in determin- 
ing particle velocities. 

The dispersion relation for electromagnetic waves in 
an isotropic and homogeneous photon gas with refrac- 
tive index n is uj — kc/n, where n? = 1 -\- 5 and 
5 = 4A^/ 3 (Bialvnicka-Birula and Bialvnicki-Birulil 1 1 970t 
iMarklund et aL .200 3) (see Eg. (l29t). Thus, the refractive 
index in this case is always larger than one, and a particle 
may therefore have a speed u exceeding the phase velocity in 
the medium. The Cherenkov condition u > c/n for emis- 
sion of radiation can thus be satisfied. This condition can 
also be expressed in terms of the relativistic gamma factor 

= (1 — u^/c^)~^/^, namely 6^"^ > 1. We will here as- 
sume that a particle with charge Ze, satisfying the Cherenkov 
condition, moves through an equilibrium radiation gas. The 
energy loss at the frequency uj per unit length of the path of 
the charged particle is then 

-duj^ Vt^ ^hujduj, (119) 



ds 



c (7^-1) 



and the number of quanta N emitted per unit length along the 
particles path is 



dN , Z'^a (572 - 1) , 

——duj — — dto. 

ds c (7"^ — 1) 



(120) 



Since 5 is normally much less than one, we need a large 
gamma factor to satisfy the Cherenkov condition. Subse- 
quently, for 57^ = 1, we have 

U = NhjjJe, and N = Z^La5/\e, (121) 

respectively, where we have used the Compton frequency as a 
cut-off. Here L is the distance traveled by the charge. 

At the present time, the cosmic microwave background 
has an energy density of the order S' ^ 10^^^ J/m'^, i.e. 
5 ^ 10^''^, i.e. the gamma factor has to be 7 > lO'^^ for 
the Cherenkov condition to be satisfied. Thus Cherenkov ra- 
diation is not likely to occur in todays radiation background. 
In fact, it is well known that the cosmic rays contain non- 
thermal hadrons, of which some are protons, that can reach 
gamma factors 10^^, but larg e r values are improbable due 
to the GZK cut-off fOreisenl 1 19661 IZatsepin and Kuzminl 
As a comparison, we may consider the situation at 



the time of matter-radiation decoupling. Since (^emitted = 
<o'received(2^/2.7)'*, whcre the temperature T is given in Kelvin, 
we have S" ~ 10^^ J/m^ at the time of decoupling (T « 
8000 K), implying 5 ^ 1G~^^. Thus, the limiting value on 
the gamma factor for the Cherenkov condition to be satisfied 
is 7 > 10^'* — 10^^, still out of reach for high energy cos- 
mic rays. However, as we demonstrate below, the situation 
changes drastically for earlier processes at even higher T. In 
particular, we will focus on the era with 10^ K < T < lO^^K 
when the required 7-factors range from 7 ~ 10^ to 7 > 3. 

The effect presented above is naturally compared with in- 
verse Compton scattering. Setting Z — 1, the cross-section 
for this scattering is u sa irrlml/M'^^, where the classical 
electron radius and M is the charged particle mass. We thus 
obtain a collision frequency v = c,y/a, where ^ is the num- 
ber density of the photons. Comparing this frequency with 
the frequency Vch — {jMc)^^dU/dt, we note that fast par- 
ticles are mainly scattered due to the Cherenkov effect when 
V < i/ch, i-e. 



1 < 
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(122) 



Here T is the temperature of the photon gas, J\f = 
[30C(3)a/fcB7r'*]r3, £ = aT^, ks is the Boltzmann con- 
stant, a = n^k%/15h^c^ « 7.6 x IQ-i^ J/m^K^ and Tch = 
(2025^(3) /UTT^a)meC^/kB « 10^^ K using the polarization 
averaged effective action charge A = {8k + 14k)/2 = 11k. 
Thus, for a single fast proton to be scattered mainly due to 
the Cherenkov effect, we need T > T^h x 10"^ - 10^ K, 
well within the limit of validity of the theory for photon- 
photon scattering. We note that at radiation gas temperatures 
around 10^^ K the quantum vacuum becomes truly nonlinear, 
and higher order QED effects must be taken into account. 

For the early universe considered above, a moderately rel- 
ativistic plasma is also present, which means that collective 
charged particle interactions can play a role. We take these 
plasma effects into account by introducing the plasma fre- 
quency ujp. The photon dispersion relation is « k^c^{l — 
S)+ojp. Thus, the Cherenkov condition is satisfied for charged 

particles with relativistic factors 7 > l/y'j — uj^/k'^c^. For 

the temperatures where the Cherenkov radiation starts to dom- 
inate over inverse Compton scattering, T ^ 10^ — 10^" K, 
we have ujp ^ 10^^^^^ rad/s, and thus Cherenkov radia- 
tion is emitted in a broad band starting in the UV range, 
u) ^ 10^^ rad/s, and continuing up to the Compton frequency 
- 8 X 10^0 rad/s. 

The Cherenkov radiation emitted during the era when T ^ 
10^ K will be redshifted due to the cosmological expansion. 
Thus, the present value of the cut-off frequency will be ap- 
proximately 2 X 10^^ rad/s, i.e. in the short wavelength range 
of the microwave spectrum. However, we do not expect direct 
detection of this radiation in the present universe, since the 
process is only expected to be of importance long before the 
time of radiation decoupling. Still, there are possible impor- 
tant observational implications due to the Cherenkov mecha- 
nism presented here. As shown by the inequality (I122> . the 
effect will be more pronounced for massive particles with a 
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given gamma factor, and protons are therefore expected to 
be more constrained than electrons by the QED Cherenkov 
emission, hi particular, J122> puts stronger limits than Comp- 
ton scattering for supra-thermal protons observed today to be 
relics of the early universe. In fact, it seems rather unlikely, 
given the inequality (I122> . that such protons could survive 
during the T = 10^ - 10^*' K era. 



b. Unmagnetized plasmas Pair production and pair plasmas 
play an important role in th e dynamics of the environments 
surrounding pulsars (see, e.g. lArendt an d Eilek ( 2002'): lAsseol 
( [2003) : Beskin et aL(.1993l) '). Charged particles will attain rel- 
ativistic energies close to the pulsar magnetic poles and radi- 
ate 7-ray photons. This, together with the super-strong mag- 
netic field present around these objec ts (Beskin et al. . 19931) . 
is beheved to produce a pair plasma dTsai and Erbed 11973) . 
Thus, nonlinear QED effects are already known to be an im- 
portant ingredient for pulsar physics. Since the pair plasma 
gives rise to radio wave emissions, and because of the large 
energy scales involved, pulsar atmospheres are likely to host 
other QED effects as well, such as vacuum nonlinearities in 
the form of photon-photon scattering. 

As presented by Stenflo et al.. (.2005.) . for circularly polar- 
ized electromagnetic waves propagating in a cold multicom- 
ponent plasma rather than in vacuum, the wave operator on 
the left-hand sides of Eqs. (|5) and dlOt is replaced by 



1 _2 
□ ^ v^ + -^ 



+fc^ (123) 



where the sum is over particle species j, we have assumed 
that the EM-fields vary as exp{ikz — iujt), the relativistic fac- 
tor of each particle species is 7j = (1 + Qj Eq /mjC^uj^y/^ , 
where Eo denot e s the absolute value of the ele c tric fi eld am- 
plitude dStenflol Il976t IStenflo and Tsintsad ze. 1979). Due 
to the symmetry of the circularly polarized EM waves, most 
plasma nonlinearities cancel, and the above substitution holds 
for arbitrary wave amplitudes. Here tjpj = {nojOj /eonijy^^ 
is the plasma frequency of particle species j and noj denotes 
the particle density in the laboratory frame. 

Next, we investigate the regime o;^ <C k^c^. From Fara- 
day's law and the above inequality we note that the dominat- 
ing QED contribution to Eq. ( I10> comes from the term pro- 
portional to S^B. Combining Eqs. (I10> and (I12> . noting that 
= Bq = k^E'^juP' is constant for circularly polarized 
EM waves, and using uP' <^ k^(?, i.e. M w AneQC^B^B and 
|M| ti>|P|/fc, we obtain from ilOi the nonlinear dispersion 
relation 
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This low-frequency mode makes the particle motion ultra- 
relativistic even for rather modest wave amplitudes. For 
electrons and positrons in ultra-relativistic motion (7^ ^ 
1) with equal densities no and elementary charge ±e, we 
thus use the approximation J2j^pj/lj ~ leriQCOJ / e^E^ = 
2ujp[ijj/ue){E„ii/ Eq) (see (^ and (0), where Up = 



(e^no/eome)^/^. The dispersion relation ( I124> then reduces 
to 
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{Eo/Ecm)ikcuJe/2ujLOp)kc 
(125) 

We note that the ratio uje/(^p is much larger than unity for 
virtually all plasmas, i.e. for electron densities up to ^ 10"^^ 
m^'^. In some applications, such as in pulsar astrophysics, 
it is convenient to re-express the dispersion relation in terms 
of the relativistic gamma factor using Eq/E^^x ~ {'^/^e)l- 
Thus, we obtain 
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from (I125t for the wavelength A — 2'K/k. 
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c. Magne tized plasmas Followin g iMarklund et alJ (12005 dll 
(see also iMarklund et alJ (l2004bll '). for a circularly polarized 
wave Eq = i?o(x ± iy) exp(iA:z — iuit) propagating along a 
constant magnetic field Bo = BqT,, the electromagnetic in- 
variants satisfy 
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Thus, Eq. ( I43t can be written as 



= 0. 
(127) 
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in the Lorentz gauge, and □ = dad"^. For circularly polar- 
ized electromagnetic waves propagating in a magnetized cold 
multicomponent plasma, the four current can be 'absorbed' 
in the wave operator on the left-hand side by the replacement 
(as in the previous section) □ —D{uj, k), where D is the 
plasma dispersion fun ction, given by (see, e.g. IStenflol(il97^ : 
IStenflo and Tsintsadza a979.) ) 
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Here the sum is over the plasma particle species j, ujcj — 
qjBo/rrij and uipj = {nojq^ / e^rfij)^/'^ is the gyrofrequency 

and plasma frequency, respectively, and 7^ = (1 + v'^Y^'^ is 
the gamma factor of species j, with Vj satisfying 
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Here n^j denotes particle density in the laboratory frame and 
rrij particle rest mass. 

The dispersion relation, obtained from Eq. ( I128> . reads 
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We note that as the plasma density goes to zero, the effect due 
to photon-photon scattering, as given by the right-hand side 
of Eq. ( I13U . vanishes, since then llP' — k'^c^ = 

Next, we focus on low-frequency {oj <C kc) mode propaga- 
tion in an ultra-relativistic electron-positron plasma (7e ^ 1), 
where the two species have the same number density uq. 
Then, Eq. ( I131> gives 
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(132) 

For background magnetic field strengths Bo in the pulsar 
range ^ 10'' - 10^° T, cBq < i?ciit, and we therefore drop the 
term proportional to Bq in Eq. ( 1132k Next, using the normal- 
ized quantities = ujLJe/uJp, K — (4a/457r)^-^/'^/ccu;e/^^p 
and E = (4a/457r)i?o/-Ecrit, the dispersion relation ( I132> 
reads 
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The dispersion relation ( I133> describe three different modes, 
two with + polarization and one with — polarization. We note 
that for K <^ 1, the disp ersion relation ( I133t agrees with 
that of iStenflo and Tsints adze ( 1979), whereas in the opposite 
limit K 1, the QED term in ( I133t dominates. For the given 
density, the latter regime applies, except for extremely long 
wavelengths (> 10^ m), and thus we note that QED effects 
are highly relevant for the propagation of these modes in the 
pulsar environment. For small K there is only one mode, but 
two new modes appear for K > 2.6. Thus for large K, appli- 
cable in the pulsar environment, there are three low-frequency 
modes (oj ^ kc) that depend on nonhnear QED effects for 
their existence. 

The effects of the quantum vacuum on electromagnetic 
wave dispersion also allows for nonlinear effects, such as 
wave steepening, shock front formation, and soliton propa- 
gation ( Marklund et al., 2005e) . 



d. Magnetohydrodynamic plasmas When analysing low fre- 
quency magnetised plasma phenomena, magnetohydrody- 
namics (MHD) gives an accurate and computationally eco- 
nomical description. Specifically, a simple plasma model is 
obtained if the characteristic MHD time scale is much longer 
than both the plasma oscillation and plasma particle colli- 
sion time scales, and the characteristic MHD length scale is 
much longer than the plasma Debye length and the gyroradius. 
These assumptions will make it possible to describe a two- 
component plasma in terms of a one-fluid description. The 
one-fluid description means a tremendous computational sim- 
plification, especially for complicated geometries. Moreover, 
if the mean fluid velocity, the mean particle velocity, and the 
Alfven speed are much smaller than the speed of light in vac- 
uum, the description becomes non-relativistic and simplifies 
furth er. 

He vl and HernauistI d 19991) considered the propagation of 
MHD modes, including the effects of photon-photon scat- 
tering and an axion field. Following iThompson and BlaesI 



( 11998b . Heyl and Hernquist start with the Lagrangian [see Eq. 

^ = i^QED + \ataQ'S ifo + + iaeo6'^, (134) 

where the field invariant is defined by (0, and 9 is the axion 
field, that acts as a Lagrange multiplier for the MHD condition 
= 0. The modified Maxwell's equations can be derived 
from Eq. ( I134t . They become [cf. Eq. ( I15H 
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Given a background magnetic field, these equations allow 
for both fast and Alfven modes. The fast modes will suf- 
fer the same type of shock wave formation as presented by 
Hevl and Hernquist ( 1997b), in the absence of the MHD ef- 
fects. A single Alfven mode will not experience the ef- 
fects of photon-photon scattering due to the absence of 
self-interactions. This is not true for the case of counter- 
propagating Alfven modes for which photon-photon scatter- 
ing introduces higher order corrections to their propagation. 



IV. APPLICATIONS 

A. Measuring photon-photon scattering 

Classically, electromagnetic waves only interact indi- 
rectly, via scattering, by passing through a suitable 
medium such as a nonlinea r optical fibre dHasegawaL 1 19751: 
iKivshar and Agrawalll2003T) . To some extent, this is still true 
in QED. One may view the quantum vacuum as a medium 
through which photons scatter off virtual charged particles, 
predominantly electron-positron pairs, producing nonlinear 
effects similar to the ones found in nonlinear optics. How- 
ever, since the nonlinear effects enter the effective Lagrangian 
through the Lorentz invariants, a plane wave will not self- 
interact, and more sophisticated techniques are needed in or- 
der to excite the nonlinear quantum vacuum. In this section, 
such means will be reviewed with the aim of establishing 
methods for direct detection of low energy elastic real photon- 
photon scattering. 

The concept of elastic photon-photon scattering is the- 
oretically well-established. Furthermore, the scattering of 
virtual photons is routinely observed in particle accelera- 
tor environments, and is thus well confirmed in experi- 
ments. Moreover, inelastic photon-photon scattering is also 
experimentally well-confirmed, but this is not the case for 
elastic photon-photon scattering (although experiments have 
been made where it in principle would have been possi- 
ble to make modification such that a direct measurement 
of elastic pho t on-ph oton scattering could have been made 
dBamber et all fl99^ ). Thus, as a fundamental test of QED 
and its predictions about the properties of the quantum vac- 
uum, an experiment on the latter type of scattering may be 
considered an important issue. 
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(a) (b) 



FIG. 26 Feynman box diagrams for (a) Delbriick scattering and (b) 
photon splitting, respectively. Here each cross denotes external field 
legs, e.g. an atomic Coulomb field or a strong background magnetic 
field. 



Closely related to photo n-photon scattering is Delbriick 
scattering ( Delbriick', 1933) and photon splitting ( I Adler etall 
[l970; Adler, 1971; Chistvakov et al., 1998), see Fig. |26| 
(cf. Fig. |4]for a comparison with photon-photon scatter- 
ing). Delbriick scattering is the elastic scattering of pho- 
tons in a Coloumb field, e.g. an atomic nucleus, mediated 
by virtual electron-positron pairs, while photon splitting is 
the down conversion of a photon into two photons of lower 
frequency through an external field, e.g. a strong magnetic 
field. These processes contain external fields mediating the 
interaction between the photons, making the cross section 
larger than for pure photon-photon scattering. In fact, us- 
ing high-Z atomic targets, Delbriick scattering for high en- 
ergy photons has been detected (Jarlskog et al., 1973|), and 
photon splitting, although not detected in a laboratory en- 
vironment, is assumed to be prominent component of many 
astrophysical environments, such as magnetars and soft 7- 
ray re peaters ( Adler and Shubert'. ' 19961 FBaring and Harding! 
I1997I i200k .Harding et al.. 1997). In fact, the splitting of 
photons in the atomic Coulomb field has been reported by 
lAkhmadal iev et al.i (2002), where good agreement with the 
calculated exact Coulomb field cross-section was obtained. 

The effects of photon scattering also manifest them- 
self in the anomalous magnetic moments of the e l ectron 
and the muon (.Bailev et al.. 197 9. Berestetskii et all Il982t 
ICalmet et all ll977HRodionovtl2004 . Even so. the detection 
of direct light-by-light scattering of real photons remains elu- 
sive, even though considerable efforts have been made in this 
area. The possibility to detect low energy photon-photon scat- 
tering would open up for new tests of QED, since fermion loop 
diagrams could give gauge invariant tests of the fermion prop- 
agator, as well discerning between QED and other theories 
predicting or postulating photon properties. Thus, photon- 
photon scattering can both produce interesting effects, as de- 
scribed in previous sections, as well as produce important tests 
for fundamental physical theories. 



1 . Pair production in external fields 

The case of inelastic photon-photon scattering deserves 
some attention in this context. Here the aim is, to some extent. 




FIG. 27 Feynman diagram for the pair creation process. 



anti-matter production on a large scale (see Fig. I27> . There 
are a number of ways, both experimentally confirmed as well 
as schemes suggeste d on numerical or t heoretical grounds, to 
produce and store (Oshima et al., 2004) positronium and an- 
timatter, e.g. laser generated relativistic superthermal elec - 
trons interacting with high-Z materials dLiang et all Il998h . 
the trident process in conjunction with ultra-intense short 
laser pulses in plasmas ( Berezhiani et al., 1992), pair pro- 
ductio n by circularly polarized w aves in plasmas dBulanov 
12004'), laser-thin foil interactions dShen and Mev er-ter-Vehnl 
l200 1b: Shen and Yu, 2003), using Bose-Einst ein co nden- 
sation traps (Greaves et all Il994t TSurko et all Il989^ (see 
ISurko a nd Greaves (2004) for an overview), and using 
fuUerenes (Oohara and Hatakevama, 2003). The formation of 
anti-plasmas and long lifetime trapping of antimatter is cur- 
rently intensely studied, and could shed light on the funda- 
mental laws of nature, e.g. giving new CPT and Lorentz 
invariance tests (Bluhmetal., 1999; Bluhm, 2004), or pro- 
ducing an annihilation laser (Mills, 2002). Since electron- 
positron pairs also constitute a unique type of plasma, promi- 
nent in e.g. the pulsar magnetosphere, the formation of large 
collections of pairs in the laborato r^^will fu rther enable the 
study of astrophysical conditions ( (Greaves and Surkot [1993), 
which we so far have only been able to observe over astro- 
nomical distances, and without control over the physical pa- 
rameter range. The anti-matter production in most laboratory 
applications rely on the plasma being cold. However, as laser 
powers approach the Schwinger critical regime, we will see 
an increased interest in using these for producing high tem- 
perature pair plasmas as well, and for exciting the quantum 
vacuum. 

After the laser was introduced, it was realized that fu- 
ture laser systems could be tools for fundamental physics 
research, an d the parr creation process was reconsidere d by 
Reissl(ll 962^. 'Nikishov and Ritus ('l964a'b', '1965', '1967"), and 
Narozhn v et al. ( 1965). Thus, the mechanism behind the pro- 
duction of electron-positron pairs from electromagnetic fields 
and photons is well-known, and was first directly observed 
by Burke et al. ( 1997) at the SLAC faciUty. Since this ob- 
servation, schemes making use of the next generation laser 
systems has therefore been in the focus of research efforts by 
using the immense laser intensities for producing, not neces- 
sarily cold, pair plasmas in the laboratory. Indeed, as reported 
bv Gahn et al. 12000) . femtosecond table-top lasers can indi- 
rectly create positrons due to electron acceleration in plasma 
channels. 

Since the pair production from the nonlinear quan- 
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turn vacuum formally depe nds crucially on the invariant 
I E p — I B p being positive tSchwingeiil 1 95 ll) . schemes with 
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(by the same argument, strong static magnetic fields does 
not excite the quantum vacuum, unless perturbed). The pair 
pro duction rate per un it volume at the one-loop level is given 
bv ( ISchwingeilll95lh 
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2 oo ^ 
n=l 
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for a uniform electric field E. Here the sum is over the 
real poles in the imaginary part of the integral \\2>\ . Thus, 
the pair creation rate is vanishingly small in most circum- 
stances. The electron-nucleus electric field (although not uni- 
form) requires a nucleus charge of the order a^^ for vacuum 
breakdown, and such nuc lei are unlikely to exist in an y other 
state than a transient one ([Greiner et all[l98 5: Milonnil l 19941 
iReinha rdt and Greiner, '1977^. However, the situation may be 
different for l aser fields, where ultra-shor t high intensity fields 
are available. iBrezin and Itzvksonl ( ll97(]l) derived the pair cre- 
ation rate for varying fields and generalized the pair creation 
rate ( I136t . In the low frequency limit (i.e. uj ^ uje, where LOe 
is the Compton frequency), their expression coincides with 
Eq. ( I136> . taking into account only the first term in the sum. 
Thus, Eq. ( I136t can be used, with good accuracy, to predict 
the pair production efficiency of different processes, even if 
the fields are alternating. 

In all the cases above, the derivations of the pair creation 
rate rely on the assumption of an electric field dominating 
over the magnetic field. In plasmas, the phase velocity v 
can ex ceed the velocity of light. This was used by Bulanov 
12004") to analyze pair production in the field of a circularly 
polarized electromagnetic wave in an underdense plasma. 
Since for a circularly polarized wave c|B| = (fcc/w)|E| = 
(c/t;)|E|, we see that jEp - c^jEp > Q. Thus, the condi- 
tion for pair creation according to Schwinger (1951) is satis- 
fied, and positrons are therefore predi cted to be produced in a 
laser-plasma environment. Moreover. I Avetisvan et alJ ( Il99ll) 
solved the Dirac equation perturbatively to find the produc- 
tion of electron-positron pairs by inelastic multi-photon scat- 
tering in a plasma. They found the probability distribution for 
transverse electromagnetic pertur bations in the plasma, and 
used this ( Avetissian et al., '2002^) to investigate pair produc- 
tion due to nonlinear photon-photon scattering from oppo- 
sitely directed laser beams. Analytical results for the num- 
ber of particles created on short interaction time scales were 
found. Fried et al. (2001) investigated the possibility for pair 
production via crossing laser beams, and concluded that laser 
intensities has to reach 10^^ W/m'^ before this could be used 
as a means for electron-positron generation. 

Pair production may possibly also be achieved without 
the intervention of a plasma or other dispersive media. Ac- 



cording to ' Narozhnv et aD('2004a'l^, focused and/or counter- 
propagating laser pulses can interact via the nonlinear quan- 
tum vacuum as to produce rea l electro n-positron pairs. The 
prediction of N arozhnv et alJ (f2004b") is that pair creation 
for colliding pulses is expected for intensities of the order 
10^6 W/cm^, which is two orders of magnitude lower than for 
single pulse generation. Moreover, INarozhnv et alJ ( l2004ah 
claim that the effect of pair creation puts an upper theoreti- 
cal limit on laser focusing, since the electromagnetic energy 
will be dissipated into fermionic degrees of freedom for high 
enough intensities. 

As intense fields create electron-positron pairs, the parti- 
cle density increases. If intense photon beams can be sus- 
tained for long enough times, this will create a pair plasma. 
In this case, the effects of this plasma on the electromag- 
netic field need to be taken into account. The back-reaction 
of pa ir crea tion on the electromagnetic field was considered 
bv iKluger et alJ ( 1199 1) in 1 + 1 dimensions. Starting from 
a semi-classical approximation, a kinetic model taking pair 
production into account using an emissive term in the elec- 
tron equation of motion was presented. From a numerical 
analysis of the governing equations it was found that high 
enough field intensities will induce plasma oscillations. Due 
to the realization that the right conditions for pair creation by 
lasers could soon be at our disposal, the problem of back- 
reaction and the dynamics of the interaction of the electron- 
positron plasma on the photons has produced an increas- 
ing number of publications over the years. AUdiofer et aU 
(120011) . iProzorkevich et al . (2000), and Rob erts et alJ 
have similarly developed self-consistent schemes where a col- 
lisionless plasma is coupled to the time-dependent electric 
field, via Maxwell's equations and the pair creation source 
term. In their application to the X-ray free electron laser, they 
arrived at plasma behaviour reminiscent of a modulational in- 
stability, and suggested necessary and sufficient conditions to 
generates a pair plasma using the XFEL. Collisions in the 
plasmas created due to intense electromagnetic fields may also 
be taken into account using a quantum kinetic description with 
a pair creation source term ( Bloch_et al., 1999, 2000). 

A so mewhat different sc heme using intense lasers was sug- 
gest bv iLiang et alJ ( Il998ft . Letting two intense laser pulses 
impinge on the surface of a thin foil made of a suitable mate- 
rial, e.g. gold, plasma formatio n takes place. The jitter energy 
for a large fraction 50% ( Wilks et all Il992l) 1 of the pro- 
duced plasma electrons is suggested to exceed the pair cre- 
ation threshold 2meC^. Thus, in this scheme the pair creation 
is a result of the thermal plasma, instead of direct laser i nterac- 
tion w ith the quantum vacuum. Similarly, Helander an d Ward! 
( 2003) suggested that runaway electrons in tokamak plasmas 
could have the same effect. Since electrons with sufficient 
energy experience a decreasing plasma friction force as the 
energy increases, such particles will in effect be accelerated 
to very high energies until direct collisions with plasma parti- 
cles occur. The typical runaway electron energy is > StoeC^, 
and these collisions could therefore trigger positron produc- 
tion, as the electrons loose their energy via brehmsstrahlung 
in the Coulomb field, the so called Bethe-Heitler process. The 
number of positrons in a facility such as JET was estimated to 
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'--^ 10 — 10 , a very large number compared to other labo- 
ratory positron production methods. 

The predicted pair production rates normally assume spa- 
tially uniform electromagnetic fields, which is often in good 
agreement with experimental parameters. However, recently 
oriented crystals have become an important tool in study- 
ing effects of quantum electrodynamics in strong fields, such 
as spin effects in el ectron energy lo s s and crystal assisted 
pair production (see iKirsebom et alJ ( 1200 ll) and references 
therein). In these experiments, the fields may not be consid- 
ered uniform, and the models described above can ther efore 
only p artially account for the observed effects. iNitta et al.l 
(l2004l) remedied this shortcoming by using the trial trajectory 
method jKhokon ov and Nitta, 2002), based on the method de- 
veloped by iBaier and Katkovl ( 119681) . Previous attempts to 
analyze the experimenta l resul ts were based on numerical 
schemes, but iNitta et all ( l2004) found an analytical expres- 
sion for the pair creation rate in a inhomogeneous field, in 
good agreement with the observed pair creation rate. Further- 
more, the inhomogeneous case also displayed pair creation for 
low amplitude fields, where the uniform field treatment effec- 
tively gives a zero pair creation rate. This result could be of 
interest in the case of strongly magnetized stars, which have a 
characteristic dipole behavior. 

Since waves in vacuum are described in terms of their be- 
havior along the null coordinates u = z — ct and v — z + ct, 
it is of interest to generalize the Schwinger results of pair 
creation to t he case of field s depen ding on u and/or v. This 
was done by 'Tomaras et al. (2000), and later generalized by 
lAvan et aL (2003) to a more complicated coordinate depen- 
dence. Furthermore, the momentum spectrum of the produced 
pairs was derived for arbitrary time dependent gauge fields by 
[Dietrich, C2003) . via the exact solution of the equation of mo- 
tion for the Dirac Green's function. 



2. Laser induced pair creation 

The production of anti-matter is of great importance for a 
variety of experimental tests of fundamental issues in physics, 
e.g. Lorentz invariance tests, as well as being of interest in its 
own right. Moreover, there can also be a test of the nonlinear 
properties of QED, since high energy photons may create mat- 
ter and anti-matter out of the quantum vacuum. This process 
is well established as a model for pair production in the vicin- 
ity of neutron stars, and corresponds to the imaginary part of 
the full Heisenberg-Euler Lagrangian, and can thus be inter- 
preted as energy being dissipated from bosonic to fermionic 
degrees of freedom. 

The implications of pair creation was understood very early 
on in the history of QED, but the direct creation of electron- 
positron plasmas from photons has long escaped experimental 
efforts. Thus, an important piece in our view of the quan- 
tum vacuum had long eluded the attempts of detection. How- 
ever, with the rapid advances in laser intensity, the prospects 
for performing a successful experiment in pair creation using 
laser sources took a turn for the better As described below, 
inelastic photon-photon scattering, where two real photons 




FIG. 28 Nonlinear Compton scattering, as given in (7371 . 
(Reprinted with permission from lBamber et alJ fl999,) .) 



gives rise to a real electro n-positron pair, has now been exper - 
imentally confirmed dBamber et allll999l:lBu'rke et al.Lll997^ ■ 
and holds the promise of further elucidating our picture of the 
nonlinear quantum vacuum (see also Meverhofer ( 1997)). 

In nonlinear Compton scattering, multi-photon absorption 
by an electron results in the emission of a single high-energy 
photon according to (see Fig. |28j 

e + ntj^e' + 7. (137) 

The effect ( f737l was first measured by Bula et al ]H996), us- 
ing a GeV electron beam and a terawatt laser source, ob- 
tained by chirped-pulse amplification. In the experiment, 
up to four laser photons interacted with a single electron. 
The high-energy photons produced by nonlinear Compton 
scattering can be used in the laser assisted production of 
a pair plasma. The usage of laser produced photons for 
the electron-positron pair production was sugg ested long be- 
fore lasers reached the ne cessary intensities ( Narozhnv et aQ, 
ll965HNikishov and Rituslll964a,b.. 1965.. 1967.: .Reiss>. 1962) 
(the direct production of pairs by photons requires hiu > 
2mgC^ in the center-of-mass system). By re-colliding the 
high-frequency photons with the original laser photons, ac- 
cordirig to the Breit-Wheeler ^ process (iBethe and Heitled 
ll934lBreit and Wheeled fl93l 

7 + (138) 

the production of electron-positron pairs can be achieved in a 
laboratory environment. This can be compared to the trident 
process 

e + nw^e'e+e". (139) 

While the multi-photon process (I138 |l require s n > 4 with ex- 
perimental values used bv lBurke et alJ (Il997h . the trident pro- 
cess requires n > 5 with the same experiment al data. The 
two-step process ( I137t and ( I138t was used by iBurke et'all 



^ iBreit and Wheeled <1934) considered the single photon scattering uji + 
UJ2 — > e+e~, thus somewhat different from the multi-photon process dis- 
cussed here. 
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( Il997h in the first reported laser production of electron- 

p ositron pairs . 

iBula etalJ(ll996h reported on the observation of the effect 
of nonhnear Compton scattering ( I137> . where the scattered 
electrons were detected using a 46.6 GcV electron beam in 
conjunction with a 1054 nm and a 527 nm laser with focal in- 
tensity ^ 10^* W/cm^. This process can also be understood 
in terms of a plane wave interaction with an electron. For a 
weak electromagnetic field with amplitude E, the maximum 
speed attained by an electron (initially at rest) due to the pass- 
ing of a plane wave is 



eE 



(140) 



where is the rest mass of the electron and uj is the fre- 
quency of the plane wave. As the field strength increases, 
higher order radiation effects becomes important as u^ax ^ c, 
which in terms of light quanta can be interpreted as multi- 
photon absorption by the electron, with the release of a single 
distinguishable light quanta as a result, i.e. the process jl37t . 
In this sense, nonlinear Compton scattering becomes impor- 
tant as the parameter 
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approaches unity. Here, the four-vector potential A'^ satisfies 
the Lorentz gauge. 

For an electron with initial energy Eq, the absorption of n 
photons of the frequency uj at an angle 9 between the electron 
and laser beam, results in the minimum electron energy 
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where s — 2£qLlj{1 + cos 9) is the scattering parameter and 
m-eff = m(l + rfy/^ gives the effective mass. With the ex- 
perimental parameters used by Bula et al. ( 1996), the inten- 
sity parameter becomes rj ^ 0.6. Linear Compton scattering 
(rj 1, n — 1) would then result in w 25.6 GcV at 
9 = 17°. Since the spectrum of multi-photon Compton scat- 
tered electrons extends belo w 25.6 GeV, it w as possible to 
identify the nonlinear effects dBula et allll99d) . 

In the same way, as the intensity parameter rj approaching 
unity signifies the onset of the nonlinear Compton effect, the 
pai-ameter ( Bamber et aU.1999: Burke et aL.1997.) 
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characterizes the strength of the vacuum polarization, as it 
contains both information of the photon frequency as well 
as the background field strength, the two important parame- 
ters for vacuum breakdown. Here, Fab is the Maxwell ten- 
sor of the background electromagnetic field, and pa is the 
four-momentum of the probe photon. As T approaches unity, 
the pair production rate according to the p rocess ( I138> be- 
comes significa nt (Burke etal., 1997; Naro zhnv et all 1 19651 
iNikishov and Ritu^ .1964a.k .1965> .1967.) . For the case of 
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FIG. 29 Calculated positron energy spectra for a 30 GeV photon 
interacting with a 527 nm laser beam. In panel (a) the polarization is 
parallel while in panel (b) the polarization is perpendicular, n gives 
the number of p hotons involved in th e interaction. (Reprinted with 
permission from lBamber et alJ ^99^-) 



single particle (n = 1) Breit-Wheeler scattering, laser wave- 
lengths of 527 nm would require single photon energies of 
111 GeV in order for significant pair production to occur, 
while for the multi-photon Breit-Wheeler process the photons 
of the same wavelength colliding with backs cattered photon s 
with energies 29 GcV gives T « 0.5?7 (Burke et al.Lll997l) . 
Thus, for large enough 77, the pair production rate would yield 
a detectable level of electrons and positrons (see Figs. 1301 and 
13 1 >. with a we l l-defin ed energy spectrum (Fig. I29l l. 

burke et'aP ( Il997h for the first time presented the results 
of a successful measurement along the lines presented above. 
The signal consisted of ^ 100 positrons above the background 
value using a 46.6 GcV electron beam an d a 527 nm Nd:glas s 
laser with focal intensity ^ 10^* W/cm^ dMeverhofeJ 1 1 997^ . 



3. Other mechanisms for pair production 

'Narozhn v et al ] (|2004a") considered pair production in an 
electromagnetic field created by two counter propagating laser 
pulses, and showed that pair production can be experimen- 
tally observed when the intensity of each beam is similar to 
10^6 W/cm^, three orders of magnitude lower than that of a 
single pulse. However, the cross-section for the Schwinger 
process at optical frequencies (or below) is so small at any 
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FIG. 30 Positron production rate per Compton scatterer as a func- 
tion of the intensity parameter r;, as given by <141> . The solid line 
is the numerical estimate from the two-step process <137> and <138> . 
while the dashed line re presents the trident process <139t . The mea- 
surements performed bv lBurke et alJ (1997) are given by the dots in 
the plot. (Reprinted with permission from Burke et al. (1997).) 
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FIG. 3 1 Pair production rate, as compared between the multi-photon 
Breit-Wheeler process <138> and the trident process (El, as a 
function of T given by jl43> . (Reprinted with permission from 
iBamberet al]<1999h .') 



laser intensity that this effect is insignificant dMittelmanl 
[1987). 

Production of pairs is also possible in the Coulomb field 
of a nucleus via virtual photons ("tridents"), which is a 
dominant energy loss mechanism at high energies. In a 
trident Bahba process, high energy electrons, with kinetic 
energies exceeding the pair production threshold 2meC^, 
can produce electron-positron pairs by scattering in the 
Coulomb potential of the nucleus . In the past , some au- 
thors ( Bunkin and Kasakoy, lilQt IScharer et all Il973l) had 
presented a preliminary discussion about pair production by 
relativistic electrons accelerated by intense laser, while oth- 
ers (Berezhiani et al., 1992) presented a detailed investigation 
of pair production due to scattering of high energy electrons 
produced in strong wake fields driven by intense short laser 
pulses. This was found to be an efficien t mec hanism for a 
"pair factory". Recentlv. iBerezhiani et alJ (12005) carried out 
computer simulations of laser plasma dynamics in overdense 



plasmas and showed that an intensive production of pairs by 
the drive motion of plasma electron takes place d ue to the tri- 
dent process. Furthermore, Bulanov et alJ (l2005l) have shown 
that electromagnetic waves could be damped due to electron- 
positro n pair production (see also Mikheev and Chistvakovl 
(l200lh for a discussion on the process in a strong magnetic 
field). 



4. Laser experiments on photon-photon scattering 



T he evolution of 
ing dMo urou et all 



iTaiima and Mouroa 



laser intensity is truly astound- 
19911 12005 : Perrv and Mourou. 1994 
20021) (see Fig. [l), and with the event 



of the X-ray free electron laser, a new domain in experimen- 
tal physics will open up. There have been an interesting set 
of both suggested and performed experiments using lasers of 
previous and current intensities. Note that one of the major 
obstacles in these investigations have been residual gas com- 
ponents in the vacuum environment. However, depending on 
the problem of study, the means for inhibiting the residual 
gas to have a detrimental effect on the measurement varies. 
In high intensity laser experiments on elastic photon-photon 
scattering, the electron expulsion at the leading edge of the 
laser pulses will in fact make the generation of background 
radiation weaker (at a vacuum of 10^^ torr), and particle ef- 
fects would therefore ha ve a negligible effect in these experi- 
ments (iLundstrom et all 12005). This is contrast to weak field 
experiment, such as cavity environments, where the effects 
due to residual gas may be significant. However, it is possible 
to design the mode interaction such as to produce a unique 
signature of photon-photon interaction, thus making it possi- 
ble, in prin ciple, to detect the scat tering by the proper filtering 
techniques dEriksson et alll2004l) . 



a. Vacuum birefringence The concept of vacuum birefrin- 
gence is well know n and has bee n theoretically explored in 
many publications dAdlej. 1 197 it jAdler and Shubert, 199e" 
| l966t iHevlandHernauisll UJ 



ErbeJ. 



and Hernauisd 11997a: .Klein and Niga 



1964albl) . The birefrin gence of the vacuum manifests itself as 



the difference in the refractive index between the propagating 
ordinary and extraordinary modes ( Rikken and Rizzo. .20000 . 
Thus, although a very difficult high precision experiment, 
this difference may in principle become measurable in strong 
enough background magnetic (or electric) fields. This idea 
has been exploited in the PVLAS set-up IBakalov et allll994 
iMelis sinoi. l2002h . for which the difference dBakalovetal 
1998) 
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is to be measured. Here |Bo| is given in Tesla. 

A linearly polarized laser beam is sent through the static 
field Bo, with E • Bq = |E||Bo| cos 6*. Due to the birefrin- 
gence of the magnetized vacuum, as given by An, the beam 
will attain an ellipticity 
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over a propagation distance L, where A is the wave length 
of the radiation. The change in ellipticity is prop osed as a 
measu rement of the birefringence of vacuum. Bak alov et al.l 
41998) also presented a detailed discussion of noise sources as 
well as a rather detailed description of the actual experimen- 
tal setup. Current superconducting magnets can reach field 
strength up to 5 — 25 T, and could in principle yield detectable 
changes in the polarization state of a laser beam traversing 
it. Unfortunately, the strong magnetic fields generate forces 
within the detection equipment which may interfere with the 
ellipsometric measurement. Moreover, magnetic fields cannot 
be shielded in any efficient way, and this is there fore a prob- 
lem that is likely to persist d Cameron et al.Ln"993l) . 

Another approach towards measuring vacuum birefrin- 
gence is by using laser interferometry. This technique 
can reach astonishing accuracy and sensitivity, and is cur- 
rently the most prom ising me t hod o f choice in gravita- 
tional wave detection dSaulsonl Il994h . Using laser inter- 
ferometry for detecting light-by-light scattering through vac- 
uum birefringence rests on the same principle as described 
above, but replacing the strong magnetic field by ultra-short 
laser pulses (Boer and van Holten, 2002: Luiten and PeterseiiL 
l2004btlPartovilll993l) . Since laser beams, i.e. laser light with 
typical pulse length much larger than its wavelength, has a 
very low energy density compared to the strongest laboratory 
magnetic fields, one has instead to resort to ultra-short highly 
focused laser pulses. Such configurations could indeed result 
in magnetic field components of the order 10^ T, orders of 
magnitude larger than quasi-stati onary magnetic fields pro - 
duced by superconducting coils dLee and Fairbanks! l2002h . 
Due to the degree of focusing of the pulse, interaction of the 
strong field with the detector can be almost eliminated. On 
the other hand, the ultra-short time- and length scales require a 
very high resolution in the detection. An experimental sugges- 
tion along these lines was put forward bv'Luiten and PetersenI 
f2004 ab), usin g a high precision biref ringence measuring 
technique ( Hal l et allboOOl) . iLuiten and Petersen ( 2004b) ar- 
gue that this technique may be used to construct a table-top 
detector of vacuum polarization using current state-of-the-art 
optical techniques. The set up consists of two concentric res- 
onant cavities with an interaction cross-section. One of the 
cavities acts as the vacuum polarizer, while the other cavity 
supplies the test photons for which the ellipticity is to be de- 
tected (see Fig. I32t . Depending on the Fabry-Perot resonator 
reflectivity, the integration time was estimated. With a reflec- 
tivity R = 99.97% the necessary operation time of the de- 
vice would be 2.6 years, R ~ 99.994 % yields 1.7 days, and 
R = 99.997% gives 2.5 h, using a 20W 200 fs laser and a 
resonator of length 3 m. 

IJeah-Sheng et all ll2004 have built and tested a 3.5 m 
Fabry-Perot inteferometer with a precision ellipsometer for 
QED tests and axion search, along the lines of the vacuum 
birefringence te st presen ted above. Note that the results pre- 
sented by Jea h-Sheng et al. (2004) are for a prototype detec- 
tor, and, although promising, a measurement of the vacuum 
polariz ation has not been performed (see also Sheng-Jui et al. 
( l2003h V Moreover, the PVLAS collaboration has recently 
claimed dZavattini et all 12005.) to have measured the dichro- 
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FIG. 32 The interferometric set up for detection of vacuum bire- 
fringe nce. The figure reference in the above set up is figure 1 
in .Luiten and PetersenI t2004h) . (Reprinted with permission from 
iLuiten and PetersenI <2004d) .') 



ism of a magnetized vacuum, an effect possible due to inter- 
action between light pseudoscalars and photons. 



b. Harmonic generation iDing and KaplanI lll989l) suggested 
that the nonlinear vacuum could be given measurable prop- 
erties by the possible harmonic generation of radiation in 
an external field. The work attracted lots of attention, and 
some questi ons as to whether the re s ult was correct or not 
was raised ("Ding and Kaolai^, Il990t iFord and Steel Il990t 
iRaizen and Rosenstein, 1990). While some of the critiques 
were flawed, the main problem in the work of Ding and 
Kaplan was the assumption of a constant background field 
(Ford and Steel, 1990). It is well-known that a homogeneous 
and time-independent background field cannot transfer mo- 
mentum to the photons, and such a field is therefore not ca- 
pable of driving a frequency upshift as sugg ested by Ding 
and Ka plan (see McKenna and Platzman ( 1963)). However, 
iDing and Kaplan ( 199 2) showed that a spatially inhomoge- 
neous background field could indeed result in higher harmon- 
ics. This idea was further developed by Kapla n and Ding ! 
(UOOO), where Maxwell's equations were analyzed with a 
weakly varying background magnetic field. 

Using the slowly varying amplitude approximation 
(|Hasegawa, 1975), Kaplan and Ding (2000) showed that the 
envelope of the electromagnetic field satisfies the second har- 
monic generation equation 
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where 
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is the second harmonic generation background source term. 
Here D^^) and H^^) are derived from the Heisenberg-Euler 
Lagrangian (|5} from the nonlinear field combinations giv- 
ing rise to terms proportional to cxp{2iujt). Kaplan and Din^ 
(2000) use Eq. ( I146I I to study the evolution of 2-D Gaussian 
beams propagating in an external non-constant magnetic field, 
giving the estimated output power. Moreover, a discussion of 
more complicated background magnetic field geometries, e.g. 
the magnetic quadrupole case, was given. Considering /iin 
lasers with focal intensities ^ 10^^ W/cm^ generating a pulse 
propagating through the backg round magneti c field strengths 
lO'' T [which in the paper bv .Kaplan and Dingi (k200Qi) is sug- 
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gested to be produced by explosive mechanisms], a rough es- 
timate gives a production of 85 photons/day by second har- 
monic generation. However, in this estimate, temporal effects, 
which may be of importance in the next generation ultra-short 
intense lasers, have been omitted, and could yield alterations 
in their estimates. 



c. Four-wave interactions In the second harmonic generation 
presented above, the interaction of photons is mediated by a 
background magnetic field. However, crossing electromag- 
netic waves would similarly interact and yield new modes of 
different frequencies. One of the more prominent modes in 
such a mechanism is given by the four-wave interaction me- 
diated mode satisfying resonance condition between the fre- 
quencies and w avevectors (i.e. photon energy and momentum 
conservation) (,Rozanovlll993ll) . It is therefore not a surprise, 
given the evolution of laser powers and frequencies, that the 
search for photon-photon scattering using resonant four-wave 
interactions has caught the attention of researchers in this area. 
This approach has also come furthest in the experi mental at- 
tempts to detect elastic scattering among photons ernard", 
[1298, 1999, 2000; Bernard et al., 2000; MouHn and Bernard. 

iMoulin et alJ ( Il996h presented experiments on light-by- 
light scattering performed in the optical regime. With this, 
they managed to put new experimental upper limits on the 
photon-photon scattering cross-section. Unfortunately, no 
scattering was detected, but stimulated the continued research 
along the lines of four-wave interactions as an experimental 
tool for probing the quantum vacuum. Using the resonance 
conditions** uj4 = uji + uj2 ^ (^3 and k4 = ki + k2 — ka, 
one may, in general, derive a set of wave interaction equations 
for slowly varying amplitudes a,- , i = 1, . . . , 4, of the form 
dWeiland and WiIhelms'sonLll977l) 

^ = Ca.akal (148) 

given any type of media through which the waves may inter- 
act. Here the coupling constants C depend on the the interac- 
tion in question, as well as on the physical parameters of the 
system around which the waves are modulated. In the case 
of a nonlinear quantum vacuum, the coupling constant will 
depend on k of the Lagrangian (|5}- 

The coupling constants may be interp reted in terms of the 
nonlinear susceptibility of the vacuum. iMoulin and Bernard 
(1999) considered the interaction of three crossing waves, 
characterized by their respective electric field vectors E^, pro- 
ducing a fourth wave E4. Starting from Maxwell's equations 
with the usual weak field limit Heisenberg-Euler third order 
nonlinear corrections (see the Lagrangian (|5}), they derive the 



equation 

(149) 

for the driven wave amplitude E4, where the overall harmonic 
time dependence exp(— icji) has been factored out. Here x'^' 
is the third order nonlinear susceptibility given by 

^(3) ^ ^ JL « 3 X 10-41 X if mVv2, (150) 

where K is a dimensionless form factor of order unity. The 
value of K depends on the polarization and propagation direc- 
tions of the pump modes, and reaches a maximum of K = 14 
for degener ate four-wave mixing ( Moulin and Bernard, 1999). 
iMouhn and Bernard ( 1999) furthermore discuss the influence 
of a non-perfect vacuum, where the susceptibility of the gas 
will introduce a threshold, in terms of a critical g as pressure, 
for the nonlinear QED effect to be detected. [Bernard et alJ 
( 12000) and Bernard (2CQQ) recently presented experiments on 
four-wave mixing in vacuum, improving previous attempts 
by nine orders of magnitude, although no direct detection of 
photon-photon scattering was achieved. Experiments along 
the same lines as described for four-wave mixing above can 
also be used for a large number of other, non-QED, test, 
such as axio n*^ search (Bernard, 1999; Bradlev et al., 2003 
iDupavs et al.Lf2005l) . Thus, progress of low-energy QED ex- 
periments could also prove to be useful for, e.g. dark matter 
searches. 

There are more recent proposals for detection of 
photon-photon scattering using four-wave interactions. 
Lundstromet al. (2005|) has done more detailed calculations 
concerning experimental constraints, in particular for the As- 
tra Gemini laser (o perational in 2007) at the Rutherford Ap- 
pleton Laboratory dCCLRCl ^005), as well as non-perfect 
vacuum problems etc., and concluded that it will be feasible 
to detect elastic scattering among photons if using a high rep- 
etition rate high intensity laser system. 

5. Cavity experiments 

As mentioned in Sec. lIII.A?! the effects of photon-photon 
scattering on cavity EM fields is to produce new wave modes. 
The new modes excited in the cavity will (approximately) sat- 
isfy the cavity dispersion relation. Thus, by varying the cavity 
cross-section, the pump modes may be filtered out, leaving the 
new modes for detection. The treatment of cavit y mode inter- 
action in the qua ntum vacuum \yas des cribed bv lBrodin et alJ 
(I200l[l2002h and lEriksson etaP (l2004 . 

If no damping or dissipation is present, Eqs. ( I69I I yield a 
linear growth of the vector potential amplitude A3 of mode 



The interaction between modes of different frequencies gives rise to several 
new modes, but the resonance conditions and time averaging, mimicking 
the act of detection over certain timescales, yield the desired equations. 



' Axions are bosons which were introduced in order t o explain the absence of 
CP symmetry breaking in QCD 1 Peccei and Ouinil lTgTTtlWeinbereLI 1 97il 
IWilczekLll978l) . and the axion is still to be detected. 
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3. In order to gain an understanding of the saturation level, 
we make the following modification to Eq. (I69> . Let d/dt 
d/dt — {ijJj,/2ttQ), where Q is the cavity quaUty factor A 
steady state amplitude 



90^ i?e'„. 



(151) 



is thus obtained. Here 073 is the frequency of the mode gen- 
erated by the third order QED nonlinearities. The number 
of excited photons in the cavity mode can be described by 
N ^ {eq J d^r)/hw'i,. Using the saturation value M5l\ 
for the vector potential of mode 3, the number of photons 
generated by the nonlinear interaction of two cavity modes 
is given by the expression 
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where the c oupling constant K^y\ can be found in 
lEriksson et al.l (I2004) . V = na^Zo, a is the cylindrical cav- 
ity radius, zq is the cavity length, and /Ja is a zero for the 
Bessel function Ji corresponding the the generated mode sat- 
isfying the resonance condition i6li . We note that the num- 
ber of generated photons depends on a large number of pa- 
rameters, and one thus needs to specify the cavity geometry 
etc. i n order to obtain an es timate of the magnitude of the ef- 
fects. lEriksson et al.l (l2004l) choose the following wave mode 
numbers: (£1, ^2, ^3) = (3, 15, 21) (fulfilling £3 ^ 2£i + I2), 
^2=^3= 3.83, corresponding to the first zero of Ji, and 
/?i = 7.01 corresponding to the second zero. This gives 
zo/a — 9.53 through the frequency matching condition i6l\ 
and determines the frequency relations to a'3/w2 = 1-26 and 
oJ^li jJA = 1.12. Su bstituting these values gives K^^yx = 3.39 
(see lEriksson et al.l ( i2004.) ). The remaining key parameters 
are th e quahty factor and the pump field strength. iLiepd 
has shown that it is possible to reach intense cavity 
surface fields, of the order IA2I ^ 0.01 - 0.03 Vs/m, 
with quality factors as high as Q = 4 x 10^" at temperatures 
of 1 K. Thus, in this case Eq. M52\ gives 
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For a cavity wall temperature of 0.5 K, the number of thermal 
photons is iVthemai ~ \ / [ex.\){hjjj, / ksT) — 1] sa 7, where ks 
is the Boltzmann constant, which is thus lower than iVgED- In 
order to reach further accuracy in the measurement, a cavity 
filtering system can be set up, so that the pump modes may 
be reduced or eliminated. It can furthermore be shown that 
the nonlinearities of the cavity wa lls will no t generate modes 
swamping the QED photons (.Eriksson et alii2004) . 



B. Laser-plasma systems and the X-ray free electron 
laser 

The X-ray free electron laser promises new and exciting 
applications for a coherent electromagnetic source. The ap- 
plications range from probing astrophysical conditions in the 



laboratory to new possibilities to do molecular biology. Could 
the XFEL also provide insight into quantum electrodynamical 
nonlinearities, such as photon-photon scattering? If affirma- 
tive, this would enhance our understanding of the quantum 
vacuum, as well as providing new prospects of testing fun- 
damental properties of physical laws, such a Lorentz invari- 
ance and symmetry breaking. Indeed, it has been stated that 
the facilities at DESY and SLAC would be able to produce 
electron-positr on pairs directly from the vacuum ( Rin gwaldl 
l20QlalbLl2()()l . due to the estimated focal intensities at these 
sources. If this scenario is demonstrated, it is likely that 
the effects of elastic photon-photon scattering would come 
into play at an even earlier stage. Due to the possible ef- 
fects of scattering amongst photons, such as photonic self- 
compression and collapse, it is therefore of interest to include 
such effects into the analytical and numerical models used in 
predicting the behaviour of these systems. Furthermore, the 
creation of a pair plasma in the laboratory could be affected 
by new low-frequency modes from nonlinear quantum vac- 
uum effects, thus altering the properties of energy transfer 
within such plasmas, as well as providing indirect tests for 
QED. XFEL will also gives the opportunity to do laboratory 
astrophysics in a new parameter regime, making the quantum 
vacuum more accessible. 

However, it is not necessary to enter the new regime of 
XFEL in order to facilitate tests of QED and Lorentz invari- 
ance, as well as doing laboratory astrophysics. Such effects 
as U nruh radia t ion ( ^ Unruh, 1976) and the related Hawking ef- 
fect (iHawnEi can hopefully b e investigated usi ng the 
next generation laser-plasma systems fBinghaml |2003), such 
as the high repetition- rate A stra Gemini laser (to be opera- 
tional 2007) fcCLRd I2OO5I) . In such regimes, it will also 
be of interest to investigate QED effects, such as photon- 
photon scattering. As seen in the previous sections, the in- 
troduction of plasma dispersion allows for new electromag- 
netic wave modes in both unmagnetized and magnetized plas- 
mas (Marklund et al., 2004b, 2005d), when nonlinear quan- 
tum vacuum effects are included. For example, at the labo- 
ratory level, current high laser powers are able to accelerate 
particles to highly relativistic speeds. Furthermore, pulse self- 
compression in laser-plasma systems may play an important 
role in at taining pow er levels well above current laser limits 
(see, e.g. iBulanovet a l. (2003); Sh orokhov et alJ ( l2003l) ). As 
the intensities approach the Schwinger limit in future laser- 
plasma setups, effects of pair-creation and photon-photon 
scattering have to be taken into account ( Bulanov et alil20()3t 
iBulano v( l2004 . Laser-plasma systems can have electron den- 
sities of the order 10^^ m~ ^, and laser intensities can be 
close to 10^^ - 10^^ W/cm^ ('Bingham'. 2003": 'Mourou et all 
.1998) . Moreover, as stated by Bulanov et al. ( 2 003) , laser 
self-focusing in plasmas could come close to E^nt, at which 
pair creation is likely to follow. In fact, it has been estimated 
that the National Ignition Facility would be able t o produce 
pairs by direct irradiation of a deuterium pellet dLee et all 
Il995h . On the other hand, the creation of laboratory electron- 
positron plasmas is alr eady a feasible task ( Greaves et aQ 
1994.' Surkoetallll989() . as is the usage of these p lasmas for 
making pair plasma experiments fcGreaves and Surka 1 199i. 
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Thus, the possibiUty to study laser-plasma interactions in pair 
plasmas could be a reality in the nearby future. The currently 
available positron densities in the laboratories are well be- 
low th ose of regular laser plasma systems, but according to 
iLee et al. (1995) there is a possibility of reaching densities of 
order 10^^ m~^. Using no ^ 10^^ m^'^, and the field inten- 
sity 10^^ V/m (due to laser self-compression ('Bulanov et al[ 
|2P03)) at the wavelength 0.3 x 10"*^ m, we find from Eq. ([125} 
that uj w 7.8 X 10^ rad/s, i.e. the frequency is in the LF band, 
as an example of QED effects in laboratory plasmas. 

In combination with plasma particle expulsion due to elec- 
tromagnetic wave trapping, the possibility of cat astrophic col- 
lapse due to p hoton-p hoton collisions arises (.Bulanov et all 
|2003; Marklun detail E004al This scenario it highly inter- 
esting, since it would make 3-dimensional ultra-intense soli- 
tonic structures possible in vacuum bounded by a plasma or 
wave-guide structure, a truly exciting prospects. This may 
even prove a valuable tool for intense pulse storage, if a suc- 
cessful cavity nonlinear QED experiment is performed. 



C. Astrophysical importance 

The implications of the QED vacuum is well-known 
within astrophysics, and the pair plasma in pulsar sur- 
roundings is partly dependen t on mechan isms which 
has no classical counterpart (Asseo, "2003"). Further- 
more, photons splitting ( Adler et al., 1 970; Adler, 1971; 
iBialvnicka-Birula and Bialvnicki-Birulal Il970l) supports the 
notion that strongly magnetized objects, such as neutron stars 
and magnetars, could be used as probes of nonlinear QED ef- 
fects. 

However, most of the effects discussed within astrophysi- 
cal applications concerning QED deals with single-photon ef- 
fects, and thus do not take collective effects into account. It 
is well-known from plasma physics that collective effects al- 
ter the charged particle behaviour in non-trivial and important 
ways. In fact, it would not be possible to understand most 
plasma effects without resorting to a collective description. 
The analogy betwe en the quantum and a p lasma system has 
been stated before d Dittrich and Giesl bOOOh. and is both use- 
ful and imaginative. Thus, in line with this, it is likely that col- 
lective quantum vacuum effects could yield crucial informa- 
tion about astrophysical systems, where extreme energy levels 
are reached. Even kinetic effects, such as Landau damping, 
could play a role in the dynamics of photons in the vicinity of 
strongly magnetized objects. This could prove a new realm of 
photon kinetics, and the a pplications to astrop hysical sources, 
such as magnetar quakes ('K ondratvevtl2002h . are of interest 
for future research directions. 

Especially strong magnetic field effects due to the quantum 
vacuum is of great interest in astrophysical applications. Since 
the earth-based magnetic field strengths are very limited, and 
are likely to remain so for the foreseeable future, magnetars 
and similar objects offers a unique perspective on the quan- 
tum vacuum (Barin2 and Hardina, 2001; Erber, 1966). Pulsar 
magnetospheres exhibit extreme field strengths in a highly en- 
ergetic pair plasma. Ordinary neutron stars have surface mag- 



netic field strengths of the order of 10^ — 10^ T, wh ile mag- 
netars can reach 10^" — 10^^ T dKouveliotouL 1 1 998i) . coming 
close to, or even surpassing, energy densities eoE^^^ corre- 
sponding to the Schwinger limit. Such strong fields will make 
the vacuum fully nonlinear, due to the excitation of virtual 
pairs. Photon splitting c an therefore play a s i gnificant role in 
these extreme systems dBaring and Harding 1200 it iHardingt 
119911) . 

Neutron stars have surface magnetic field strengths of the 
order of 10*^-10^ T , while magnetars can reach lO^^'-lO" T 
( lKouveliotoulll998h . coming close to energy densities corre- 
sponding the Schwinger critical limit eoi^crit; here, the quan- 
tum vacuum becomes fully nonlinear. Single-particle QED 
effects, such as photon splitting can play a significant role 
in the understanding an d interpretation of observations from 
these extreme systems (f Baring and Harding 1200 it iHardingl 
Il99ll) . In fact, the pair plasma creation in pulsar environ- 
ments itself rests on nonlinear QED vacuum effects. The 
emission of short wavelength photons due to the accelera- 
tion of plasma particles close to the polar caps results in a 
production of electrons and positrons as the photons propa- 
gate through the pulsar intense magnetic field ( Beskin et all 
Il993h . The precise density of the pair plasma created in this 
fashion is difficult the estimate, and the answer is model de- 
pendent. However, given the Goldreich-Julian density ncj = 
7 X 1015(0.1 s/P)(B/10®T) m-3, where Pis the pulsar pe- 
riod and B the pulsar magnetic field, the pair plasma den- 
sity is e xpected to satisf y nn = Mng j . M b eing the mul- 
tipUcity jBeskin et all Il993t iLuo et all l2002l) . The multi- 
plicity is determined by the model through which the pair 
plasma is assumed to be created, but a moderate estimate is 
M = 10 dLuo et all 12002). Thus, with these pre-requisites, 
the density in a hot dense pair plasma is of the order 10^^ 
m^"^, and the pair plasma e xperiences a relativistic factor 
^ 10^ — 10^ ('Asseql l2003l) . We may use these estimates 
to obtain estimates for particular QED processes in plasmas. 
For example, inserting the above values in ( I126t . we obtain 
A ^ 10^1^ — 10^11 m. On the other hand, the primary beam 
will ha ve n p noj and 7 10^ - 10^ (Asseo, 2003), at 
which ( I126t yields A ~ 10^* — 10^^ m. Thus, in this case we 
obtain short wavelength effects. 

The field of laboratory astrophysics ties the experimental 
domain of laser-plasma systems to areas of r esearch where w e 
so far have been restricted to observations (l HEDLAil200l . 
Interesting studies, such as shock front formation relevant to 
supernova explosions, could in principle be achieved in fa- 
cilities such as NIF. However, the scales of the astrophysical 
event and the laboratory setup differs by orders of magnitude. 
Thus, it is reasonable to ask if it is possible to apply laboratory 
findings to astrophysical events. .R yutov e t aL (2000^ con- 
sider the prospects of investigating MHD phenomena of rele- 
vance for supernova hydrodynamics. From self-similarity in 
the governing system of equations and boundary conditions, 
as well as from dimensionless variables (such as the magnetic 
Reynolds number Rcm) they argued that the laboratory re- 
sults could be translated to astrophysical settings (however, 
Rcm in the laboratory cannot reach the extreme values of su- 
pernova ejecta but can reach values much larger than 1). Sim- 



43 



ilarly, 'Budil et al." ('200(/) discussed the applicability of peta- 
watt lasers to radiative-hydrodynamics relevant to, e.g. super- 
nova rem nant evolutio n. The testbed experimental results pre- 
sented by Budil et al. ( 2000) indicated that the results could be 
useful in calibrate models of radiation hydrodyna mics in su- 
pernova remnants (see also lShigemori et al]ll200fll) \ Thus, the 
use of high intensity lasers for probing astrophysical phenom- 
ena, in particular as tool for testing and calibrating simulations 
of certain events, has undergone rapid development over the 
last decade. For testing QED effects within astrophysical sys- 
tems the relevant dimensionless parameters are the frequency 
compared to the Compton frequency, the field strength over 
the Schwinger critical field strength Q, as well as the sign of 
the relativistic invariant c^B^ — E^. As can be seen by the 
second of these requirements, the laboratory experiments of 
today will at most be weakly nonlinear, whereas the effects 
in astrophysical systems, such as magnetars, can be strongly 
nonlinear However, the combined effect of laser-plasma dy- 
namics and vacuum nonlinearities would yield unique signa- 
tures, and could be probes of more exotic phenomena in astro- 
physical plasmas. One such example in the testing of the Un- 
ruh effect (Chen and Taiima, 1999; Unruh, 1976) as a means 
of understanding the Hakwing effect (iHawkind. 119741) . 



the formation of ultra-high intensity pulse trains, due to self- 
compression and pulse splitting, is a truly exiting prospect. 
The fact that such configurations are within laboratory reach, 
using the next generation laser-plasma facilities, makes the 
predicted effects and their connection to astrophysical events 
even more interesting and may open up new possibilities for 
basic and applied research in the future. 
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V. CONCLUSION AND OUTLOOK 

The possibility of simulating astrophysical events in a lab- 
oratory enwonm^ntj^^ durin g the last decade, progressed 
dCheni I2OO3I iRemingtonl Eo0 5). Apart from the astrophys- 
ical tests, laser-plasma systems also provide an opportunity 
to test certain aspects of fundamental physics, e.g. the prop- 
erties if the quantum vacuum, via strong fields. Strong (~ 
10 — lOOMV/m) coherent electromagnetic fields can nowa- 
days be produced in superconducting cavities (Graber, 1993), 
and fields within plasmas could come close to the Schwinger 
limit ([1). Moreover, QED effects are part of many astrophys- 
ical phenomena, such as pair cascading, and thus laboratory 
astrophysics has a natural connection to investigations of the 
quantum vacuum. 

Here, we have reviewed the implications of QED correc- 
tions to classical electrodynamics and the propagation of elec- 
tromagnetic waves and pulses. In particular, QED corrections 
on photon-plasma interactions were described. The modifica- 
tions introduced by the nonlinear quantum vacuum were con- 
sidered for, e.g. coherent and incoherent pulse propagation. 
Analytical, perturbative, and numerical ways of analyzing the 
governing equations were presented. Moreover, the properties 
of nonlinear collective effects were presented, such as three- 
dimensional pulses collapse and the formation of light bullets. 

The application of the results can be seen both in an astro- 
physical context as well as in a laboratory setting. For ex- 
ample, in magnetar environments (Kouveliotou, 1998) pho- 
ton splitting ( Adler. Il97ll) is important and it is believed to 
give a plausible explanation for the radio silence of magnetars 
([Harding, 1991 ). On the other hand, collective effects, such as 
the ones presented here, could give valuable insight of QED 
phenomena in astrophysical environments. In the laboratory. 
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